Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,2540GA32.6% 85.4 / 18.8 49S29L (TCA→TTA) emrKtripartite efflux pump membrane fusion protein EmrK
Reads supporting (aligned to +/- strand):  ref base G (26/7);  new base A (0/16);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 7.32e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTG  >  NC_000913/2483010‑2483504
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           
tttATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACg                                                                                                                                                                                                                                                       >  2:551113/1‑251 (MQ=255)
       gTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaataa                                                                                                                                                                                                                                                 >  2:554873/1‑250 (MQ=255)
          tGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCACTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaataaata                                                                                                                                                                                                                                              >  2:3989/1‑250 (MQ=255)
                aCAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAATATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTAc                                                                                                                                                                                                                                      >  2:498481/1‑251 (MQ=255)
                  aaTATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCg                                                                                                                                                                                                                                    <  2:290611/251‑1 (MQ=255)
                     atTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGcc                                                                                                                                                                                                                                  >  2:144830/1‑250 (MQ=255)
                     atTTGCCAGATTATTTTTAGATTTATTGAGTGCGATAGTGGAATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGCACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATGCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATACCATTGCCCAATAGGCATAGGCACCTGAAAACGCAATACATAAAACTACCGCCa                                                                                                                                                                                                                                 <  2:32726/251‑1 (MQ=255)
                       ntGCCAGATCATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGCAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTACTTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAAt                                                                                                                                                                                                                               >  2:339831/2‑251 (MQ=255)
                                             aTTGAGTGCGATAGTGGCATCACTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTGGTTCGTATTTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAAAAATTGAATCTGCATCCCCCGGGACATAGGCGTCTTCTGTACTAATCATGTCTTCTATTCCCATTGACCAATAGGCATAGCCCCCTAAAACCGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTg                                                                                                                                                                                                          <  1:3989/250‑1 (MQ=255)
                                               tGAGTGCGATAGTGGCATCATTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCGGAGACTTGTGCAGAAATTGGATCTGCATCCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATCCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                                                                                                                        <  1:144830/250‑1 (MQ=255)
                                                gagggCGATAGTGGCATCATTTTTATCCAGTGAAACTAAAATGTCCCCTTGTCAACCGTAGTTCGTTTTTTTATGATTAACGACAGTGACCTTCCCTGAGACTTGTGCAAAAATTGGTTTTGCTTTCCCCGTAACATAGGCGTCATCTGTACTAATCTTTTCTTCTAATTCCTTTGACCAATAGGAAAAGGCACCTAAAACCGAAATAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTa                                                                                                                                                                                                       <  2:711025/246‑1 (MQ=255)
                                                gAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                                                                                                                        >  1:343084/1‑249 (MQ=255)
                                                 ggTGCGATAGTGGCATCAGTTTTACCCAGTGAAACTAAAATGTCCCCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGTTCTGCATTCCCCGTGACATAGGCTTCATCTGTCCTAATCATGTCTTCAAATTCCATTGCCCAATAGGCAGAGCCACCTGAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                                                                                                                                     <  1:335030/250‑1 (MQ=255)
                                                  gTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTCATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                                                                                                                                     >  1:290398/1‑250 (MQ=255)
                                                  gTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATa                                                                                                                                                                                                                              >  1:474697/1‑225 (MQ=255)
                                                                   gTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                                                                                                                                   >  2:537641/1‑251 (MQ=255)
                                                                    ttttATCCAGTGAAACTAAAATGTCACCTTGTTGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                                                                                                                                   >  1:556301/1‑250 (MQ=255)
                                                                     tttATCCAGTGAAACTAAAATGGCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTa                                                                                                                                                                                  >  2:125584/1‑250 (MQ=255)
                                                                          ccAGAGAAAATAAAATGTCACCGTGTCGAACGTATTTCGTATCTTTATGATAAACGGCAGTGACTCTCCCTGGGGTTTGCGCAAAAATTGGAGCTGCTTTCCCCGTGAAATAGGCGGAAGATGAACTAATCATTTCTTATAATTCCATTGACAAATAGGTAAAGGCCCGTGAAAACGTAATAAAAAAAATTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTg                                                                                                                                                                             <  2:65366/250‑1 (MQ=255)
                                                                           cAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGtt                                                                                                                                                                           >  1:688425/1‑251 (MQ=255)
                                                                              tGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTcc                                                                                                                                                                         >  2:685366/1‑250 (MQ=255)
                                                                                 aacaaaaaTGTCCCCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGCCACTCCCTGAGCCTTGTGCAGAAATTGGATCTGCATCCCCCGTGACATAGGCGTCATTTGTACTAATAATGTCTTCAAATTCCATGGACCAATAGGAATAGCCCCCTAAAAACGAAAAAAAAAAAACTCCCGCCAATAAAGaa                                                                                                                                                                                                                         <  2:733361/195‑1 (MQ=255)
                                                                                 aaCTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGaa                                                                                                                                                                                                                         >  1:733361/1‑199 (MQ=255)
                                                                                  attAAAATGTCCCCTTTTCGAACGAATTTGGTATTTTTATGATTAACGAAAGTGACACTCCCTGAGATTTGTGAAAAAATTGGATCTGAATCCCCCGTGACATAGGCGTCATTTGTAATAATCATGTCTTATAATTCCATTGACCAATAGGCATAGGCACCTAAAACCGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtat                                                                                                                                                                    <  1:300689/249‑1 (MQ=255)
                                                                                    aaaaaTGTCCCCTTGTAGAACGTAGTTCGTATTTTTATAATTAACGACAGTGACACTACCTGAGTTTTGTGCAGAAATGGGATCTGCATCCCCCGTGAAATAGGCGTCATCTGTACTATTAATGTCTTCTAATCCCATTACCCAATAGGAATAGGCACCGAAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtatta                                                                                                                                                                  <  1:372039/250‑1 (MQ=255)
                                                                                                 tGTCGAACTTATTTCGTATTTTTATATTTAACGACAGTCACACTCCCTGAATTTTGTAAAAAATTTGGATCTGCTTTCCCCGTGAAATAGGCGTCATCTGTACAAATCATTTCTTCTATTCCCATTGCCCAATAGGCAAAGACACCTAAAAACAAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtcn                                                                                                                                                      <  2:343084/250‑2 (MQ=255)
                                                                                                  ggcGAACGTAGTTCGTATCTTTATGATTACCGACAGTGACCCTACCTGAGACTTGTGCAGAAATTGGATCTGCATCCCCCGTAACATAGGCGTCATCTGTACTAATCATGTCTTCTATTTCCATTGCCCAATAGGAATAGGCACCTAAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCGCATTtctc                                                                                                                                                     <  1:119320/248‑1 (MQ=255)
                                                                                                           aGTTCGTTTTTTTATGATTACCGCCAGTGACCCTCCCTGAGACTTGTGCAGAAATGGGATCTGCTTTCCCCGTGCCATAGGCGTCATTTGTCCTAACCTTGTTTTCTAATTCCATTGACCAATAGGCATAGCCCCCTAAAAACCCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGaaa                                                                                                                                           <  1:554873/251‑1 (MQ=255)
                                                                                                                 tATCTTTATGATTAACGACAGTGACACTTCCGGAGATTTGTGAAGAAATGGAATCTCCTTCCCCCTTGACATGGGCGTCATCTGTACTATTCATGTCTTCTAATTCCATGTACAAATAGGCATAGCCCCCTAAAACCGAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAt                                                                                                                                      <  1:551113/250‑1 (MQ=255)
                                                                                                                  aTCTTTATGATTAACGACAGTGACACTCCCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTCCTAATCATTTCTTCTAATCCCATTGACCAATAGGCATAGGCACCTGAAACCGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtg                                                                                                                                     <  1:537641/250‑1 (MQ=255)
                                                                                                                  aTCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtg                                                                                                                                     >  1:24142/1‑250 (MQ=255)
                                                                                                                          gATTAACGACAGTGACACTCCCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCAATTGACAAATAGGAATAGGCACCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTtcttct                                                                                                                                                                                                           <  1:641238/172‑1 (MQ=255)
                                                                                                                          gATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTtcttct                                                                                                                                                                                                           >  2:641238/1‑172 (MQ=255)
                                                                                                                                              ccTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAttcttc                                                                                                         >  2:37124/1‑250 (MQ=255)
                                                                                                                                                  agaCTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATc                                                                                                    >  1:419796/1‑251 (MQ=255)
                                                                                                                                                   gactggggaaaaaaTTGGATCTCCATCCCCCGTGACATGGGCGTCATCTATAATAATCATGTCTTCAAATTCCATGAACCAAAAGGAATAGGCACCTAAAAACGAAAAAAAAAAACTTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCa                                                                                                   <  2:438358/240‑1 (MQ=255)
                                                                                                                                                       tgtgCAGAAATTGGATCTGCATTCCCCGTGACATAGGCTTCATCTGTACTAATCATGTCTTCTATTTCCATTAACAAAAAGGCATAGCCCCCTAAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAac                                                                                                <  1:498481/250‑1 (MQ=255)
                                                                                                                                                         ggcaaaaaTTGGATCTCCATCCCCCGTGAAATAGGCGTAATCTGTACTAATAATGTCTTCTAATTCCATTGACCAATAGGAATAGGCACCTGAAAACGAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAAcaca                                                                                             <  1:389005/246‑1 (MQ=255)
                                                                                                                                                           aaGAAATTGGATCTGCATCCCCCGTAACAGAGACTTCATCTGTACAAATCATTTCTTCTAATTCCATTGACCAATAGGAAAAGACCCCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGt                                                                                           <  2:24142/250‑1 (MQ=255)
                                                                                                                                                                     aTCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGGATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTg                                                                                 >  1:377400/1‑251 (MQ=255)
                                                                                                                                                                         gCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGcac                                                                              >  1:496863/1‑250 (MQ=255)
                                                                                                                                                                                               cATCTGTACTAATCATCTCTTCTAATCCCATTGACAAATAGGCAAAGGCCCCTAAAAACGCAATAAATAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAgg                                                        <  2:731349/250‑1 (MQ=255)
                                                                                                                                                                                                         aaTCATGTCTTCTAATTCCATTGACCAATAGGTATAGGAACCTGAAAACGGAATAAATAAAACTACCGCTAATAAAGAAAAGTATTATCGTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCAAATTATCACAGGAACGTTGGAGGACTAATTTaa                                              >  2:1893/1‑250 (MQ=255)
                                                                                                                                                                                                         aaTCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGAATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTaa                                              >  1:564369/1‑250 (MQ=255)
                                                                                                                                                                                                                           cATTGCCCAATAGGCATAGGCCCTTAAAAACTAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAaaca                           <  2:81237/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGt          >  1:63761/1‑251 (MQ=255)
                                                                                                                                                                                                                                               caCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTaa                                              >  2:57119/1‑212 (MQ=255)
                                                                                                                                                                                                                                               caCCTAAAAACGCAATAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTaa                                              <  1:57119/212‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCtt     >  1:271540/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    gAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTg  >  1:658307/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           
TTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTG  >  NC_000913/2483010‑2483504

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: