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.

GGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCA  >  NC_000913/4003887‑4004380
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                             
ggTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGGAGTAAAACAACTTT‑AAGTACGCTGGCGCTAT‑AACCCAAAAAAAAAAAACGGc                                                                                                                                                                                                                                                           >  1:471473/1‑250 (MQ=255)
 gTGATTGTTCACCGACCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCaa                                                                                                                                                                                                                                                         <  2:388940/251‑1 (MQ=255)
    aTTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGTCCGTAGTAAATCAACTTT‑AAGTACGCCGGCGCTAT‑ACCCCAAAAAAAATAAACGGccaa                                                                                                                                                                                                                                                        >  1:73371/1‑246 (MQ=255)
    aTTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAagag                                                                                                                                                                                                                                                      <  2:636685/251‑1 (MQ=255)
          cACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                                                                                                                                                                                                                                 <  2:311100/250‑1 (MQ=255)
          cACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                                                                                                                                                                                                                                                <  1:518647/251‑1 (MQ=255)
            ccGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                                                                                                                                                                                                                              <  1:120851/251‑1 (MQ=255)
               cccAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCgg                                                                                                                                                                                                                                                            <  1:178219/234‑1 (MQ=255)
                          aGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCGAGTAAATAAAATTT‑AAATAAGCTGGGGCCAT‑ACCCCCAAAAAAATAAACCGCCAGAGGGAGGAATA‑AGCCCTCCcgcgg                                                                                                                                                                                                                               >  1:830037/1‑250 (MQ=255)
                                     gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑cat                                                                                                                                                                                                                  <  1:179160/251‑1 (MQ=255)
                                     gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑cat                                                                                                                                                                                                                  <  1:208457/251‑1 (MQ=255)
                                     gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAGTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCTCTGGCGGCAAATGCTCATCAGCACCTCCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGTGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑CCCCCAAATAAAATAAGCGTCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑cat                                                                                                                                                                                                                  <  2:407842/251‑1 (MQ=255)
                                           aaTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATAAACTTG‑AAGTACGCTGGAGCCAT‑AACCCAAATAAAAAAAACGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑GTTGttttttatccctcg                                                                                                                                                                                                             >  1:474907/1‑241 (MQ=255)
                                             tCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAGATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAGCTTG‑TAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCAGC‑gg                                                                                                                                                                                                          <  1:727114/250‑1 (MQ=255)
                                               aGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGGAGGAAAACAACTTT‑AAGTACGCCTGCGCTAT‑ACCCCCAATAAAAAAAACAGCAAAAACGGGGAATA‑AACCCCGCCGCG‑CTTGTTTTT‑TATCCATC‑ggggg                                                                                                                                                                                                       >  2:284221/1‑249 (MQ=255)
                                                gCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGGCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGTAGTAAAGCAGCTTG‑AGGTACGCTGGCGCTAT‑TCCCCAAATAAAATAAGCGGGAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCAGC‑GGGAAt                                                                                                                                                                                                      <  1:703256/251‑1 (MQ=255)
                                                    ggCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCTCCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTc                                                                                                                                                                                                  <  2:830037/251‑1 (MQ=255)
                                                    ggCAGAGACCGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTCAAAAAACTTG‑AGATACGCTGGGCCCAT‑ACCCCAAAAAAAATAATCGGCCAGCGGCAGGAATA‑AGCCCTGCCGCGGTTTTTTTTG‑C‑TCCATC‑ggggagaacc                                                                                                                                                                                                  >  1:568795/1‑244 (MQ=255)
                                                       agagACTTCCAGCATAAAAATTTTCTGTGGCGGCTGAATCAGCGTTTTTAAATAGGACCACGGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCAtt                                                                                                                                                                                               <  2:677999/251‑1 (MQ=255)
                                                       agagACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCAtt                                                                                                                                                                                               <  1:189510/251‑1 (MQ=255)
                                                          gaCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGGTGTAAAACAACTTG‑AAGTACGCTGGCGCCAT‑ACCCCCAAAAAAAAAAAAGGCAAAAGAGGGGAATA‑CGCCCCTCCGCG‑TTTGTTTTT‑CCTCCAAC‑GGGGATACACAttttt                                                                                                                                                                                            >  2:41748/1‑251 (MQ=255)
                                                                    aTAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGTAAATCAACCTG‑GAGTAAGCGGGCGCCAT‑AACCCAAATAAAATAAGCGGCAAAAGCGGGGAAAA‑AGCCCCGCCGCG‑GTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTTATTTGTAAc                                                                                                                                                                                  >  1:76792/1‑251 (MQ=255)
                                                                      aaaaaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGAAGTAAAACAACTTG‑TAGTAAGCTGGGGCTATCACCCCAAA‑AAAATAAAAGGGCACAGTGGGAAATAACTCCCCCCCCGG‑TTTTTTTTG‑‑ATTCATC‑GGGAATACTCATTTTTTATTGGTAAca                                                                                                                                                                                 >  2:154140/1‑250 (MQ=255)
                                                                                 gtgGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAATCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTTGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCTCG‑TTTGTTTTG‑CATCCATC‑GTGAATACTCATTTTTAATTGGTAACAGCAG‑TTTAc                                                                                                                                                                       <  1:224153/248‑1 (MQ=255)
                                                                                   ggCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAAcagca                                                                                                                                                                              <  1:389545/240‑1 (MQ=255)
                                                                                   ggCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCGAGGAAATAAACTTG‑AAGTACGCTTGCGCTCT‑ACCCCCAATAAAAAAAAAGGGAAGAGCGAGGAAAA‑CACCCTGCCGCG‑GTTGTTTTTTCATC‑TTC‑GGGGATACTCATTTTTTAtttgttacaaga                                                                                                                                                                              >  2:389545/1‑230 (MQ=255)
                                                                                          gAATCAGCGTTTTTAAATAGGCCCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCAGC‑GGGAATACTCATTGTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGt                                                                                                                                                           <  1:638188/251‑1 (MQ=255)
                                                                                            aTCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAGTAAATCAACTTG‑AAGTAAGCTGGCGCTAT‑ACCCCCAATAAAATAAGAGGGAAGAGGGAGGAATA‑CGCCCTTCCGCG‑GTTGTTTTT‑CATTCATC‑GGGGATACTCATTTTTAATTGGTAACAGCAG‑TTTTCCTGCTTTTATGttt                                                                                                                                                         >  1:85577/1‑249 (MQ=255)
                                                                                                    ttttAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCGAGTAAATCAACTTT‑AAGTAAGCTGGCGCTAT‑ACCCCCAATAAAAAAAACGGCCAGAGCGAGTAATA‑AGCCCCGCCGCG‑TTTGTCTTG‑CATCCATC‑GGGGATACTCATTTTTAATTGGTAACAGCAG‑GTTACCTGCTTTTATGTCTTCAACCt                                                                                                                                                  >  1:842706/1‑250 (MQ=255)
                                                                                                        aaaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCat                                                                                                                                             <  2:767699/251‑1 (MQ=255)
                                                                                                             ggACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATAAACTTG‑AAGTAAGCTTGCGCTAT‑ACCCCAAATAAAAAAAACGGCCAGAGCGAGGAATA‑CGCCCTGCCGCG‑TTTGTTTTT‑CATCCATC‑GGGGATACTCCTTTTTAATTTGTAACAGCAG‑TTTACCTTCTTTTTTTTCTTTAAACTACCCTATAgg                                                                                                                                        >  1:709814/1‑251 (MQ=255)
                                                                                                               aCCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCACAGTAAATCAACCTG‑AAACACGCTTGGGCTCA‑ACCCCCAACAAAAAAAAAGGCCAGAGAGAGGAATA‑ATCCCCGCCGCG‑GTTGTTTTG‑TATTCCTC‑CGGGATACTCAATTTTATTTGGTAAAAACAG‑ATTTCCTTCTTTTTTTTCTTCAACCCACCATATTggtt                                                                                                                                      >  2:454881/1‑250 (MQ=255)
                                                                                                               aCCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAAGTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAggtg                                                                                                                                      <  2:127559/251‑1 (MQ=255)
                                                                                                                           aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGccc                                                                                                                                                                                                                                       <  1:336564/146‑1 (MQ=255)
                                                                                                                           aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACAACTTT‑AAGTACGCCTGCGCCAT‑ACCCCCAAAAAAAAAAACGGCAAGAGCGGGTAATA‑CGACCCGCCCCG‑GTTGTTTTT‑CATCCATC‑GGGGATACTCATTTTTAATTTGTAACAGCAG‑TTTTCCTGCTTTTTTTTCTTCAAACTACCCTATAGGTGGAGGGGGCAccg                                                                                                                          >  1:31890/1‑249 (MQ=255)
                                                                                                                            aaTGCTCATCCGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGGAGTAAATCAACTTGGAAG‑ACGCCGGGGCTTT‑AACCCAAATAAAAAAAACGGCAAGAGCGAGGAATA‑CGCCCTGCCGCG‑GTTGTTTTG‑CATTCATCAGGGAA‑ACTCATTTTTTATTTGTAACAGCCG‑TTTACATGCTTTTTTTTCTTTAACCTACCCTATAAGTGGCGGTGGCACTgg                                                                                                                         >  1:128873/1‑251 (MQ=255)
                                                                                                                            aaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAAATTT‑AAGTAAGCTGGCGCTAT‑ACCCCCAATAAAATAAAAGGCAAAAGCGGGAAATA‑AGCCCCACCGCG‑TTTTTTTTG‑CATCCATC‑GGGGATAATCATTTTTTATTGGGAACAACAG‑TTTTACTGCTTTTTTTTTTTCAAACTTCCCTATACGTGGCGGGGGCACTgg                                                                                                                         >  1:660406/1‑251 (MQ=255)
                                                                                                                                  catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAat                                                                                                                     <  1:284221/249‑1 (MQ=255)
                                                                                                                                   atcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGCAATAAAACAACTTG‑AAGTAAGCTGGCGCCAT‑ACCCCCAAAAAAATAAAAGGCACGAGAGAGTAATA‑CGCCCCGCCGCG‑GTTGTTTTG‑CACTCATC‑GGGAAAACTCATTTTTTATTGGTAACAGCAG‑TTTACATGCTTTTATGTCTTCTACCTACCATATatgtgg‑ggtggggtccgggcaactt                                                                                                                 >  2:542592/1‑237 (MQ=255)
                                                                                                                                        caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAATAAATACACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAAAAAAATAAACGGCCAGAGCGAGGAATA‑CGCCCCGCCGCG‑ATTTTTTTT‑CCTCCATC‑AGGGATACTCCTTTTTTATTTGGAACAACAG‑TTTTCCTGCTTTTTTTTCTTCCACCTACCATATAAGTGGGGGTGGCACCTGCAATATaaaact                                                                                                             >  2:5641/1‑248 (MQ=255)
                                                                                                                                          ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGt                                                                                                           <  1:346781/251‑1 (MQ=255)
                                                                                                                                                         aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAt                                                                                            <  1:678380/251‑1 (MQ=255)
                                                                                                                                                         aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAt                                                                                            <  1:106442/251‑1 (MQ=255)
                                                                                                                                                          cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTCTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAtt                                                                                           <  2:842706/251‑1 (MQ=255)
                                                                                                                                                          cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTG‑AAGGAAGCTGGCGCTAT‑ACCCCAAATAAAAAAAACGGCAAGAGCGAGGAATA‑CGCCCTGCCGCG‑GTTGTTTTT‑CATCCCTC‑GGGAATACTCCTTTTTTATTGGTAACAGCAG‑TTTACCTGCTTTTATTTCTTCAACCTACCATATAGGTGGCGGGGGGACTGGGAATATAAAGCTGTTCTTCATTGTGTAAtt                                                                                           >  1:844531/1‑251 (MQ=255)
                                                                                                                                                             ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAAAGAATCAACTTA‑AAGTACGCTGGGGCTAT‑ACCCCCAATAAAATAACCAGCAAGAGCGAGTAAAA‑CGCCCTGCCGCG‑TTTGTTTTG‑CCTCCATC‑GGGGATAATCATTTTTAATTTGTAAAAGAAG‑GTTACCCCCTTTTTTTTCTTTAACCTAACATATAGGTGGGGGAGGCACTGGCTATATAAAGCtattttttattttttaattta                                                                                         >  1:261375/1‑230 (MQ=255)
                                                                                                                                                                             gTCAGGATTTCATCGGCGGGCACGTAGTAAAAAAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCCAATAAAAAAAACGGCAAGAGCGAGGAATA‑CGCCCTGCCGCG‑TTTGTTTTT‑CATCCATC‑GGGAAAAATCATTTTTTATTTGTAACAACAGGTTTAC‑TTCTTTTATTTTTTCAACCTACCCTATAGGTGGCGGGGGCACTGGCAATATAAAGCTGTTCTTCATTTTGTTATTCAAAGCGggcgacggcggc                                                                        >  1:255227/1‑249 (MQ=255)
                                                                                                                                                                                   atttcatcggcgccaaagcagcaaaTCAACTTG‑AAGTACCATGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGAGAGTACTA‑CGCCCTGCCGCG‑TTTGTTTTT‑CATCCAAC‑GGGAAAACTCATTTTTTATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTAACATATAGGCGGCGGTGGCACTGGCAATATAAAGATGTTCTTTCTTGTGTAATTCAACGCGGGCGACGGCGACCTTAtt                                                                  >  2:244322/23‑251 (MQ=255)
                                                                                                                                                                                        aTCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTg                                                             <  2:31890/251‑1 (MQ=255)
                                                                                                                                                                                                        taAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATT‑CCTGCACGCAACAGGctact                                             <  2:844531/250‑1 (MQ=255)
                                                                                                                                                                                                           aTCAACGTG‑AAGTACGCTTGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTTCCTCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACtagt                                          <  1:448920/251‑1 (MQ=255)
                                                                                                                                                                                                              aaCTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGc                                       <  1:471184/251‑1 (MQ=255)
                                                                                                                                                                                                              aaCTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGc                                       <  1:542592/251‑1 (MQ=255)
                                                                                                                                                                                                                  tG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CTCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAg                                     <  2:474014/249‑1 (MQ=255)
                                                                                                                                                                                                                       gTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAAc                               >  1:134124/1‑251 (MQ=255)
                                                                                                                                                                                                                         aCGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACg                              <  2:85577/250‑1 (MQ=255)
                                                                                                                                                                                                                          cGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCt                            >  2:29349/1‑251 (MQ=255)
                                                                                                                                                                                                                            cTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCg                          >  1:503128/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaataaaataAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAggt                                                                                                                                       <  2:480564/127‑1 (MQ=255)
                                                                                                                                                                                                                                            aaataaaataAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAggt                                                                                                                                       >  1:480564/1‑127 (MQ=255)
                                                                                                                                                                                                                                                    taaGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATc   >  1:42869/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      aGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCa  >  1:412097/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                             
GGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATA‑CGCCCTGCCGCG‑TTTGTTTTG‑CATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCA  >  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: ‑