Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,621,1180AC40.0% 47.0 / 12.1 47V7G (GTG→GGG) eamAcysteine/O‑acetylserine exporter EamA
Reads supporting (aligned to +/- strand):  ref base A (24/3);  new base C (1/17);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.67e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGGACAGACAAATGAACTTATCCCTACGACGCTCTAC  >  NC_000913/1620879‑1621357
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
cagcagGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTA‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAACGCCAACCCCC‑C‑ATCcttt                                                                                                                                                                                                                                      >  2:72046/1‑251 (MQ=255)
 agcagGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGATAATCCATACCCCACAAGCAAATTCAGTGGGACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAACGACCAGCATAAAGCGCACACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAACATTTAGCCCCCATA‑CG‑ACCACTACCAGGAGCGCCAACACCC‑C‑ATCTTTTc                                                                                                                                                                                                                                     >  1:461013/1‑251 (MQ=255)
  gcagGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTcg                                                                                                                                                                                                                                    >  1:421297/1‑251 (MQ=255)
   cagGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTcgc                                                                                                                                                                                                                                   <  1:45140/251‑1 (MQ=255)
     ggCATACCGAAGTTAATGGCACAAAAAAGAAAAGCCAACTGCGCAAAACAGATGGTTAATCAATACCCCAGCAGCCAATTCAGTGGTACTTTC‑GGTCGTGCGACAGAAAAGATAGCCGGAAAAGCGACAAGCATCAAGCGCAAACCGGCCAGCAT‑CAGCGGAGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCCCTACCAGAAGCGCCAACACCC‑A‑ATCTTTTcgcg                                                                                                                                                                                                                                  >  2:153101/1‑250 (MQ=255)
      gCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAACCTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCTTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CATCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGAc                                                                                                                                                                                                                                <  1:518057/251‑1 (MQ=255)
        aTACCGACGTGAATGGCACAAAAAAGAAAAGCAACATGCCCAAAACTGATGGGTAAGCCATACCCCAGCAGCAAATTCAGTGGTACTTGC‑GGCCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCAGAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTGGATGACCCCAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACTCCC‑C‑ATCTTTTCGCGatat                                                                                                                                                                                                                              <  2:192113/251‑4 (MQ=255)
             gAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTAc                                                                                                                                                                                                                         >  2:132670/1‑251 (MQ=255)
                gTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACCAAATTTAGCCCCCCTA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGc                                                                                                                                                                                                                      >  2:5043/1‑251 (MQ=255)
                   aaTGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATCAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCct                                                                                                                                                                                                                    <  1:317643/250‑1 (MQ=255)
                     tGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACACCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAACAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGAGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTTGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCCACCGCCTCTg                                                                                                                                                                                                                 >  2:426346/1‑251 (MQ=255)
                          caAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAAttt                                                                                                                                                                                                            >  1:15417/1‑251 (MQ=255)
                               aagaaaagCAAAATGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTAGCGACATTCTACCGCCTCTGACTTGCATCt                                                                                                                                                                                                       >  2:91475/1‑251 (MQ=255)
                                   aaagCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTAACAGGAGCGCCAAAACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTATGAATTTCAACttttt                                                                                                                                                                                                   >  1:233342/1‑250 (MQ=255)
                                    aagCAAACTGCGCAAAACTGATGGTTAACCCATACCCCAGCAGCAAATTCAGTGGTCCTTGC‑GGTCGTGCGAAAAAAAAGATAGCCGGAAAAGCGACCAGAAAAAAGCGAAAACCGGCCAGCAC‑CAGCGGTGGCATGTTATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑CCCCCTCCCAGTGGCGCCACCCCCC‑C‑ATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                                                                                                                  <  1:175762/251‑1 (MQ=255)
                                     agCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                                                                                                                  >  1:113053/1‑250 (MQ=255)
                                             gcgcAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                                                                                                                          >  1:369397/1‑250 (MQ=255)
                                             gcgcAAAACTGAGGGTTAATCCATACCCAACCAGCAAATTCAGTGGTACTTTC‑GGTCGGGCGAAAAAAAAGATAGCCGAAAAAGCGCCCAGCAAAAAGCGCAACCCGGCCAGAAT‑CAGCGGTGGCATGTTAGGAAGCCCCCCTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTCCAAGTAGCGCCAACCCCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                                                                                                                         <  1:194041/251‑1 (MQ=255)
                                              cgcAAAACTGATAGTTAATCCATACCCAAGCAGCAAATTCAGTGGTACTTTC‑GGCCGGGCGAAAAAAAAAATAGCCGAAAAAGCGACCAGCATAAAGCGCAAACCGGCCACCAT‑CAGCGGTGGCATGTTATGAACCCCCACTTTGATGACCAAAAAATTTAGCCCCCATA‑CG‑ACCCCTACCAGTAGCGCCAACCCCC‑C‑ATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCaa                                                                                                                                                                                        <  1:86577/251‑1 (MQ=255)
                                                      tGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                                                                                                                >  1:52526/1‑251 (MQ=255)
                                                               acccATACCCCAGCAGCAATTTCAGGGGT‑CTTTCCGGTGGGGCGACAAAAAAGATAGCCGAAAAAGCGACCAGCAAAAAGCGAAACCCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCCCAAAATTTAGCCCCCATC‑CG‑CCCACTCCCAGTAGCGCCAACCCCC‑C‑TTCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                                                                                                                       <  1:8656/249‑1 (MQ=255)
                                                                gtccagaggtgtgctctcccgatcgggTACGTTG‑GGTCGTGCGACAAAAAAGATAGCCGGAAAACCGACCAGAATAAAGCGCAACCCGGCCAGCAC‑CAGCGGTGGCTTGTTATGAGGCCCCTCTTTGATGACCACAAAATTTAGCCCCCATC‑CG‑ACCCCTACCAGTAGCGCCAACCCCC‑C‑TTCTTTGCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                                                                                                                      <  1:118968/226‑1 (MQ=255)
                                                                     aCCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                                                                                                                 >  1:33808/1‑251 (MQ=255)
                                                                                        agggaCTTTC‑GGTCGTGTGACAAAAAAGATAGCCGGAAACCCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGGGGCAAGTTATGAAGCCCCCCTTTGAAGACCACAAAAGTTAGCCCCCATA‑CG‑ACCACTACCAGAAACGCCAACACCC‑C‑TCCTTTTCGCGACATTCAACCGCCTCTGAATTTCATCTTTTGTAATCAATCAACTAAGCTGAAATTTCTTTTCTTTAACCATTGATTCGTTattc                                                                                                                                              >  2:531667/5‑248 (MQ=255)
                                                                                         ggaaCTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCCACTTAGCTGAATTTACttt                                                                                                                                                                      >  2:118968/4‑226 (MQ=255)
                                                                                               tttgggTCGTGGCACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAC‑CAGCGGTGGCATTTTATGAACCCCCACTTTGATGCCCACAAAATTTACCCCCCAAA‑CG‑CCCCCTCCCAGTCGCCCCACCCCCC‑C‑ACCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCgg                                                                                                                                         <  1:390192/246‑1 (MQ=255)
                                                                                                cc‑GGTCGGGCGACAAAAAAGAAAGCCGAAAAGGCAACCAGCATAAAGCGCAACCCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCTCTTTGATGACCACAAAATTTACCCCCAATC‑CG‑ACCACTCCCAGTAGCGCCAACCCCC‑C‑ATCTTTCCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                                                                                                                       <  1:140799/249‑1 (MQ=255)
                                                                                                         gCGAAAAAAAAGATAGCCGGAAAAGCACCCAGCAAAAAGCGAACACCGGCCCG‑ATGCAGGGGTGGCTTGTATTGAACCCCCATTTGGATGACCAAAAAATTAACCCCCCATA‑CG‑CCCCCTACCATTAGCCCCAACTCCC‑C‑ATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCa                                                                                                                              <  1:338662/251‑1 (MQ=255)
                                                                                                             aaaaaaaGATAGCCGGAAAAGCGACCAGCAAAAAGCGCAAACCGGCCAGCAC‑CAGCGGGGGCATGTTATGAAGCCCCCCTTGGATGACCACAAAATTTAGCCCCCATC‑CG‑ACCCCTACCAGTAGCGCCAACCCCC‑C‑ATTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTaa                                                                                                                          <  1:281626/250‑1 (MQ=255)
                                                                                                                       ggCCGGAAAAGCGCCCAGCAAAAAGCGCAAACCGCCCATCAC‑CAGCGGTGGCATGTTATGAACCCCCCTTTTGATGCCCAAAAAATTTACCCCCCATA‑CG‑ACCACTCCCAGTAGCGCCAACCCCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                                                                                                                 <  2:15417/249‑1 (MQ=255)
                                                                                                                        gCCGGAAAAGCGACCAGCAAAAGGCGAAAACCGGCCAGCAC‑CAGCGGTGGCATTTTATGAAGCCCCACTTGGATGCCCACAAAATTTAGCCCCCATA‑CG‑ACCACTACAATTAGCTCCAACCCCC‑C‑ATCTTTTCGCGACATCCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                                                                                                               <  1:511535/251‑1 (MQ=255)
                                                                                                                              aaaGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt                                                                                                         >  1:472046/1‑251 (MQ=255)
                                                                                                                                aGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAACTTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCACATAAGTCGCtat                                                                                                       >  2:316230/1‑251 (MQ=255)
                                                                                                                                       aGCATCAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTGGCCCCCTTA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTTTACGCAATCCACTTAGCTGAATTTACTTTTCTTTAACAGTTGTTTCGTTTTACGCCGGTTACGTCGGCCTTAATGCGCAAATTAGTCGCtttccttccg                                                                                                >  2:324499/1‑242 (MQ=255)
                                                                                                                                            aaaGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtttt                                                                                           >  1:329291/1‑251 (MQ=255)
                                                                                                                                                   aaaacaggccataagcAGCGGGGGCATGTTATGAAGCCCCCCTTACGTCACCACAAAAAT‑CGCCCCCCTAACG‑TCCTGTACCAGTA‑CGGCAACCCCCTC‑ATCTTTTCCAGACATTTTTCCGCCTCTGAATTTCATCTTGTGTAAGCAATTAACTTAGCTGAAGTTCCTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGGCAt                                                                                     <  1:35881/236‑1 (MQ=255)
                                                                                                                                                     cccGGCCAGCAC‑CAGCGGTGGCATTTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑CCCACTACCAGTAGCGCCAACCCCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCt                                                                                  <  1:111459/250‑1 (MQ=255)
                                                                                                                                                                   ggCGGTGGCATGTTGTGAAGCCCCCTTTTGAGGACCACAAAATTAACCCCCCA‑A‑CGCACCCCCACCAGAAGCGCCAACCCCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTTAATTTGCTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCAATGCTATTTCTTTACAtc                                                                     <  2:446940/250‑1 (MQ=255)
                                                                                                                                                                    gCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGACAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCCAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACAtct                                                                    >  2:56797/1‑251 (MQ=255)
                                                                                                                                                                     cGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACAtctc                                                                   >  1:121643/1‑251 (MQ=255)
                                                                                                                                                                      ggTGGCATTTTATGAACCCCCCCTTTGATGACCAAAAAATTTACCCCCCATA‑CG‑ACCACTACCAGAACCGCCAACCCCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATctct                                                                  <  1:132670/251‑1 (MQ=255)
                                                                                                                                                                               tgaggaacccacgcttggaggaCCACAAAATGTAGCCCCTATA‑CG‑CCCACACCCAGTAGCGC‑AACCCCC‑CGTTCTTTCCGAGATTGTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTGTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTACGGATTTTTGCCATGCTATTTCTTTACAACTCTaaaacaaa                                                          <  2:152281/230‑1 (MQ=255)
                                                                                                                                                                                 tataaGCCCCACTGTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTCCCAGTAGCCCCAACCCCC‑C‑TTCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTaaaacaaaac                                                        <  1:128308/247‑1 (MQ=255)
                                                                                                                                                                                              ttgaggccaaaaaaaTTTACCCCCCATC‑CG‑CCCACTACCAGTAGCGCCAACCCCC‑C‑ATCTTTCCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCa                                           <  1:6714/239‑1 (MQ=255)
                                                                                                                                                                                                gatgaCCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTATGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATCAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTg                                        >  2:218064/1‑251 (MQ=255)
                                                                                                                                                                                                   gaCCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTATTCAATCAACTTAGCTGAATTTACTTTTATTTAACAGTTGATTCGTTCGTCGCCGGTTACGAGGGCATTAATGCCCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAACCAAAACATATCGAACCGCATTGCCg                                     >  2:78829/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGGACAGACAAATGAACTTATCCCTACGACGCTCTAc  >  2:409364/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
CAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTC‑GGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCAT‑CAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATA‑CG‑ACCACTACCAGTAGCGCCAACACCC‑C‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGGACAGACAAATGAACTTATCCCTACGACGCTCTAC  >  NC_000913/1620879‑1621357

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