Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA37.9% 409.5 / 316.6 264intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (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 = 1.86e-01

Only 100 of 277 total aligned reads displayed.

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070
                                                                                                    |                                                                                                            
      tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                        <  1:549120/100‑1 (MQ=255)
         aGTATTGCTGGCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                     <  1:847356/100‑1 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                  <  1:1482907/100‑1 (MQ=255)
               gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                               >  1:3010286/1‑100 (MQ=255)
                ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                              >  2:2197720/1‑100 (MQ=255)
                  gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                            >  2:1721599/1‑100 (MQ=255)
                    cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                          <  2:1545906/100‑1 (MQ=37)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                         <  2:886869/100‑1 (MQ=37)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                        <  2:324655/100‑1 (MQ=255)
                       aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                       <  2:3674621/100‑1 (MQ=255)
                       aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                       <  1:306446/100‑1 (MQ=255)
                       aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                       <  1:1670167/100‑1 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                             <  1:3380800/100‑1 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                             >  1:829657/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                           <  2:960325/100‑1 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                           <  1:2353725/100‑1 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                           >  1:431454/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                           >  1:1864015/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                           >  2:1017670/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                           <  2:74745/100‑1 (MQ=255)
                                   gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                          <  1:2405646/100‑1 (MQ=255)
                                     agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                        <  1:2851927/100‑1 (MQ=255)
                                      tgcgtgcaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                       >  1:1053604/7‑100 (MQ=255)
                                        cgtgcaGCGCCATACGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                      >  2:2403165/5‑100 (MQ=17)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                   >  1:2747808/1‑100 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                   >  2:467330/1‑100 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                   >  1:2468664/1‑100 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                 <  1:1149802/100‑1 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                 <  1:1954244/100‑1 (MQ=255)
                                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGa                                                    <  2:2264495/100‑1 (MQ=255)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACGTTGC‑‑AGAAtttt                                            <  1:2545522/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                        <  1:1017674/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGACGGATGCGGCGTTAACGCCTTATCCGGACTACAAAGTCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                        <  2:3840933/100‑1 (MQ=255)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑Gt                                       >  2:3597139/1‑100 (MQ=255)
                                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAgg                                    <  2:2897720/100‑1 (MQ=255)
                                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTACAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAgg                                    >  2:2447239/1‑100 (MQ=255)
                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGc                                   >  2:659254/1‑100 (MQ=255)
                                                                   cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGc                                   >  1:2617192/6‑100 (MQ=18)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGcc                                  >  1:2056876/1‑100 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGGATTTTC‑‑GTAGGCCt                                 >  2:3735165/1‑99 (MQ=255)
                                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGAt                              <  2:1889246/100‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATa                             <  1:1209518/100‑1 (MQ=255)
                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGc                          <  1:3653156/100‑1 (MQ=255)
                                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCg                         <  2:249966/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                     <  2:3361887/100‑1 (MQ=255)
                                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcgc                   >  1:2189719/1‑100 (MQ=255)
                                                                                   cttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAgcgc                   >  1:3261556/10‑100 (MQ=14)
                                                                                    gtttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  2:1128730/9‑100 (MQ=18)
                                                                                     cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  2:59665/8‑100 (MQ=17)
                                                                                      ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:3033894/7‑100 (MQ=18)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                <  1:2422553/100‑1 (MQ=255)
                                                                                       tttATGCAGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  2:1937997/100‑1 (MQ=255)
                                                                                       tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAGTTCGTGCAAATTCAATAAATTGC‑‑AATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:355080/6‑100 (MQ=18)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  2:1476382/100‑1 (MQ=255)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  1:1381135/100‑1 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:584944/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:2372414/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:1264774/5‑100 (MQ=21)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAACATGCAAATTCAATAGATTGC‑‑AGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:1999707/5‑100 (MQ=17)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:2789865/4‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1724317/4‑100 (MQ=18)
                                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             <  2:2441015/100‑1 (MQ=255)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2124592/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGC‑‑AATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc            >  1:2737539/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCAAATTCAACATATTGC‑‑AATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:218065/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCAtaa               >  1:1441892/3‑95 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1962226/3‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc            >  1:3761027/3‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑‑TCGATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1304816/3‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:3308462/3‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1351576/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:298374/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AGGAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc            >  1:685846/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTAC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc            >  2:2729459/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:3414957/3‑100 (MQ=21)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:1572263/1‑100 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:2714300/1‑100 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:945137/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:471890/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:3232206/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  2:2658378/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1902340/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2517076/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:1994075/2‑100 (MQ=18)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:3505543/1‑100 (MQ=255)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:3370183/1‑100 (MQ=255)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:3632710/1‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:3563343/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:745626/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag          >  1:2640866/1‑99 (MQ=18)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  1:2819595/1‑100 (MQ=255)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  2:3532196/1‑100 (MQ=255)
                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt        >  2:1614247/1‑100 (MQ=18)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATTGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGTATTtttg     <  2:2044357/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:1053600/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  2:1261191/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  1:2321858/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  2:429465/100‑1 (MQ=255)
                                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  1:2862474/1‑100 (MQ=255)
                                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  <  1:958040/100‑1 (MQ=255)
                                                                                                    |                                                                                                            
AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070

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