BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_012967 1,329,516 C→T 100% H33Y (CAC→TAC)  topA → type I DNA topoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129671,329,5160CT100.0% 49.9 / NA 19H33Y (CAC→TAC) topAtype I DNA topoisomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/11);  total (8/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGC  >  NC_012967/1329484‑1329547
                                |                               
gTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                              <  2:536308/36‑1 (MQ=37)
   gTGACTACGTGGTGAAATCCAGCGTCGGGTACATcc                           >  2:638003/1‑36 (MQ=25)
    tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                          <  2:917386/36‑1 (MQ=37)
        tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                      >  1:961844/1‑36 (MQ=37)
        tACGTGGTGAAATCCAGCGTCGGGTACATTCGCGAt                      >  1:770289/1‑36 (MQ=16)
          cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                    <  2:741418/36‑1 (MQ=37)
             ggtgAAATCCAGCGTCGGTTTCATCCGCGGTTTGcc                 >  2:196985/1‑36 (MQ=16)
             ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                 <  1:442916/36‑1 (MQ=38)
              gtgAAATCCAACGTCGGTTACATCCGCGATTTGCCg                <  1:56434/36‑1 (MQ=37)
                 aaaTCCAGCGTCGGTTACAACCGCGATTTGCCGAcc             >  1:904040/1‑36 (MQ=25)
                   aTCCAGCTTCGGTTACATCCGCGATTTGCCGACCAg           <  2:859986/36‑1 (MQ=37)
                   aTCCACCGTCGGTTACATCCGCGATTTGCCGACCAg           <  2:523755/36‑1 (MQ=37)
                     ccAGCGTCGGTTACATCCGCGATTTTTCGACCAGTg         >  2:466689/1‑36 (MQ=19)
                     ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg         <  1:589067/36‑1 (MQ=39)
                     ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg         <  1:502295/36‑1 (MQ=39)
                        gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt      <  2:583720/36‑1 (MQ=39)
                         cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc     >  1:636767/1‑36 (MQ=39)
                         cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc     <  2:39012/36‑1 (MQ=39)
                           tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg   <  1:576345/36‑1 (MQ=39)
                            cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc  >  2:886622/1‑36 (MQ=39)
                                |                               
GTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGC  >  NC_012967/1329484‑1329547

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: