Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 30.3% intergenic (‑195/+113) 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,9850TC30.3% 506.2 / 264.0 313intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (109/109);  new base C (86/9);  total (195/118)
Fisher's exact test for biased strand distribution p-value = 7.56e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-01

Only 100 of 315 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                            
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                             >  1:2252936/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                            >  1:3246108/1‑100 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCTTTAACGCCTTacc                                                                                                           <  2:2898226/100‑3 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                           >  2:2335483/1‑100 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGAGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                          >  1:3039823/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                        >  2:2086148/1‑100 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                      <  1:3732519/100‑2 (MQ=37)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                      <  2:3968923/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                     <  2:4047473/100‑1 (MQ=255)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                     <  2:646497/100‑3 (MQ=37)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACa                                                                                               <  2:3414368/100‑1 (MQ=255)
             cgttGCACGCAGAGGTACAGCGCCCCCGGGTATATAGATTGTTGTCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACa                                                                                               <  1:3858042/96‑1 (MQ=21)
             cgcTGTTCGCAGAGGCACAGCGCTACCGGGTATACAGATTGGTGACAAAGTGCGCTTTGTTTATTCCGGATAAGGCGTTATCGCCGCATCCGGCAT‑‑aaac                                                                                               <  1:1002361/100‑4 (MQ=21)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaaa                                                                                             >  1:1268005/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaaa                                                                                             >  2:1877786/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaaa                                                                                             >  2:1919390/1‑100 (MQ=255)
               cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaaaa                                                                                            <  2:3086585/100‑1 (MQ=255)
                tGCAGGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAAt                                                                                           <  1:2319343/100‑1 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAAt                                                                                           <  2:1767500/100‑1 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAAt                                                                                           <  2:1124123/100‑1 (MQ=255)
                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGt                                                                                        >  1:4011915/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                                       <  1:1867245/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                                       <  1:2980139/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                                       <  1:147500/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGc                                                                                      <  2:1994244/100‑1 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAAtt                                                                                 >  1:2217753/1‑100 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAAtt                                                                                 >  2:2262042/1‑100 (MQ=255)
                          gacacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAAtt                                                                                 >  2:3654339/3‑100 (MQ=255)
                           gtgcaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGTTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAAGCGTGCa                                                                                     <  1:3399583/92‑1 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCa                                                                               >  2:4241280/1‑100 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAAt                                                                             <  2:1035893/100‑1 (MQ=255)
                               aGCGCCACCGGGTCTACAGATTGTTGACAAAGTGCGCCTTGTTTATACCAGATTCGGCGTTCACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTTAATa                                                                            >  1:462295/1‑100 (MQ=17)
                                     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGc                                                                     <  2:2502954/100‑1 (MQ=255)
                                             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                          <  2:3308850/100‑1 (MQ=255)
                                              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                                         >  1:577399/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                      <  2:2109352/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                      <  2:3196924/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGAGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                          >  1:407070/1‑97 (MQ=255)
                                                   gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAg                                                   >  2:2848951/7‑99 (MQ=17)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                                <  1:2871173/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                                <  1:644276/100‑1 (MQ=255)
                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAggagat                                              >  2:1686246/1‑96 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTg                                              <  2:1370459/99‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATaa                                          >  1:3012156/1‑98 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                       <  1:1978695/100‑1 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                       >  2:2175300/1‑100 (MQ=255)
                                                               tGCGCGTTGTTTATGCCGGATTCGGAGTTATCGGCTTATCTGGCCT‑‑ACAAAATCGTGCAAACTCAAGACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCgg                                      <  1:2624705/100‑1 (MQ=18)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                                    <  2:3472903/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                   <  1:1911322/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                   >  2:4162666/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                   <  1:1567754/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:1531116/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  2:547377/1‑100 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:2047092/100‑1 (MQ=255)
                                                                    nctttgtaggCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGTCCTGCAGATTCAATATATTGCAGGTACTGC‑‑GTAGGCCTGATAAGCGTAgc                                   >  1:3130058/10‑100 (MQ=9)
                                                                    attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATAAACAAAATCGTGCAAATTCAATAAATTGCAGAATTTCC‑‑GTAGGCATGCTAAGCGGa                                     >  1:3543556/10‑98 (MQ=2)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                <  1:3805391/100‑1 (MQ=255)
                                                                      tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAACATATTGCAATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                               >  2:27686/8‑100 (MQ=18)
                                                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                               >  1:1600225/1‑100 (MQ=255)
                                                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                               >  1:2553322/1‑100 (MQ=255)
                                                                        tgtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCT‑‑ACAAGAACGTGCAAATTCAACAAATTGTGATGCATCCT‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  2:115490/6‑100 (MQ=21)
                                                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  2:2557763/5‑100 (MQ=17)
                                                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCAT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:679164/5‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:251099/4‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:3502878/4‑100 (MQ=21)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:3457958/4‑100 (MQ=21)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAGTCGTGCAAATTCACTACATTGAGGAATTTAC‑‑GTAGGCCTGATAAGGGGAGCGCATCAgg                           >  1:861178/1‑100 (MQ=17)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                           <  1:1429360/100‑1 (MQ=255)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                           >  2:184129/1‑100 (MQ=255)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:384718/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:3253923/3‑100 (MQ=18)
                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          <  1:1191859/100‑1 (MQ=255)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑CTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:2564641/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2640927/2‑100 (MQ=18)
                                                                             gtcggataaggcgTTCACGCCGCATCCGACCC‑‑ACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCaa                        >  1:3884023/10‑99 (MQ=21)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:4083402/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  1:1646114/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:2704278/1‑99 (MQ=18)
                                                                               cggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                      >  1:1786424/8‑100 (MQ=17)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  2:2957040/7‑98 (MQ=21)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:3491711/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1127275/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  2:1312661/7‑96 (MQ=25)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                     >  1:280916/7‑96 (MQ=21)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  2:494128/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  1:1320148/7‑100 (MQ=17)
                                                                                ggataaggcgTTAACGCCGCACCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                     >  2:3754633/7‑100 (MQ=18)
                                                                                 gataacgcgTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCCGAAACCCT‑‑GTAGGCCTCAACAGCGGCGCGGATCACGc                          >  1:2085452/7‑94 (MQ=9)
                                                                                  aTGCGGCGTTAACGTCTTATCCGGCCT‑‑ACAAAATCGGGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:2836547/100‑1 (MQ=255)
                                                                                   taaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAAGCAGTtttgc                  >  2:1971434/4‑100 (MQ=17)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                  >  2:2425261/1‑100 (MQ=255)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                  <  1:813497/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATACGCGGAGCGCATCAGGCATTtttgcct                <  2:2605978/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  1:2889977/100‑1 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            >  2:771379/1‑100 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  2:871169/99‑1 (MQ=255)
                                                                                            accGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt         >  2:2888554/3‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       >  2:1487354/1‑100 (MQ=255)
                                                                                               gCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      >  2:1515926/1‑100 (MQ=255)
                                                                                                 cTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc    <  2:1605103/100‑1 (MQ=255)
                                                                                                   |                                                                                                            
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

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