Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9700GA37.9% 336.0 / 282.5 264intergenic (‑180/+128)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (72/92);  new base A (87/13);  total (159/105)
Fisher's exact test for biased strand distribution p-value = 8.35e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-02

Only 100 of 276 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                       
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTAATGCCGGATg                                                                                                         >  2:2545499/1‑100 (MQ=255)
    gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                     <  1:3766443/100‑1 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                   <  2:2194208/100‑1 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGTACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                              >  1:3676795/1‑100 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                            <  1:1482907/100‑1 (MQ=255)
              ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                           <  2:2409021/100‑1 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                        >  2:1854025/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                      >  1:2710255/1‑100 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                   <  2:886869/100‑1 (MQ=37)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                 <  2:3674621/100‑1 (MQ=255)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                 <  1:1670167/100‑1 (MQ=255)
                         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                <  1:968000/100‑1 (MQ=37)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                               >  2:2258523/1‑99 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                              <  2:101957/100‑3 (MQ=25)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                       >  1:829657/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                     <  2:960325/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                     >  1:431454/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                     <  2:2123197/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                     <  2:74745/100‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                    <  1:2405646/100‑1 (MQ=255)
                                       tgcgtgcaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                 >  1:1053604/7‑100 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                               <  1:689852/100‑1 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                             >  2:2490258/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                             >  1:557721/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                             >  2:467330/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                             >  1:2747808/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                             >  1:2468664/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                           <  1:1149802/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                           <  1:1954244/100‑1 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                      <  1:3172799/100‑1 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                    <  1:2843988/100‑1 (MQ=255)
                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t                                           <  2:1674616/100‑1 (MQ=255)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACGTTGCA‑GA‑At‑ttt                                         <  1:2545522/100‑1 (MQ=255)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑ttt                                         >  1:3019903/1‑100 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTc                                        <  2:15635/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCg                                       <  1:1017674/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCg                                       <  1:2506486/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGACGGATGCGGCGTTAACGCCTTATCCGGACTAC‑AAAGTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCg                                       <  2:3840933/100‑1 (MQ=255)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGt                                      >  2:3597139/1‑100 (MQ=255)
                                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAgg                                   >  2:2446772/1‑100 (MQ=255)
                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGc                                  >  2:659254/1‑100 (MQ=255)
                                                                   cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGc                                  >  1:2617192/6‑100 (MQ=18)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GG‑AT‑TTTCGTAGGCCt                                >  2:3735165/1‑99 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATa                            <  1:1209518/100‑1 (MQ=255)
                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGc                         >  2:1945300/1‑100 (MQ=255)
                                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGAtgcgct                        >  1:204003/1‑95 (MQ=255)
                                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGa                      <  1:1854034/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                    <  2:3361887/100‑1 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                   <  1:2945038/100‑1 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                   <  1:949760/100‑1 (MQ=255)
                                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcgc                  >  1:2189719/1‑100 (MQ=255)
                                                                                   cttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAgcgc                  >  1:3261556/10‑100 (MQ=14)
                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCa                 <  2:1095845/100‑1 (MQ=255)
                                                                                     tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCAt                >  2:3265691/8‑100 (MQ=18)
                                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCAt                >  2:2294694/1‑100 (MQ=255)
                                                                                     cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CCACGTAGGCCTGATAAGCGTAGCGCAt                >  2:59665/8‑100 (MQ=17)
                                                                                      ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc               >  2:3033894/7‑100 (MQ=18)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc               >  2:1356618/1‑100 (MQ=255)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc               <  1:2422553/100‑1 (MQ=255)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg             <  1:1381135/100‑1 (MQ=255)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg             <  2:1211273/100‑1 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:3571225/5‑100 (MQ=17)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  1:2372414/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAg             >  1:1264774/5‑100 (MQ=21)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TAT‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  1:1999707/5‑100 (MQ=17)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:2789865/4‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1724317/4‑100 (MQ=18)
                                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg            <  2:2006929/100‑1 (MQ=255)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:2124592/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:601315/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:1982260/3‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑TCTGGTAGGCCTGATAAGCATCGCGCATCAGGc           >  1:2737539/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:1506232/3‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑AT‑A‑‑TGACGTAGGCCTGATAAGCGTAGCGCATCAGGcg          >  1:3399402/3‑99 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CCACGTAGGCCTGATAAGCGTAGCGCAtaa              >  1:1441892/3‑95 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1304816/3‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:3308462/3‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1351576/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:298374/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTACAGGG‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGc           >  2:2729459/3‑100 (MQ=18)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc           >  1:1572263/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:1903404/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  1:234876/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:471890/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:3440470/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:3232206/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa          >  2:2658378/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1902340/2‑100 (MQ=17)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:3505543/1‑100 (MQ=255)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:3370183/1‑100 (MQ=255)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  2:3632710/1‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:3563343/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:745626/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCac         >  1:903681/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGCag         >  1:2640866/1‑99 (MQ=18)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt        >  2:3532196/1‑100 (MQ=255)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt        >  1:2819595/1‑100 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  2:290028/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  2:429465/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  1:2321858/100‑1 (MQ=255)
                                                                                                   |                                                                                                       
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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