Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,8740GA35.7% 75.4 / 24.9 56A56T (GCT→ACT) ychSputative uncharacterized protein YchS
Reads supporting (aligned to +/- strand):  ref base G (16/20);  new base A (13/7);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGC  >  NC_000913/1286580‑1287132
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                  
gTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAGCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTAGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCACCGCTTGGGCAGAATATTTAATTGCGGATTCGTTTGGGAA‑TTCCGGGACtttt                                                                                                                                                                                                                                                               >  2:379559/1‑300 (MQ=255)
gTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGCAGCTCCTCACCCTCCCCCAACTCAGCAGAATAGATGATTGCGGCTTCGTTTGAGAA‑TTCCGGGGCtttt                                                                                                                                                                                                                                                               >  1:610482/1‑300 (MQ=37)
gTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCtttt                                                                                                                                                                                                                                                               >  1:168805/1‑300 (MQ=255)
gTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGAG‑TTCAGGGGCtttt                                                                                                                                                                                                                                                               >  1:656575/1‑300 (MQ=255)
    gggtttaGTTGTCGGATGCGCTTCGCTTATCCGACCTACATGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGaa                                                                                                                                                                                                                                                            <  2:644783/293‑1 (MQ=255)
      ggttaGTTTTCGTATGCGCTTCGCTTATCCGACCTACAGGTGAGGATATTGTAGTCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGTTTTTGATTATTTAAGGAATCGATGGTGTTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGT‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGaa                                                                                                                                                                                                                                                            <  2:311370/293‑1 (MQ=25)
           gTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTT‑GGGAAGTTCAGGGACTTTTGAAAGTGAt                                                                                                                                                                                                                                                      >  1:476299/1‑298 (MQ=18)
                   aTGCGCTTCGCTTATCCGACCTACAGGTGAGGATAGTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGTAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGTCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                                                                                                                                                                                                                                            <  2:97116/300‑1 (MQ=255)
                          tcgcttATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCGACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTT‑GGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                                                                                                                                                                                                     >  2:564462/1‑300 (MQ=17)
                          tcgcttATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTT‑GGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                                                                                                                                                                                                     >  2:451924/1‑300 (MQ=17)
                               tatccgatcggcgggtgaggaTATTGTGGGCCCGTTAGACTTGGTGCCATCGGGCAAAATTAGCTCGGAGTAATGGGTTTTGATTATTTAAGGGAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                                                                                                                                                <  2:2990/284‑1 (MQ=18)
                               tATCCTACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGTAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCTG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                                                                                                                                                <  2:330893/300‑1 (MQ=255)
                               tATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGTTGGAAAG‑TTCAGGGACTTTTGAAAGTGAtggtgttgtata                                                                                                                                                                                                                                           >  2:134713/1‑282 (MQ=18)
                               tATCCGACCGACAGGGGAGGATATTGTATGCCCGTTAAACGTGGTGCCATCGGGCAAAATTAACTCGGAGTAATAGGTTTTGATTATTTGAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                                                                                                                                                <  2:482383/300‑1 (MQ=14)
                                    gACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                           <  1:17188/300‑1 (MQ=255)
                                     agcgaCAGGGGAGGATATTGTAGGCCCGGTAAACTTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                           <  2:705395/295‑1 (MQ=255)
                                       cTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGNAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTNGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAA‑TTCAGGGACTTTTGAAAgctatctcttctac                                                                                                                                                                                                                                             >  2:523043/1‑266 (MQ=18)
                                       cTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCACGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCt                                                                                                                                                                                                                          >  1:705446/1‑298 (MQ=255)
                                       cTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgcct                                                                                                                                                                                                                         >  1:705395/1‑297 (MQ=255)
                                             tggaggaTATGGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCGCGGAGGAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                                                                                                                                                                                                                      <  2:457398/295‑1 (MQ=255)
                                             gggAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTT‑GGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcgaaccttcgaagt                                                                                                                                                                                                                      <  1:451924/296‑15 (MQ=14)
                                              ggaggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTAATTGCGGATTCGTTGGGGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcgccctttcga                                                                                                                                                                                                                         >  2:49488/1‑281 (MQ=17)
                                              ggaggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt                                                                                                                                                                                                                   <  1:545612/298‑1 (MQ=255)
                                                 ggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCANAATATTTGATTGCGGATTCGTT‑GGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcgcaccttcgacgtccctgac                                                                                                                                                                                                               >  2:543491/1‑278 (MQ=17)
                                                          gggCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAATCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTc                                                                                                                                                                                                     <  1:120455/299‑1 (MQ=255)
                                                          gggCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTc                                                                                                                                                                                                     <  1:120433/299‑1 (MQ=255)
                                                           ggCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTcc                                                                                                                                                                                                    >  2:120455/1‑299 (MQ=255)
                                                           ggCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTAGCCCTTCGGGTcc                                                                                                                                                                                                    >  2:120433/1‑299 (MQ=255)
                                                                ggTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAAGACGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCACGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTAATCCTTCCCCGCTTGGGCGGAATATTCGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCtttt                                                                                                                                                                                                                                                               >  1:260349/1‑236 (MQ=21)
                                                                    aaCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAG‑TTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                                                              >  1:170403/1‑297 (MQ=255)
                                                                           tgccagagggtaaaattgactcgtaggggtaGGTTTTTGTTATTTAAGGAAGTGATGGTGGTGGGGGAAGGACTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGTAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                                                                                                                                    <  2:31646/270‑1 (MQ=14)
                                                                                   ggCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                                                                                                                            <  1:430937/300‑1 (MQ=255)
                                                                                         aTTAGCTCGGAGTAACAGTTTTTGATTATTTAAGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                                                          <  2:703874/296‑1 (MQ=255)
                                                                                                                             actgatgttgtttgggtaatggttACTCAGCTCTTCGCGCTCCGCTCTTCGTTGT‑GTTGCCTGCGTCAACGCCCTCCCGCTTGTGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATGTGATTGCGGATTCGTTTGAGAA‑TTCCGGGACTTTTGAAAGTGATGGTGTTGGGGGAACGATTACTCAGCGCTGCACGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGATGCTCGAGCCAAACCCTAGTCGAAGCTTCTCACCCTTc                                                                                                                                   <  2:166069/277‑1 (MQ=14)
                                                                                                                                  tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                                             <  1:56425/300‑1 (MQ=255)
                                                                                                                                           gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGTAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                                                                                                                                    <  1:190586/236‑1 (MQ=255)
                                                                                                                                           gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                                                                                                                                    >  2:190586/1‑236 (MQ=9)
                                                                                                                                                ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                                               >  1:342097/1‑300 (MQ=255)
                                                                                                                                                ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTCGGTCGAAGc                                                                                                                                               >  2:576861/1‑268 (MQ=255)
                                                                                                                                                ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGTTTGGGAG‑TTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcgccaac                                                                                                                                                                                                                             >  2:576868/1‑183 (MQ=11)
                                                                                                                                                   ttACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                                            <  1:379272/300‑1 (MQ=255)
                                                                                                                                                     aCTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTACACGCCTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAACTCTGCGCCCCTCGccc                                                                                                                                                                                                             >  2:232222/1‑201 (MQ=14)
                                                                                                                                                                            gtctg‑gTCGTTGCCGGCGGTAACGCTCTTTCGCGTGTGTTCAAGTCGAACTTTGGTCGAAGCTTCTCACCCTTCCCCGCTTGTGCAGAATATGTAATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTGTGAA‑GTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                                                  <  2:71721/296‑1 (MQ=17)
                                                                                                                                                                             tCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCATGGa                                                                                  <  1:21126/300‑1 (MQ=255)
                                                                                                                                                                                g‑gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                                                                                                                                                                                                                    >  1:464128/1‑151 (MQ=255)
                                                                                                                                                                                        gCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTAATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCGCCCTTCGGGCCGTTGCCTGTGGCAACGCTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGATTCTCATCCTTCCCCGccgagttcaaactttcaacctct                                                                                                        >  1:224085/1‑246 (MQ=18)
                                                                                                                                                                                                       ctctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTAATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                              >  1:170770/1‑295 (MQ=255)
                                                                                                                                                                                                              tGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCATGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                   <  1:346751/299‑1 (MQ=255)
                                                                                                                                                                                                                      aaGTCGAACCTTGGTCGAAGCTTCTCACCCTTCCCCGCTTGGTCAGAATATGTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgc                                           <  1:401042/299‑1 (MQ=255)
                                                                                                                                                                                                                       aGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAATTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                         <  1:322797/300‑1 (MQ=255)
                                                                                                                                                                                                                       aGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTAATCCTTCCCCGCATGGGCAGGATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                         >  1:216227/1‑300 (MQ=255)
                                                                                                                                                                                                                                  ggTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGANTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGTGCTGCGCGCTTc                               >  1:499847/1‑299 (MQ=255)
                                                                                                                                                                                                                                          cTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCACCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                                                            >  1:39820/1‑198 (MQ=255)
                                                                                                                                                                                                                                                   cTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc             >  1:492431/1‑300 (MQ=255)
                                                                                                                                                                                                                                                    ttCCCCGCTGGGGCAGAAGGTTTGATTGCGGATTCGTGTGAGAA‑TTTCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTTTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                                                              <  2:185350/202‑1 (MQ=14)
                                                                                                                                                                                                                                                              gggCATAATATGTGATTGCGTATTCGTTTTAGAA‑TTTCGGGGTTTGTGAAATTGATGGTGGTGGGGGAAGGACTACTCATCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc  <  1:530024/300‑1 (MQ=17)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                  
GTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG‑GTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGC  >  NC_000913/1286580‑1287132

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