Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 25.9% intergenic (‑258/+50) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0480GT25.9% 325.6 / 100.9 158intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (63/50);  new base T (35/3);  total (103/55)
Fisher's exact test for biased strand distribution p-value = 2.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

Only 100 of 163 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                         |                                                                                                  
 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATCAGCGGa                                                                                                   >  2:1217667/1‑100 (MQ=255)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                  >  1:587145/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                <  1:3040692/100‑1 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 <  1:405163/99‑1 (MQ=255)
     tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                                               >  1:746255/1‑100 (MQ=255)
       tttgtaggCCGGATAAGTCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTACAGAGAT‑‑TGC‑‑GTAGGCCTGATAAGCGTAg                                                                                                  >  2:1462092/8‑95 (MQ=11)
       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                             >  1:2836239/8‑100 (MQ=18)
        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATc                                                                                            >  2:2399942/7‑100 (MQ=18)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCTTGCAAATTCAATACATTGCAGAAAT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                            >  1:1652348/1‑100 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGGATT‑‑TTC‑‑GTAGGCCTGATCAGCGGAGCGCATc                                                                                            <  1:2044564/100‑1 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGAATCa                                                                                           >  1:577411/1‑100 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑ATC‑‑GTAGGCCTGATAAACGGAGCGCATCa                                                                                           >  2:2268461/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  1:1815373/6‑100 (MQ=17)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          >  2:3020255/5‑100 (MQ=18)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          <  1:2399970/96‑1 (MQ=18)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATCCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATAAgg                                                                                         >  2:1289242/1‑100 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑CCT‑‑GTAGGCCGGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2050239/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:85464/2‑100 (MQ=14)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑TAT‑‑GTAGACCTGATAAGCATAGCGCATCAGGCa                                                                                       >  1:1212855/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGAGAC‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                                                       >  1:401355/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATA‑‑CGT‑‑ATAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:321438/2‑100 (MQ=17)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:370551/1‑100 (MQ=255)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACACAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  1:25540/1‑99 (MQ=17)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                                                     >  1:494082/1‑97 (MQ=21)
               ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                      >  1:449408/9‑99 (MQ=21)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  1:146806/9‑98 (MQ=17)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                     >  2:257069/9‑98 (MQ=18)
               ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                     >  1:2028258/1‑100 (MQ=255)
               ccGGATCCGGCGTTAACGCCTTATCCAGCCTACAAAATTGTGCAAAATCAATTGATTGCACAGCT‑‑AGC‑‑GTAGGCCTGATGCGCTACGCTCATCAGGc                                                                                        <  2:1483766/97‑1 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                  >  1:2576818/7‑98 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:131475/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                   >  2:756655/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                   >  2:1708710/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACT‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                   >  2:1239319/7‑96 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                                                  >  1:1066030/7‑98 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCATACAAATCATGCAAATCCAAAATATTGCAGGAGA‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACAttta                                                                                   >  1:2756742/7‑98 (MQ=18)
                   aTGGGGCGTTATCGCTTTATAGGGCGAAAAAAATCGTGGAAATTCACTACATTGCAGAATT‑‑TTC‑‑GGAGGGCTAATAAGTGGAGCGCATCAGGCATTtttg                                                                                 <  2:1974961/100‑1 (MQ=16)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAATTCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:89473/100‑1 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                               >  1:895262/1‑100 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              <  1:2913962/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           >  2:636789/1‑100 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:2555109/98‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         <  2:475307/100‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         >  2:1710077/1‑99 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        <  2:1381428/100‑1 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  1:2420009/100‑1 (MQ=255)
                               cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCa                                                                     >  1:2333888/1‑100 (MQ=255)
                                gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                                    <  1:1116001/100‑1 (MQ=255)
                                 ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCCTTTGCCATCatt                                                                   >  2:3045457/1‑98 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             >  1:1980526/1‑100 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             <  1:2202677/100‑1 (MQ=255)
                                          gCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑TTAGGCCTGATAAGCGTAGCGCATCATGCATTTTTGCCTTTGCCATCAGTCtt                                                                <  2:2656694/94‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  2:1831134/100‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  2:1626755/100‑1 (MQ=255)
                                               cAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtaga                                                     >  1:2102142/1‑97 (MQ=255)
                                                   aTCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                                                 >  2:919909/1‑100 (MQ=255)
                                                    tCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                                                <  2:388500/100‑1 (MQ=255)
                                                     cGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                               >  1:1187900/1‑100 (MQ=255)
                                                      gTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                              >  2:2082608/1‑100 (MQ=255)
                                                      gTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                              >  2:2680400/1‑100 (MQ=255)
                                                       tGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCTCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:997315/100‑1 (MQ=255)
                                                       tGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             >  1:2595544/1‑100 (MQ=255)
                                                       tGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  2:2148835/100‑1 (MQ=255)
                                                       tGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:1549065/100‑1 (MQ=255)
                                                        gCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             >  2:1093481/1‑99 (MQ=255)
                                                        gCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            >  2:1334593/1‑100 (MQ=255)
                                                        gCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            >  2:2894226/1‑100 (MQ=255)
                                                        gCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  1:1476580/100‑1 (MQ=255)
                                                          aaaTTCAATACATTGCAGAAGT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAGCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                          <  2:2741684/100‑1 (MQ=255)
                                                             ttCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                       <  2:1690613/100‑1 (MQ=255)
                                                              tCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTATATGGCTTTTAATTATTCGATagag                                      >  1:632627/1‑100 (MQ=255)
                                                                 aTACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg                                   >  2:1666298/1‑100 (MQ=255)
                                                                    cATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                >  1:68886/1‑100 (MQ=255)
                                                                     aTTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                               >  2:1113891/1‑100 (MQ=255)
                                                                        gCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  1:259708/100‑1 (MQ=255)
                                                                        gCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  1:260192/100‑1 (MQ=255)
                                                                         cAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  1:276113/1‑100 (MQ=255)
                                                                          aGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  2:2333861/100‑1 (MQ=255)
                                                                          aGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  2:2391597/100‑1 (MQ=255)
                                                                            gAATT‑‑TTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         <  2:68886/100‑1 (MQ=255)
                                                                            gAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         <  2:2233571/100‑1 (MQ=255)
                                                                              aTT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       >  2:2302553/1‑100 (MQ=255)
                                                                              aTT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       >  1:2223029/1‑100 (MQ=255)
                                                                              aTT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       >  2:1879276/1‑100 (MQ=255)
                                                                                     tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   <  2:2521330/100‑1 (MQ=255)
                                                                                           agtccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:429393/97‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:950906/100‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:701614/100‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:468000/100‑1 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  2:1048013/1‑100 (MQ=255)
                                                                                               cTGCTAAGCGCAGCCTATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           <  2:1944670/100‑1 (MQ=39)
                                                                                                tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  2:1132205/1‑99 (MQ=255)
                                                                                                tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          >  1:2870006/1‑100 (MQ=255)
                                                                                                    aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      <  2:284376/100‑1 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  1:2970680/1‑100 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  2:2172312/1‑100 (MQ=255)
                                                                                                       cGGAGCGCATCAGGCATTTTTGCCTTTTCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   >  2:3008719/1‑100 (MQ=255)
                                                                                                        ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  >  2:374413/1‑100 (MQ=255)
                                                                                                        ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  1:1816693/100‑1 (MQ=255)
                                                                                                         gAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  1:1786123/99‑1 (MQ=255)
                                                                                                         |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: