Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 24.4% intergenic (‑188/+120) 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,705,9780AC24.4% 469.1 / 145.0 235intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (74/87);  new base C (44/2);  total (146/89)
Fisher's exact test for biased strand distribution p-value = 1.33e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-02

Only 100 of 240 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077
                                                                                                   |                                                                                                          
 aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                           >  1:85964/1‑100 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                       <  2:3606580/100‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                     <  1:564723/99‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                    >  1:2170085/1‑100 (MQ=255)
         ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                   >  1:1994346/1‑100 (MQ=255)
          tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                  <  1:3196072/100‑1 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                 <  2:830180/100‑1 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                 >  2:590897/1‑100 (MQ=255)
           gacGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                 >  2:2579282/3‑100 (MQ=255)
            ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                <  1:3371274/100‑1 (MQ=255)
             cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                               <  2:2024779/100‑1 (MQ=255)
             cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                               <  1:2364407/100‑1 (MQ=37)
             cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                               <  2:2198857/100‑1 (MQ=37)
              ggATGGCGCTGCACGCAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                              <  1:367745/100‑1 (MQ=25)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                              <  1:1830443/100‑1 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                              <  2:654093/100‑1 (MQ=37)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                             <  1:1439088/100‑2 (MQ=37)
                   gcgcTGCACGCAGAGGCACAGCGTCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                        <  1:511058/100‑1 (MQ=255)
                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATCCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                        <  1:2011271/100‑1 (MQ=255)
                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                      >  1:935600/1‑100 (MQ=255)
                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                      >  2:3026780/1‑100 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                >  1:2084578/1‑100 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTg                                                                                >  2:3473030/1‑100 (MQ=21)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                               <  2:1254493/100‑1 (MQ=255)
                               gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑caa                                                                           >  1:798010/1‑100 (MQ=255)
                                 ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAAtt                                                                         <  2:1885534/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                                       >  1:2202194/1‑100 (MQ=255)
                                   cacaGCGCCACCGGGTATAATAATTGTTGCCAAAGTGCGCATTGTTTATGCCGGATGCGGCGTTAACGCCCTATCCGGCCTAC‑AAAATCGTGccattt                                                                          >  1:2654461/1‑93 (MQ=17)
                                     caGCGTCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAAt                                                                     <  2:1266067/100‑1 (MQ=255)
                                          ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACAtt                                                                <  2:2182920/100‑1 (MQ=255)
                                           cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTg                                                               >  2:3005161/1‑100 (MQ=255)
                                            aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGc                                                              <  2:247163/100‑1 (MQ=255)
                                            aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGc                                                              <  2:1647400/100‑1 (MQ=255)
                                            aCCGGGTATACACCTTGTTGACAAAGTCCGCTTTGTTTATCTCGGATGCGTCCTTAACGCCTTATCCGTCCTAC‑AAAATCGTGC‑AAATTCATTACATTGc                                                              <  1:313699/100‑1 (MQ=17)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑g                                                           <  1:2032926/100‑1 (MQ=255)
                                               gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑Ga                                                          >  2:1992425/1‑100 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                               <  1:2800359/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                               >  1:3224683/1‑100 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                               >  1:1271249/1‑100 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                               <  2:492490/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                               <  1:1283176/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                               >  2:3218815/1‑100 (MQ=255)
                                                          ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAgg                                           >  1:3510846/1‑100 (MQ=255)
                                                          ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAgg                                           <  1:1917468/100‑1 (MQ=255)
                                                           tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGc                                          >  1:3800681/1‑100 (MQ=255)
                                                                aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGAt                                     <  2:1501042/100‑1 (MQ=255)
                                                                    aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                                 <  1:3751427/100‑1 (MQ=255)
                                                                        cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                             <  1:2709829/100‑1 (MQ=255)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGTTAAGCg                                >  2:2423031/1‑97 (MQ=255)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                            >  2:2063265/1‑100 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                           >  2:2643108/1‑100 (MQ=255)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                         >  2:3356334/1‑100 (MQ=255)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                         <  1:838261/100‑1 (MQ=255)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                         <  2:1240099/100‑1 (MQ=255)
                                                                             tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                        >  1:2726716/8‑100 (MQ=18)
                                                                             cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGATCGTGC‑AAATTCAATATATTGCA‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                        >  1:3669143/8‑100 (MQ=17)
                                                                              ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCC‑AAATTCAATATATTGCA‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                       >  2:2307935/7‑100 (MQ=18)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                       >  1:2451771/1‑100 (MQ=255)
                                                                                gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAATTC‑TGC‑AAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  1:1548395/5‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTAT‑AAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                    >  2:2889443/4‑100 (MQ=25)
                                                                                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                  >  2:1530808/1‑100 (MQ=255)
                                                                                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                  >  1:1945101/1‑100 (MQ=255)
                                                                                   ggtcggataaggcgTTCACGCCGCATCCGGCCCAC‑AAAGTCTTGC‑AAAATCAATAAATTGCAGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2767944/11‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGC‑TAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1737619/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1389603/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGC‑AAATTCAATATATTGC‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2116612/2‑100 (MQ=18)
                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGC‑TAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  1:1362251/1‑99 (MQ=18)
                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGC‑AAATTCAATAAATTGCA‑‑A‑AT‑CAACATGTAGGCCGGATAAGCGCAGCGCATCAGGCa                  >  1:2571710/1‑100 (MQ=18)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGC‑TAATTCAATATATTGTG‑GG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCTCATCAGGCag                 >  1:2497197/1‑99 (MQ=18)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGC‑TAATTCAATATATTGCA‑GG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                 >  1:1217899/1‑99 (MQ=18)
                                                                                     ctggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGC‑AAATTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                <  2:2915479/92‑3 (MQ=17)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGC‑AAATTCAATATATTGCA‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  1:3558594/9‑98 (MQ=17)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGT‑AAATTCAATAAATTGAT‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:105779/9‑98 (MQ=18)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTT‑AAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  1:1152417/9‑98 (MQ=21)
                                                                                     ccggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                >  2:1799584/9‑100 (MQ=18)
                                                                                     ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGT‑AAATTCAATAAATTGAT‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  1:3438258/1‑98 (MQ=18)
                                                                                      cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGC‑AAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑CTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg              >  1:3469626/8‑99 (MQ=18)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGC‑AAACTCAATATATTGCA‑GAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt               >  1:1803314/8‑97 (MQ=17)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGC‑AAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc             >  2:1475008/7‑98 (MQ=18)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGC‑AAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc             >  1:1032296/7‑98 (MQ=18)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGC‑AAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc             >  1:1381035/7‑98 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  1:2235449/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTAT‑AAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct              >  1:1859341/7‑96 (MQ=25)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGT‑AAATTCAATAAATTGAT‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  2:557841/7‑100 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTT‑AAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc              >  2:1211845/7‑96 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGA‑AAATTCAATACATTGCA‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct              >  2:732943/7‑96 (MQ=17)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGC‑AAATTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatag              >  1:2915465/7‑96 (MQ=17)
                                                                                        gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagttct             >  1:2310443/6‑95 (MQ=18)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  1:3025232/100‑1 (MQ=255)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  1:2341624/100‑1 (MQ=255)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  2:2232917/100‑1 (MQ=255)
                                                                                          tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc           <  2:1952107/100‑1 (MQ=255)
                                                                                            cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         <  1:3742334/100‑1 (MQ=255)
                                                                                             ggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGC‑AAACTCAATATATTGCA‑GAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagcgtt        >  2:881092/1‑94 (MQ=17)
                                                                                               cgTTTACGCCGCATCCGGCCTAC‑AAGATCGTGC‑AAATTCAACATATTGC‑‑CA‑CT‑CACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCAattttacattt       >  2:2714125/1‑94 (MQ=25)
                                                                                               cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      >  2:3673206/1‑100 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  2:2958997/100‑1 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc     <  1:1401484/100‑1 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      >  2:1704593/1‑99 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcn    <  1:281638/100‑2 (MQ=255)
                                                                                                   |                                                                                                          
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077

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