Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,1320GA26.0% 101.4 / 10.0 51A16V (GCT→GTT) rarDputative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (10/27);  new base A (13/0);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 3.22e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATC  >  NC_000913/4003896‑4004379
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         
tCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAg                                                                                                                                                                                                                                              <  1:184951/250‑1 (MQ=255)
                       caggtGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt                                                                                                                                                                                                                      <  1:463833/247‑1 (MQ=255)
                            gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                                                                                                                                                                                                                 <  1:267747/251‑1 (MQ=255)
                                ccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                                                                                                                                                                                                               <  1:479565/249‑1 (MQ=255)
                                 caATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc                                                                                                                                                                                                            <  1:436885/251‑1 (MQ=255)
                                      aGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAAACAACTT‑TAAGTA‑CGCTGGCGCTATACCCCCAATAAAATAAAAGGCAAGAGCGGGTAATACGCCCCCCCGCGTTTGTTTTGCCTCCATCgggg                                                                                                                                                                                                        >  1:568096/1‑249 (MQ=255)
                                        cACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa                                                                                                                                                                                                     <  2:581723/251‑1 (MQ=255)
                                           ggCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGACCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGTGAATACTc                                                                                                                                                                                                  <  1:447115/251‑1 (MQ=255)
                                                  aCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtttt                                                                                                                                                                                             <  2:122496/249‑1 (MQ=255)
                                                  aCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAAAGCGAGTAATACGCCCCGCCGCGTTTGTTTTGCATCCATCGGGAATAATCCtttttt                                                                                                                                                                                           >  1:122496/1‑250 (MQ=255)
                                                        aGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTGTTAATTGGt                                                                                                                                                                                     <  2:21618/251‑1 (MQ=255)
                                                             aaaaaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACAACTT‑GAAGTA‑AGCTGGCGCTATACCCCAAATAAAATAAACGGCAAAAGCGGGGAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGGATACTCCTTTTTTATTGGTAAc                                                                                                                                                                                  >  1:610542/1‑249 (MQ=255)
                                                               aaaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAAca‑gca                                                                                                                                                                             <  2:568096/251‑1 (MQ=255)
                                                                     tCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTa                                                                                                                                                                       <  1:226978/250‑1 (MQ=255)
                                                                        gtgGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTg                                                                                                                                                                   <  1:14078/251‑1 (MQ=255)
                                                                           gCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCtt                                                                                                                                                                <  1:252448/251‑1 (MQ=255)
                                                                                           ttttAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCt                                                                                                                                                 <  1:285196/250‑1 (MQ=255)
                                                                                           ttttAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGGAAAACAAATT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAAAAAACGGCAAGAGCGAGGAATAAGCCCTGCCGCGTTTGTTTTTCATCCATCGGGAATAATCATTTTTTATTTGTAACA‑GCAG‑TTTAACCGCTTTTATTTCTTTAACCt                                                                                                                                                 >  1:209789/1‑250 (MQ=255)
                                                                                              tAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTAcc                                                                                                                                              <  1:294570/250‑1 (MQ=255)
                                                                                               aaaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATAAACAT‑GAAGTA‑AGCGGGGGGCATACCCCCAAAAAAAAAAACGGCAAGAACAGGGAATAAGCCCTGCCCCGTTTGTTTTGCATCCATCAGGGATAATCATTTTTTATTGGTAAGA‑GCAGGTTTCCCT‑CTTTTATGTCTTTAACCTACCac                                                                                                                                            >  2:33705/1‑250 (MQ=255)
                                                                                                aaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGGGCTATACCCCCAATAAAATAAACGGCAAGAGCGGGGAATACGCCCCGCCGCGTTTGTTTTGCATCCATCGGGGATACTCATTTTTTATTGGTAAAA‑GCAG‑TTTTCCTGCTTTTTTGTTTTCAACCTACCat                                                                                                                                            >  1:239693/1‑250 (MQ=255)
                                                                                                               ggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCACATT‑GAAGAAACGC‑GGCGCTATACCCCCCATAAAAAAACAAGCAAGAGCGAGGAAAACGCCCTCCCCCGTTTGTTTTGTCTCCCTCGGGGATAATCCTTTTTTATTTGTAACACGC‑G‑TTTAACTGCTTTTATGTCTTCCACCtccccattaattggcgggggg                                                                                                                             >  2:163497/1‑230 (MQ=255)
                                                                                                                  aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                                                                                                                                                                                                             <  2:71442/169‑1 (MQ=255)
                                                                                                                  aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGGAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAACGGCAAGAGCGAGGAATACGCCCCGCCGCGGTTGTTTTGcatccat                                                                                                                                                                                                             >  1:71442/1‑169 (MQ=255)
                                                                                                                      gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCa                                                                                                                      <  1:219896/250‑1 (MQ=255)
                                                                                                                      gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAAACAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAACGGCAAGAGCGGGTAATAAGCCCTGCCGCGTTTTTTTTTCATCCCTCGGGGATACTCATTTTTTATTGGTAACA‑ACAG‑TTTACCTGCTTTTTTGTCTTCAACCTACCATATAGGTGGCGGGGGGACTGGCa                                                                                                                      >  1:611349/1‑250 (MQ=255)
                                                                                                                         catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCGAGGAAATCAACTTTGAA‑CA‑CGCTGGCGCTTTTCCCCAAATAAAATAAACGGACACAGCGAGTAATACCCCCCGCCGCGCTTTATTTGCCTCCCCCGGGCATAATCACTTTTTATTGGTTACA‑AGAG‑TTTACCTGCTTTTTTTTCTTCCACCTACCCCATACGCGGCGGTGGCACTTGCAatat                                                                                                                  >  1:621319/1‑251 (MQ=255)
                                                                                                                          atcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGAAGGAAATCAACTT‑TAAGTA‑CGCTTGCGCTATAACCCAAAAAAAAAAAACGGCAAGAGCGAGGAATAAGCCCCTCCGCGTTTGTTTTGCATCCATCGGGAAAACTCATTTTTTATTTGTAAAA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTCCCCTATAGGTGGCGGGGGCACTggtgatat                                                                                                                  >  1:98407/1‑244 (MQ=255)
                                                                                                                               caccacCATAAAGAAAAACGACCAGATCCCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTTGCGCTATACCCCAAATAAAATACGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTTTTTTGCATCCATCTGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCt                                                                                                            <  2:485901/251‑1 (MQ=255)
                                                                                                                                accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTg                                                                                                           <  1:33705/251‑1 (MQ=255)
                                                                                                                                  caccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGGA‑CGCTGGCGCTATACCCCAAATAAAATAAAAGGCAAAAGCGGGTAATACGCCCCGCCCCGTTTGTTTTTCATCCATCGGGAATACTCATTTTTTATTTGTAACA‑ACAG‑TTTACCTGCTTTTATTTTTTTCACCTACCATATAGGTGGCGGGGGGCCCGGCCATATAAAACTGtt                                                                                                         >  1:338154/1‑251 (MQ=255)
                                                                                                                                      aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGttcttc                                                                                                     <  2:239693/251‑1 (MQ=255)
                                                                                                                                                        aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAATAAATCAACTT‑GGAGAA‑CGCTGGCGCTATACCCCCAAAAAAAAAAACGGCCAGAGCGGGTAATACGCCCTGCCGCGTTTTTTTTTTATCCCACGGGGATATTCATTTTTAATTGGTAAAA‑ACAG‑TTTACCTGCTTTTTTTTTTTTAACCTTCCATATAGGTGGCGGCGGCCCTTGCAATATAAAGATGTTTTTTCTTGTGTAAtttaaagcg                                                                                   >  1:397503/1‑244 (MQ=255)
                                                                                                                                                           aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGggc                                                                                <  1:630217/251‑1 (MQ=255)
                                                                                                                                                            cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGggcg                                                                               <  2:611349/251‑1 (MQ=255)
                                                                                                                                                             gcgATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTT‑GAAGGA‑CGCCGGCGCTATACCCCAAATAAAATAAACGGCAAGAGCGAGTAATACGCCCCGCCGCGGTTGTTTTGCATCCATCGGGAATACTCCTTTTTTATTTGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGGACTGGCAATATAAAGCTGTTCTTTATTGTGGAATTTAAAGCGggcgg                                                                              >  1:269307/1‑250 (MQ=255)
                                                                                                                                                                         gATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTAt                                                                  <  1:373424/251‑1 (MQ=255)
                                                                                                                                                                                  ggcggGCCCGAAGTAAAACAACTT‑GAAGTA‑CGCTGGCGCTGTACCCCAAATAAAATAAGCGGTAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCAGCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTCTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCAc                                                         <  2:137336/251‑1 (MQ=255)
                                                                                                                                                                                  ggcggGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCAc                                                         >  1:573246/1‑251 (MQ=255)
                                                                                                                                                                                        cACGTAGTTATTCAACTT‑GAAGTA‑CGCTGGAGCTATACCACAAATCAAATAAACGGCAAGAGCGAGTAATACGCACTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTATAATTGGTAACA‑GCAG‑TTTACCTGCTGTTATGTCTTCAACCTACCATCTAGGTGGCGTTGGCACTGGCAATATAATGCTGTTATTCATTGTGTAATTTAACGCGGTCGACGGCGTCTTTATTGCCTGCACGCaaac                                                   >  2:263756/1‑249 (MQ=255)
                                                                                                                                                                                             agtaAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGct                                               <  1:417528/250‑1 (MQ=255)
                                                                                                                                                                                                    cAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACtagt                                         >  2:187317/1‑249 (MQ=255)
                                                                                                                                                                                                                   gCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTc                          >  1:168255/1‑251 (MQ=255)
                                                                                                                                                                                                                    cTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCt                                                                                                                                                                 >  2:503320/1‑115 (MQ=255)
                                                                                                                                                                                                                    cTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCt                                                                                                                                                                 <  1:503320/115‑1 (MQ=255)
                                                                                                                                                                                                                            ataCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCAACGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGCGGCACTGGCACTATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGTCGACTGCGACTTTATTGCCGGCACACAACAGGCTACTAGTAGCAGTAAACCGCGCGCCcctgcc                 >  2:2522/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ccAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGgc            >  1:145981/1‑251 (MQ=255)
                                                                                                                                                                                                                                      taaaataaGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGgcgc          >  1:149140/1‑248 (MQ=255)
                                                                                                                                                                                                                                          ataaGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCCGTAAAACGCTCGCCCCTGACTGGGCGCAGCTAAt   >  2:467655/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ataaGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAAt   <  1:157205/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATc  >  1:419796/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         
TCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTT‑GAAGTA‑CGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACA‑GCAG‑TTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATC  >  NC_000913/4003896‑4004379

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