Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,1320GA30.0% 109.2 / 19.2 62A16V (GCT→GTT) rarDputative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (9/33);  new base A (18/0);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 5.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCA  >  NC_000913/4003887‑4004380
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                             
ggTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGGAGTAAAACAACTTTAAGTACGCTGGCGCTATAACCCAAAAAAAAAAAACGGc                                                                                                                                                                                                                                                           >  1:471473/1‑250 (MQ=255)
 gTGATTGTTCACCGACCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                                                                                                                                                                                                                         <  2:388940/251‑1 (MQ=255)
    aTTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGTCCGTAGTAAATCAACTTTAAGTACGCCGGCGCTATACCCCAAAAAAAATAAACGGccaa                                                                                                                                                                                                                                                        >  1:73371/1‑246 (MQ=255)
    aTTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                                                                                                                                                                                                                      <  2:636685/251‑1 (MQ=255)
          cACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                                                                                                                                                                                                                                 <  2:311100/250‑1 (MQ=255)
          cACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                                                                                                                                                                                                                                                <  1:518647/251‑1 (MQ=255)
            ccGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                                                                                                                                                                                                                              <  1:120851/251‑1 (MQ=255)
               cccAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                                                                                                                                                                                                                                            <  1:178219/234‑1 (MQ=255)
                          aGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCGAGTAAATAAAATTTAAATAAGCTGGGGCCATACCCCCAAAAAAATAAACCGCCAGAGGGAGGAATAAGCCCTCCcgcgg                                                                                                                                                                                                                               >  1:830037/1‑250 (MQ=255)
                                     gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                                                                                                                                                                                                                  <  1:179160/251‑1 (MQ=255)
                                     gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                                                                                                                                                                                                                  <  1:208457/251‑1 (MQ=255)
                                     gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAGTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCTCTGGCGGCAAATGCTCATCAGCACCTCCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGTGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATCCCCCAAATAAAATAAGCGTCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                                                                                                                                                                                                                  <  2:407842/251‑1 (MQ=255)
                                           aaTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATAAACTTGAAGTACGCTGGAGCCATAACCCAAATAAAAAAAACGGCAAGAGCGAGTAATACGCCCTGCCGCGGTTGttttttatccctcg                                                                                                                                                                                                             >  1:474907/1‑241 (MQ=255)
                                             tCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAGATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAGCTTGTAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCAGCgg                                                                                                                                                                                                          <  1:727114/250‑1 (MQ=255)
                                               aGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGGAGGAAAACAACTTTAAGTACGCCTGCGCTATACCCCCAATAAAAAAAACAGCAAAAACGGGGAATAAACCCCGCCGCGCTTGTTTTTTATCCATCggggg                                                                                                                                                                                                       >  2:284221/1‑249 (MQ=255)
                                                gCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGGCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGTAGTAAAGCAGCTTGAGGTACGCTGGCGCTATTCCCCAAATAAAATAAGCGGGAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCAGCGGGAAt                                                                                                                                                                                                      <  1:703256/251‑1 (MQ=255)
                                                    ggCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCTCCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTc                                                                                                                                                                                                  <  2:830037/251‑1 (MQ=255)
                                                    ggCAGAGACCGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTCAAAAAACTTGAGATACGCTGGGCCCATACCCCAAAAAAAATAATCGGCCAGCGGCAGGAATAAGCCCTGCCGCGGTTTTTTTTGCTCCATCggggagaacc                                                                                                                                                                                                  >  1:568795/1‑244 (MQ=255)
                                                       agagACTTCCAGCATAAAAATTTTCTGTGGCGGCTGAATCAGCGTTTTTAAATAGGACCACGGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtt                                                                                                                                                                                               <  2:677999/251‑1 (MQ=255)
                                                       agagACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtt                                                                                                                                                                                               <  1:189510/251‑1 (MQ=255)
                                                          gaCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGGTGTAAAACAACTTGAAGTACGCTGGCGCCATACCCCCAAAAAAAAAAAAGGCAAAAGAGGGGAATACGCCCCTCCGCGTTTGTTTTTCCTCCAACGGGGATACACAttttt                                                                                                                                                                                            >  2:41748/1‑251 (MQ=255)
                                                                    aTAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGTAAATCAACCTGGAGTAAGCGGGCGCCATAACCCAAATAAAATAAGCGGCAAAAGCGGGGAAAAAGCCCCGCCGCGGTTGTTTTGCATCCATCGGGAATACTCATTTTTTATTTGTAAc                                                                                                                                                                                  >  1:76792/1‑251 (MQ=255)
                                                                      aaaaaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGAAGTAAAACAACTTGTAGTAAGCTGGGGCTATCACCCCAAAAAAATAAAAGGGCACAGTGGGAAATAACTCCCCCCCCGGTTTTTTTTGATTCATCGGGAATACTCATTTTTTATTGGTAAca                                                                                                                                                                                 >  2:154140/1‑250 (MQ=255)
                                                                                 gtgGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAATCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTTGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCTCGTTTGTTTTGCATCCATCGTGAATACTCATTTTTAATTGGTAACAGCAGTTTAc                                                                                                                                                                       <  1:224153/248‑1 (MQ=255)
                                                                                   ggCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAAcagca                                                                                                                                                                              <  1:389545/240‑1 (MQ=255)
                                                                                   ggCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCGAGGAAATAAACTTGAAGTACGCTTGCGCTCTACCCCCAATAAAAAAAAAGGGAAGAGCGAGGAAAACACCCTGCCGCGGTTGTTTTTTCATCTTCGGGGATACTCATTTTTTAtttgttacaaga                                                                                                                                                                              >  2:389545/1‑230 (MQ=255)
                                                                                          gAATCAGCGTTTTTAAATAGGCCCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCAGCGGGAATACTCATTGTTAATTGGTAACAGCAGTTTACCTGCTTTTATGt                                                                                                                                                           <  1:638188/251‑1 (MQ=255)
                                                                                            aTCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAGTAAATCAACTTGAAGTAAGCTGGCGCTATACCCCCAATAAAATAAGAGGGAAGAGGGAGGAATACGCCCTTCCGCGGTTGTTTTTCATTCATCGGGGATACTCATTTTTAATTGGTAACAGCAGTTTTCCTGCTTTTATGttt                                                                                                                                                         >  1:85577/1‑249 (MQ=255)
                                                                                                    ttttAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCGAGTAAATCAACTTTAAGTAAGCTGGCGCTATACCCCCAATAAAAAAAACGGCCAGAGCGAGTAATAAGCCCCGCCGCGTTTGTCTTGCATCCATCGGGGATACTCATTTTTAATTGGTAACAGCAGGTTACCTGCTTTTATGTCTTCAACCt                                                                                                                                                  >  1:842706/1‑250 (MQ=255)
                                                                                                        aaaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCat                                                                                                                                             <  2:767699/251‑1 (MQ=255)
                                                                                                             ggACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATAAACTTGAAGTAAGCTTGCGCTATACCCCAAATAAAAAAAACGGCCAGAGCGAGGAATACGCCCTGCCGCGTTTGTTTTTCATCCATCGGGGATACTCCTTTTTAATTTGTAACAGCAGTTTACCTTCTTTTTTTTCTTTAAACTACCCTATAgg                                                                                                                                        >  1:709814/1‑251 (MQ=255)
                                                                                                               aCCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCACAGTAAATCAACCTGAAACACGCTTGGGCTCAACCCCCAACAAAAAAAAAGGCCAGAGAGAGGAATAATCCCCGCCGCGGTTGTTTTGTATTCCTCCGGGATACTCAATTTTATTTGGTAAAAACAGATTTCCTTCTTTTTTTTCTTCAACCCACCATATTggtt                                                                                                                                      >  2:454881/1‑250 (MQ=255)
                                                                                                               aCCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAAGTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAggtg                                                                                                                                      <  2:127559/251‑1 (MQ=255)
                                                                                                                           aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                                                                                                                                                                                                                       <  1:336564/146‑1 (MQ=255)
                                                                                                                           aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACAACTTTAAGTACGCCTGCGCCATACCCCCAAAAAAAAAAACGGCAAGAGCGGGTAATACGACCCGCCCCGGTTGTTTTTCATCCATCGGGGATACTCATTTTTAATTTGTAACAGCAGTTTTCCTGCTTTTTTTTCTTCAAACTACCCTATAGGTGGAGGGGGCAccg                                                                                                                          >  1:31890/1‑249 (MQ=255)
                                                                                                                            aaTGCTCATCCGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGGAGTAAATCAACTTGGAAGACGCCGGGGCTTTAACCCAAATAAAAAAAACGGCAAGAGCGAGGAATACGCCCTGCCGCGGTTGTTTTGCATTCATCAGGGAAACTCATTTTTTATTTGTAACAGCCGTTTACATGCTTTTTTTTCTTTAACCTACCCTATAAGTGGCGGTGGCACTgg                                                                                                                         >  1:128873/1‑251 (MQ=255)
                                                                                                                            aaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAAATTTAAGTAAGCTGGCGCTATACCCCCAATAAAATAAAAGGCAAAAGCGGGAAATAAGCCCCACCGCGTTTTTTTTGCATCCATCGGGGATAATCATTTTTTATTGGGAACAACAGTTTTACTGCTTTTTTTTTTTCAAACTTCCCTATACGTGGCGGGGGCACTgg                                                                                                                         >  1:660406/1‑251 (MQ=255)
                                                                                                                                  catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAat                                                                                                                     <  1:284221/249‑1 (MQ=255)
                                                                                                                                   atcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGCAATAAAACAACTTGAAGTAAGCTGGCGCCATACCCCCAAAAAAATAAAAGGCACGAGAGAGTAATACGCCCCGCCGCGGTTGTTTTGCACTCATCGGGAAAACTCATTTTTTATTGGTAACAGCAGTTTACATGCTTTTATGTCTTCTACCTACCATATatgtggggtggggtccgggcaactt                                                                                                                 >  2:542592/1‑237 (MQ=255)
                                                                                                                                        caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAATAAATACACTTGAAGTACGCTGGCGCTATACCCCAAAAAAAATAAACGGCCAGAGCGAGGAATACGCCCCGCCGCGATTTTTTTTCCTCCATCAGGGATACTCCTTTTTTATTTGGAACAACAGTTTTCCTGCTTTTTTTTCTTCCACCTACCATATAAGTGGGGGTGGCACCTGCAATATaaaact                                                                                                             >  2:5641/1‑248 (MQ=255)
                                                                                                                                          ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGt                                                                                                           <  1:346781/251‑1 (MQ=255)
                                                                                                                                                         aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAt                                                                                            <  1:678380/251‑1 (MQ=255)
                                                                                                                                                         aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAt                                                                                            <  1:106442/251‑1 (MQ=255)
                                                                                                                                                          cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTCTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAtt                                                                                           <  2:842706/251‑1 (MQ=255)
                                                                                                                                                          cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGGAAGCTGGCGCTATACCCCAAATAAAAAAAACGGCAAGAGCGAGGAATACGCCCTGCCGCGGTTGTTTTTCATCCCTCGGGAATACTCCTTTTTTATTGGTAACAGCAGTTTACCTGCTTTTATTTCTTCAACCTACCATATAGGTGGCGGGGGGACTGGGAATATAAAGCTGTTCTTCATTGTGTAAtt                                                                                           >  1:844531/1‑251 (MQ=255)
                                                                                                                                                             ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAAAGAATCAACTTAAAGTACGCTGGGGCTATACCCCCAATAAAATAACCAGCAAGAGCGAGTAAAACGCCCTGCCGCGTTTGTTTTGCCTCCATCGGGGATAATCATTTTTAATTTGTAAAAGAAGGTTACCCCCTTTTTTTTCTTTAACCTAACATATAGGTGGGGGAGGCACTGGCTATATAAAGCtattttttattttttaattta                                                                                         >  1:261375/1‑230 (MQ=255)
                                                                                                                                                                             gTCAGGATTTCATCGGCGGGCACGTAGTAAAAAAACTTGAAGTACGCTGGCGCTATACCCCCAATAAAAAAAACGGCAAGAGCGAGGAATACGCCCTGCCGCGTTTGTTTTTCATCCATCGGGAAAAATCATTTTTTATTTGTAACAACAGGTTTACTTCTTTTATTTTTTCAACCTACCCTATAGGTGGCGGGGGCACTGGCAATATAAAGCTGTTCTTCATTTTGTTATTCAAAGCGggcgacggcggc                                                                        >  1:255227/1‑249 (MQ=255)
                                                                                                                                                                                   atttcatcggcgccaaagcagcaaaTCAACTTGAAGTACCATGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGAGAGTACTACGCCCTGCCGCGTTTGTTTTTCATCCAACGGGAAAACTCATTTTTTATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTAACATATAGGCGGCGGTGGCACTGGCAATATAAAGATGTTCTTTCTTGTGTAATTCAACGCGGGCGACGGCGACCTTAtt                                                                  >  2:244322/23‑251 (MQ=255)
                                                                                                                                                                                        aTCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTg                                                             <  2:31890/251‑1 (MQ=255)
                                                                                                                                                                                                        taAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTCCTGCACGCAACAGGctact                                             <  2:844531/250‑1 (MQ=255)
                                                                                                                                                                                                           aTCAACGTGAAGTACGCTTGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTTCCTCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACtagt                                          <  1:448920/251‑1 (MQ=255)
                                                                                                                                                                                                              aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGc                                       <  1:471184/251‑1 (MQ=255)
                                                                                                                                                                                                              aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGc                                       <  1:542592/251‑1 (MQ=255)
                                                                                                                                                                                                                  tGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACTCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAg                                     <  2:474014/249‑1 (MQ=255)
                                                                                                                                                                                                                       gTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAAc                               >  1:134124/1‑251 (MQ=255)
                                                                                                                                                                                                                         aCGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACg                              <  2:85577/250‑1 (MQ=255)
                                                                                                                                                                                                                          cGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCt                            >  2:29349/1‑251 (MQ=255)
                                                                                                                                                                                                                            cTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCg                          >  1:503128/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaataaaataAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAggt                                                                                                                                       <  2:480564/127‑1 (MQ=255)
                                                                                                                                                                                                                                            aaataaaataAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAggt                                                                                                                                       >  1:480564/1‑127 (MQ=255)
                                                                                                                                                                                                                                                    taaGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATc   >  1:42869/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      aGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCa  >  1:412097/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                             
GGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCA  >  NC_000913/4003887‑4004380

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