Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,2680TA31.2% 83.4 / 18.6 48L24F (TTA→TTTemrKtripartite efflux pump membrane fusion protein EmrK
Reads supporting (aligned to +/- strand):  ref base T (26/7);  new base A (0/15);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.56e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTG  >  NC_000913/2483021‑2483504
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
gCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATATGCGTCATCTGTACTAATCATGTCTGCTAATTTCATTGACCAATAGGCATAGGCACCTGAAAACGCAataaataa                                                                                                                                                                                                                                             <  1:394241/250‑1 (MQ=255)
   gAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACt                                                                                                                                                                                                                                         >  1:793341/1‑251 (MQ=255)
    aaCAATATTTGCCAGATTATTTTTAGCTTTAATGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTGCTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACCCAATAAATAAAACTa                                                                                                                                                                                                                                       <  2:600024/251‑1 (MQ=255)
     aCAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCATTTTTATCCAGTGAAACTAAAATGTCGCCTTGTCGAACGTAGTTCGTATCTGTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGGGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATGTACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTAc                                                                                                                                                                                                                                      <  1:913255/251‑1 (MQ=255)
     aCAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTCCAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGCCTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTAc                                                                                                                                                                                                                                      >  1:876269/1‑251 (MQ=255)
                      attTTTAGCTTTATTGGGTGCGATAGTGGCATCAGTTTTATCCAGTGAACCTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATATGCATTCCCCGTGACATCGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGt                                                                                                                                                                                                                     >  2:576068/1‑251 (MQ=255)
                          ttAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTtttt                                                                                                                                                                                                                 >  2:635522/1‑251 (MQ=255)
                              cTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTttctt                                                                                                                                                                                                             >  2:249129/1‑251 (MQ=255)
                                 tattgattcgaatattggaaTCAGTTTTATCCATTAAAACTAAAATGTCCCTTTGTCGACCGTAGTTCTTATTTTTATGATTAACGACAGGGACACTCCTGGAGCCTTGTGAAGAAATTGTATCTGCATCCCCCGTGCCAGAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCTTTGACCAATGGGCATAGCCCCCTAAAAACGAAATAAATAAACCTACCGCCAATAAAGAAAAGTATTTtcttct                                                                                                                                                                                                           <  1:668819/231‑1 (MQ=255)
                                  aTTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGACACTAAACTGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGt                                                                                                                                                                                                         >  2:470719/1‑251 (MQ=255)
                                     tccagccgtgtgctctccctatttACCCAGTGAAACAAAAATGTCCCTTTGCCGAACTTAGTTCGTATCTTTATGATTAACGACAGTGACACTCCCGGAGACTTGGGAAAAAATGGGATCTGCATCCCCCGTGACAAGGGTGTCATCTGTACTAATATTGTCTTCAAATTCCATGCACCAATAGGCAAAGGCACCTAAAAACTCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTa                                                                                                                                                                                                       <  1:948118/229‑1 (MQ=255)
                                       gTGCGATAGGGGCATCAGTTTAACCCAGTGAAACAAAAATTCCCCCTGGTCGAACGTGGTTCGTTTCTTTAGAATAAACGACAGTGACCCTCCCGGAGATTTGTGCAAAAATTGGATCGGCATCCCCCGGGACAAAGGCTTGATCTTTACAAATGATGTCTTCAATGTCCATTGACCAAGGGGAAAAGGCCCCTAAAACCGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGaa                                                                                                                                                                                                    <  1:470719/251‑1 (MQ=255)
                                              aGTGGCATCAGTTTTACCCAGTGAAACAAAAATGTCACTTTGTCGACCTTAGTTCGTATCTTTATGATTAACGACAGTGACTCTCCCGGAGACTTGTGCAAAAATTGGATCTGCATTCCCCGTAACAGAGGCGTCATCTGTACTAATCTTGTTTTCTAATCCCATTGCCCAATAGGCATAGGCACCTGAAAACGAAAAAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttt                                                                                                                                                                                              <  1:442041/250‑1 (MQ=255)
                                                 ggCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTa                                                                                                                                                                                                       >  1:754499/1‑238 (MQ=255)
                                                  gCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAtt                                                                                                                                                                                         >  2:656569/1‑251 (MQ=255)
                                                    aTCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTg                                                                                                                                                                                       >  2:434074/1‑251 (MQ=255)
                                                         ttttACCCAGTGAAAATAAAAGTGCCCCTTGTCGACATTAGTTCTATTTTTAATGATTACCGACAGTGACACAACCGGAGACTTGTGAAAAATTTTGATCTGTATCCCCCTTGACAGAGGCTTTATGTGTACTAATCTTGTCTCCTAATCAAATTGCCCAATAGGCAAAGGCACTAAAAACCGCAAAAAAAAAACTTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                                                                                                                                   <  2:801160/250‑1 (MQ=255)
                                                         ttttACCCAGTAAAACAAAAATGTCCCCTTGTAGAACGTAGTTCGTTTTTTTATGATAAACGACATTGCCACTACCTGAGACTTGGGAAGAAATTGGATCGGCATCCCCGGTGACGTGGGCGTCATCTGAATTAATCATGTTTTTTAATCCAATTCACCAAAGGGAAGAGGCACCTGAAAACCAAATAAATAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                                                                                                                                   <  1:434074/250‑1 (MQ=255)
                                                           ttATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAAt                                                                                                                                                                                >  2:30794/1‑251 (MQ=255)
                                                             aTCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTa                                                                                                                                                                                                       >  2:948118/1‑226 (MQ=255)
                                                                    tgaaataaaaagctccccTTGTCGAACGTAGTTCGATTCTTTATGATAAACGACAGTGCCACTCCGTGAGCTTTGGACAGAAATGGGATCGGAATCCCCCTTAAAATAGGCTCCATCTATACTAATCATGTCTTCTAATTCCATTACCAAATAGGCATAGACACCTAAAAACTCAAAAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTcc                                                                                                                                                                         <  2:159909/234‑1 (MQ=255)
                                                                    gAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCAc                                                                                                                                                                       >  1:541365/1‑251 (MQ=255)
                                                                    aaaaaaaaaaatgtccctttgcagAACGTAGTTCGTATCTTTGTGATTAACGACAGTGCCTCTACCTGAGACTTGTGCAGAAATTGGATCTCCATTCCCCGTGACACAGTCGTCATCTGTAATAATAATTTCTTCTAATCCCATTCACCAAGAGACAAAGGCACTTAAAAACTCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCa                                                                                                                                                                        <  2:155289/228‑1 (MQ=255)
                                                                      aaacaaaatggcctctttGTCAAACGTAGTTCGTTTCTTTATGATTAACGACAGTGCCTCTACCTGAGGTTTGTGCAAAAATGGGATCTGCATTCCCCGTGACATAGGCGTCTTCTTTTCAAATCATGTCTTCTATTCCCTTTGCCAAAGAGGCATAGGCACCTAAAAACGAAATAAATAAACCTGCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACt                                                                                                                                                                      <  1:635522/235‑1 (MQ=255)
                                                                       aCTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGAATTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGAACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtat                                                                                                                                                                    >  1:499339/1‑251 (MQ=255)
                                                                                                             aaTGATTAACGACATTCACACTACCTTAGTCTTGAGAAGAAATTGTACCTCCATCCCCCGTGAAATAGGTGTCATCTGTTCAAATCATGTCTTCTAATTCCATTCACCAATAGTCAAAGCCACCTAAAACCGCAATAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTACTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg                                                                                                                                <  2:279268/249‑1 (MQ=255)
                                                                                                                 gATTAACGAAAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGACTGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATtgaatg                                                                                                                            >  2:585423/1‑251 (MQ=255)
                                                                                                                     accgaaagggaccctcccTGAGCCTTGTGCATAAATTGGATCTCAATCCCCCGTGACATGGCCGCCTTTTGTTTAAATCATCTTTACTAATTCCATACACCAATAGGAATAGCCACCTAAAATCGTAATAAAAAAATCAACCGCCAATAAATAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAGTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAAGTTATGAATGGAATGaaaa                                                                                                                        <  1:656569/235‑1 (MQ=255)
                                                                                                                     aaCGACAGTGACACTACCTGCGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGGACTACTCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAGATAAAACTACAGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaa                                                                                                                         >  2:422010/1‑250 (MQ=255)
                                                                                                                            gTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcat                                                                                                                  >  1:926310/1‑250 (MQ=255)
                                                                                                                                      ccggagctttGTGCAAAATTGGGATCTCCTTTCCCCGTGACATAGGCGTCATCTGTACTAATCATTTTTTCAAATCCCATTAACAAATAGGAATAGGCACCTAAAAACAAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCa                                                                                                        <  1:349930/243‑1 (MQ=255)
                                                                                                                                                      aaaTTGGATCTGCATCCCCCGGGACATAGGCGTCATCTGTACAAATCATTTCTTCAAATCCAATTGCCCAATAGGCATAGGCACCTAAAAACGAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGGAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCa                                                                                        <  1:30794/251‑1 (MQ=255)
                                                                                                                                                            gATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGtt                                                                                  >  1:115860/1‑251 (MQ=255)
                                                                                                                                                                      cccccGTGACATAGGCGTCTTTTTTACTAATCATTTCTTCAAATTCCATTAACAAATAGCAATAGCCACCTAAAACCGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTa                                                                          <  1:295768/249‑1 (MQ=255)
                                                                                                                                                                                    gCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTgg                                                           >  1:168926/1‑251 (MQ=255)
                                                                                                                                                                                      gTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATCAAACTACCGCCAATCAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAgg                                                                   >  1:399863/1‑241 (MQ=255)
                                                                                                                                                                                         aTCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGcc                                                      >  2:837505/1‑251 (MQ=255)
                                                                                                                                                                                         aTCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGCAGGcc                                                      >  2:835181/1‑251 (MQ=255)
                                                                                                                                                                                              aacaaaacATGCCTTCGAATCCCAATGACCAAGAGGAATAGCCACGAAAAAAAGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCGCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAtt                                                 <  1:835181/244‑1 (MQ=255)
                                                                                                                                                                                                 tAATCATGTCTTCTAATTCAATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATACAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGCGGCCTACTTTa                                               >  2:583890/1‑250 (MQ=255)
                                                                                                                                                                                                         tcttctAATTCAATTGACCAATAGGCAAAGGCACCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGt                                       <  1:249129/250‑1 (MQ=255)
                                                                                                                                                                                                          cttctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGGAATAAATAACCCTACCGCAAATCAAAACACGTATTTTCATCTGTTAGAATGTTTATTCTTTGCATTACTCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAACGAGCATTCATTCTTCATATCCAACACAGTGCATGTGTTGCACATTATCCACGGAACGTTGGAGGCCTAATTTAAATGAGGTGc                                     >  1:887570/1‑251 (MQ=255)
                                                                                                                                                                                                              aaaatcccATTAACCAATAGGCATAGGCACCGAAAAACCAAAAAAAAAAAACTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaaa                                  <  2:168926/245‑1 (MQ=255)
                                                                                                                                                                                                                                    cataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAAc             >  1:610196/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gcaCCTAAAAACGAAAAAAAAAAATCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTgg        <  1:253413/251‑1 (MQ=255)
                                                                                                                                                                                                                                              caaaaaCGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCttt    <  1:287549/248‑1 (MQ=255)
                                                                                                                                                                                                                                               agaaaccgaaataaaaaAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACTCAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCtt     <  1:585423/240‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaaaCGCAAAAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATTAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTAGCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTg  <  1:583890/250‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaaCGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCttt    >  2:287549/1‑248 (MQ=255)
                                                                                                                                                                                                                                                       cAATAAATAAAACTACCGCCAATAAAGCAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCt                                                     >  1:815572/1‑193 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
GCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTG  >  NC_000913/2483021‑2483504

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