Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,244,2120CT25.4% 100.6 / 10.9 67G101G (GGC→GGT) yeiGS‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base C (35/12);  new base T (0/16);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 8.30e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACACGACTCCAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTGGGGAATCAAAGCGT  >  NC_000913/2243962‑2244459
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                          
gACACGACTCCAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACCCTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGc                                                                                                                                                                                                                                                            <  1:530980/251‑1 (MQ=255)
  cacGACTCCAGTACCTTAAACTGCCCGATGACGTTCAGTTGCTTTCTCCCTCCACCTCGTGAGCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCtt                                                                                                                                                                                                                                                          <  2:328478/251‑1 (MQ=255)
         ccAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCtt                                                                                                                                                                                                                                                   <  1:200146/251‑1 (MQ=255)
              cccTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGc                                                                                                                                                                                                                                             <  2:386847/250‑1 (MQ=255)
               ccTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGcc                                                                                                                                                                                                                                            <  1:462314/251‑1 (MQ=255)
                   aaaCTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGc                                                                                                                                                                                                                                        <  1:42028/251‑1 (MQ=255)
                       tGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGGG‑GTGCCCGGC‑ACCACCCCGCGCGGCGGAAA‑GGGTGCCAACG‑ACGAAGGCTACGATTTAGGCCCGGGCGCAGGCctttttttttattccacgccacc                                                                                                                                                                                                                                     >  2:599214/1‑232 (MQ=255)
                       tGCCCGATGAAGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑ATGTACCGGGG‑GTGCCAGGC‑ACCCGCCCGCGCGGCGAAAACGGTT‑CCCACC‑ACGATGGCCACGCATTAGGCCCGGGCGCAGGCCTTTATCTT‑AATGCCACGCCAcc                                                                                                                                                                                                                                    >  2:479231/1‑251 (MQ=255)
                                 cGTTCCGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑G‑TTGTAATGGTG‑ATGCCAGAC‑ACCCGACCCCCCGGCGGAAA‑GGTTGCCAACG‑ACGAAGGCTACGATTTAGGGCAGGGGGCAGGGTTTTATTTT‑AAATCCACGCCAccccccggggcg                                                                                                                                                                                                                          >  2:217657/1‑241 (MQ=255)
                                       gTATCTTTCTCCCCCCCCCTCGTGATCACACTCCGCACCCAGTGCTGTACGGGCTTTCCGGATTAACCTGCAAGGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACGGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGTCTACGATTTAGGCCAGGGCGCAGGTTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCAt                                                                                                                                                                                                                   <  1:379600/251‑1 (MQ=255)
                                                 acctccACATCGTGATCACACTCCGCAACCCGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGACAACTTCAC‑CACCCAGGCGGGGGGCCAGCGGGTAGCGGCGTGACTGGGG‑G‑GTGTACTGGTG‑ATTCCATGC‑ACCACCCCGCGCCGCGCAAC‑GGTTGCCCACA‑CCGCGGCCTACGGTGTAACCCAGGGCGCAGGCATTTATCTT‑AGTGTCACACCAccgcccacggcggcccccttcgcgatgg                                                                                                                                                                                                          >  2:90947/2‑228 (MQ=255)
                                                        ccTCGTGATCACACTCCGCCTCCAGTGCTTTAGTGGCTTTCCTGATAACCCTGCAAGGACGAGAACTTCCC‑CGCCAAGGCGGGGCCCCAGCGGGTGGCGGCGGAA‑TGGGGTA‑TTGTTCTGGTG‑ATGCCAGCC‑ACCAGCCCGCGGGGCAAAAA‑GGTTGCCAGCG‑ACGATGGCTACGATTTGGGCCGGGGCGCAGTTTTTTATTTT‑ATGGCCACACAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTa                                                                                                                                                                                                 <  2:630358/251‑1 (MQ=255)
                                                          cctggATCACCCTCCGCCACCAGTGCTGTACTGGCTTCCCGGATTAACCGGAAATGACGAGAACTTCAC‑CCCCAAGGCGGGGGCCCAGCGGGTAGCGGCGGACCTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGTTTTTATCTT‑AATGCCCCGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATc                                                                                                                                                                                               <  1:53880/247‑1 (MQ=255)
                                                              gATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTgcg                                                                                                                                                                                           >  1:336128/1‑251 (MQ=255)
                                                              gATCACACCCCGCCCCCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCCC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTgcg                                                                                                                                                                                           <  1:405253/251‑1 (MQ=255)
                                                               aTCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CAACAAGGCGGGCGCCCAGCAGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGGG‑AAGCCAGAC‑ACCAACCCCCGAGGCGACAA‑GGTTGCCAAAG‑ACGCCTGCTACGGTTTAGGCCCGGGGGACGGCTTTTCTCTT‑AAACCAACGCACCCCCCCTGGGC‑GGCGCATTACC‑GGATGTATGATTATCTgcgc                                                                                                                                                                                          >  2:505778/1‑251 (MQ=255)
                                                                  acacTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGCACTTGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCCCGCGGCGAAAA‑GGTTGCCAACG‑ACGCTGGCCACGAATTAAGCACGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCCCCCTGGGC‑GACGCATTACCCGGAT‑TATGAtttacttcccggt                                                                                                                                                                                       >  2:423533/1‑240 (MQ=255)
                                                                       ccccccccaGTGCTGTGCTGGCTTTCCGGATTACCCTACATGGCCGAGACCTCCCC‑CCCAAAGGGGGGCGCCCAGCGGGTAGCGGCGGAACGGGGG‑A‑TTGTACTGGGG‑AGGCCAGAC‑ACCAGCCCCCGCGGCGAAAA‑GGTTGCCAACG‑ACGTTGGCTACGATTTAGGCCGGGGCGCAGGTTTTTATATT‑AATGCCACCCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAAtt                                                                                                                                                                                  <  2:474657/245‑1 (MQ=255)
                                                                             ccaGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAAAACTCCAC‑CCCCAAGGCGGGGGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACGGGTG‑ATGCCAGAC‑ACCGGCCCGCGGGGCGAAAG‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCGGGGCGCAGGTTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGgc                                                                                                                                                                            <  1:599214/251‑1 (MQ=255)
                                                                               aGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGGGCCCAGCGGGTAGCGGCGGTACTGGCG‑ACTTTTA‑TGGTG‑ATACCAGAA‑ACCAGACCCCGCGGAGAAAA‑GGTTGACAACG‑ACGATGGCTACGATTTAGGGCAGGGCGCAGGACGTTATCTT‑AAAGCCCAGCCACCACCGCGGGCCGACG‑ATTACC‑GGATGTATGTTTAtaaccgccgtggatgaccggcgc                                                                                                                                                                          >  2:327604/1‑229 (MQ=255)
                                                                                          ggCTTCCCGGACTAACCTGCAATGACGAAACCTTCAC‑CCCCAAGGGGGGTGCCCAGCGGGGAGCGGCGGAACTGGGG‑A‑TTGTCCGGGTG‑ATGCCAGCC‑CCCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTCGGTTTTAGGCCAGGGCGCAGGCTTTTACCTT‑AATCCCACGAAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTc                                                                                                                                                                <  1:343245/250‑1 (MQ=255)
                                                                                            cTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑AAGATGGCTACGATTTAGGACAGGGCGCAGGCTTTTATTTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGCTGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGtagc                                                                                                                                                              >  1:124957/1‑247 (MQ=255)
                                                                                               tCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGGG‑ATGCCAGAC‑ACCCGCCCGCGCGGCGACAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCCGGCTTTTTTCTT‑AATGCCACGCCACCGCCGTGGGC‑GACGCATTACC‑GGATGGATGATTATCTGCGCGATGAAATACCGGCGCGGGTTccggcgcactt                                                                                                                                                          >  2:231734/1‑241 (MQ=255)
                                                                                                   gATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCTGGTGCCCAGCGGGTAGCGGCGGAACTTGGG‑A‑ATGTAATTGTG‑CTGCCAGAA‑ACCCGCCCGCGCGGCAAAAA‑AGTTGCCAACG‑ACGAAGGCTACGATTTTGGCCAGGGCGCAGGCTTTTTTCTTTAACCCCAC‑CCACCGCCCGGGGC‑GACGCACTACC‑GGAGGTATGACTTTCTCCGCGGTGAATTTCCCGCGCTgggtacggcgacacttaa                                                                                                                                                       >  2:603177/1‑234 (MQ=255)
                                                                                                   gATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGAACTGGTG‑ATGCCAGAC‑AACAGCCCGCGCGGCGGAAA‑GGTTGCCAACG‑ACGATGGCTACCATTTAGGCCCGGGCCCCGGCATTTATCTT‑AATGCCAAGCAACCGCCGTGGGG‑GCCGCAGTACC‑GGATGTTTGGTTATCTTCGCGCTGAATTACCGGGGCTGGTTCCGGCCCAATTTAAt                                                                                                                                                      >  2:144749/1‑251 (MQ=255)
                                                                                                       aaccgccAATGACGAGAACGT‑ACGCACCAGGGCGGGTGCCCAGGGGGTAGAGGGGGAACTGGGG‑A‑TTGGAC‑GGTGGATGCCAGAC‑CCCAGCCCGAGCGGAGAAAA‑GGTTGCCAACG‑ACGAGGGCTACGATTGAGGCCAGGGCGCGGGTTTTTGTCTA‑AATCCCACTCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTc                                                                                                                                                   <  2:69731/244‑1 (MQ=255)
                                                                                                                     agaACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTAATGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGGTTTAGGCCCGGGCGCAGGCTTTTATTTT‑AATGCCACGCAACCGCCGTGGGC‑GCCGCATTTCC‑CGGTTTATGTTTTTCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATTTCAGCGACCGCTTCGcc                                                                                                                                    >  2:420749/1‑251 (MQ=255)
                                                                                                                            tcccctccgaggggggTGCCCGGGGGGTAGCGGCGGAACTGGGG‑T‑TTGTACTGGTG‑ATGCAAGAC‑ACCAGCCCGCGTGGCGAAAG‑GGTTCCCAACG‑ACTGTGCCTACGATTTAGGCCAGGGCGAAGGTTTTTATTTT‑AATGCCACGCAGCCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTa                                                                                                                                <  2:130744/236‑1 (MQ=255)
                                                                                                                                caccaAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTCGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTc                                                                                                                          >  1:563969/1‑251 (MQ=255)
                                                                                                                                    aaggggggcccccggggggaaggggcggACCGGGGG‑T‑TTGTCCTGGTG‑AGGCCAGAC‑CCCAGCCCGGGCGGGGAAAA‑GGTTCCCAACG‑CCGAGGCCTGCAATTTAGGCCGGGGTGCAGGTTTTTATCTT‑AATCCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTc                                                                                                                      <  2:129930/228‑1 (MQ=255)
                                                                                                                                                ccggggggtggcggtgaagcggggaa‑TTGTGCTGGTC‑ATGCCAGCC‑ACCGGCCGGGGGGGGAAAAA‑GGTTGCCAACG‑AGGATGGCGACGATTTGGGCCAGGGCGCGGTCTTTTATCTT‑AATCCCACCCACCCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTc                                                                                                           <  2:242414/226‑1 (MQ=255)
                                                                                                                                                    ggggaGCGGCGGGACTTGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGATTTTATCTT‑ACTGCCACGCAACCGCCGTGGGC‑GACGCATTACA‑GGATTGATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCCGCGCCATTAGCGGTCACTCCATGGGGGGTCACGGt                                                                                                      >  2:587800/5‑251 (MQ=255)
                                                                                                                                                     ggTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGCAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTCGCGGTACCTCAATGGGTAGTCACgggg                                                                                                     >  2:19046/1‑249 (MQ=255)
                                                                                                                                                      gtggcggcgAAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑CCCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGAGGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATTTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTgc                                                                                                    <  1:479231/248‑1 (MQ=255)
                                                                                                                                                      gTAGCGGCGGAACGGGGG‑A‑TTGTACTGGGA‑ATGCCAGAC‑CCCACCCCGCGGGGCAAAAG‑GGTTGCCACCG‑ACGAGGGCTCGGATTGAGGCCGGGGCGCAGGTTTTTATCTT‑AATCCCACCCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTgc                                                                                                    <  1:616076/251‑1 (MQ=255)
                                                                                                                                                        aGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑CCCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCGCCGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATCCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTgcg                                                                                                   <  1:607606/250‑1 (MQ=255)
                                                                                                                                                         ggggcggacgggggg‑A‑TTGAACGGGTG‑ATCCCAGCC‑TCCAGCCCGCGTGGAGAAAT‑GGTTCCCACCA‑GCGATGGCGACGATTTAGGCAGGGGCGCAGGCTTTTATCTT‑AAGGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCCCGGTGCGCt                                                                                                 <  2:390622/240‑1 (MQ=255)
                                                                                                                                                         gcggggGAGCTGGGG‑G‑TTGGACTGGTG‑ATGCCAGCCGACC‑GCCCGCGCGGCAAAAT‑GGGTGCCAAGG‑ACGATGGCTTCGTTTTAGGCCAGGGCGCGGGTTTTTATTTT‑ATTCCCGCGCAGCCGCCGTGGGC‑GACGCATTACC‑GGATGTATGTTTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCt                                                                                                 <  2:438509/251‑1 (MQ=255)
                                                                                                                                                         gcggcgGAGCGGGGG‑T‑TTGTATTGGTG‑ATGCCAGAC‑ACCGCCCCGCGCGGCGAAAG‑GGTTGCCAACG‑CCGATGGCTACGATTTAGGCCAGGGCGCAGGTTTTTATCTT‑AATGCCACCCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCt                                                                                                 <  1:603957/251‑1 (MQ=255)
                                                                                                                                                         gcggagGAAATGGGG‑A‑TTGTACTGGTG‑ATGCAAGCC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCAAGGGCGCAGGCTTTTATATT‑AATGCCAAGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGGATTAACGGCGCTGGGTCAGGCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCATTGGGTGGTTACGGTgcgcg                                                                                                 >  2:580833/1‑250 (MQ=255)
                                                                                                                                                                               ccTGGTG‑ATGCCAGAC‑CCCAGCCCGCGGGGCAAAAA‑GGTTGCCAGCG‑ACGTTGGCTACTATTTGGGCCGGGGCCCGGGTTTTTATCTT‑AAGGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATc                                                                             <  1:217657/250‑1 (MQ=255)
                                                                                                                                                                                 tctggtgctgccgacccccacccaccgcgGTGAAAA‑GTTTCCCAACT‑ACGATTTCTCCGTTTTAGGCCAGGCCTCATGTTTTTACCTA‑AATGCCTCGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAAt                                                                              <  1:449112/225‑1 (MQ=255)
                                                                                                                                                                                 gggTG‑ATGCCAGCC‑CCCAGCCCGCGCGGCGAAAA‑GGTTGCAAACG‑ACGTTGGCTACGTTTTAGGCCAGGGCGCAGGTTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCg                                                                           <  1:423533/250‑1 (MQ=255)
                                                                                                                                                                                                            cggcgAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCt                                                  >  1:274282/1‑251 (MQ=255)
                                                                                                                                                                                                             ggcgAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCtt                                                 >  1:45250/1‑251 (MQ=255)
                                                                                                                                                                                                               cgAAAG‑GGTT‑CCAACGGACGATGGCTACGATTTAGGCGCGGGCCCAGGTTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTg                                               <  1:293362/251‑1 (MQ=255)
                                                                                                                                                                                                               cgAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGCGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGAGGTCACGGTGCGCTGATTATGGCGCTGCAAAATCCGGGTAAATACACCAGCGTTTTGGCCTTTg                                               >  2:371084/1‑251 (MQ=255)
                                                                                                                                                                                                                gAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGCTGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTgc                                              >  2:567777/1‑251 (MQ=255)
                                                                                                                                                                                                                gAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTgc                                              >  1:331121/1‑251 (MQ=255)
                                                                                                                                                                                                                aaaaa‑GGTTGCCAACG‑ACGAGGTCTACGATTTAGGCCGGGGCGCAGGTTTTTATCTT‑AATGCCACCCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTgc                                              <  1:367297/250‑1 (MQ=255)
                                                                                                                                                                                                                    a‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCCATGGGTGGTCACGGTGCGCTGCTTATGGCGCTGAAAAATCCGGGGAAATACACCAGCGTTTCGGCCTTTGCGcc                                           >  2:58070/1‑250 (MQ=255)
                                                                                                                                                                                                                       gTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGcccat                                        >  2:327281/1‑248 (MQ=255)
                                                                                                                                                                                                                       gTTGCCAAAG‑CCGCTGTCTACGATTTCGGCCAAGGCGCGGGCTTTTATCTT‑AATGCCACGCAAACGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGAGATGAATTACCGTCGCTGGTTCCGTCGCAATTTAATGTCAGAGACCGCTGCGCCATTGGCGGTCACTCAATGGGTGGTCACAGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCAGCCTTTgagacact                                        >  2:541789/1‑244 (MQ=255)
                                                                                                                                                                                                                         tGCCAACG‑ACGATGGCAACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACAGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTAACTCAATGGGTGGTCACGGTGCGATGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCACTtg                                      >  2:24234/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATtgt                                     >  1:349562/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaCG‑ACGCTGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGa                                   >  1:457250/1‑250 (MQ=255)
                                                                                                                                                                                                                                g‑acgaTGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATcc                               >  1:235222/1‑251 (MQ=255)
                                                                                                                                                                                                                                  acgaTGGCTACGATTTAGGCCAGAGAGCAGGCTTTGATCTT‑AATGCCACGCAACCGCCGGGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGTTCACTCAATGGGTGGTCACGGTGCGCTGATTATGACGCTGAAAAATCCGGGTAAATACACCATCGTTTCGGCCTCTGCGCCCATTGTGAATCCg                              >  2:628301/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ggtggcgaagATTAAGGCCAAGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGAGC‑GACGAATTACC‑TGATGTATGATTATCTGCGCGATGCATTACCTGCCCTGGTTCAGTCTCAATTTAATGTCAGCCACCGCTGCGCCATTAACGGTCACTCAATGGGTGGTCCCGGAGCGCTGATTATAGCCCTGAAAAATCCGGGTCAATACACCATCTCTTCTGCCTTTGCACAACTTGTGAGTCCGtt                            >  2:202045/10‑250 (MQ=255)
                                                                                                                                                                                                                                     aTGGCTACGATTTAGGCCAGGGCGCTGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCCATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAACGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTACGGCCTTTGCGCAAATTGTGAATCAGTGc                           >  2:333400/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGc                           >  1:414887/1‑250 (MQ=255)
                                                                                                                                                                                                                                        gCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGc                        >  2:367798/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cTACGAATTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCg                       >  1:404189/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGt                      >  1:412186/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aTTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCg                  >  2:198913/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tttAGGCCAGGGCGCAGGAATTTATCTT‑AATGCCAAGCAACCGACGTGGGA‑GACGCATTACC‑GGATGTATGATTATCTGCGAGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGATGCGCCATTAGCGGGCACTCAATGGGGGGACACGGTGCGCTGGTTATGGCGCTGAAAAATCCGGGTAAATACAACAGCGTTTCGGCCATTGCGCCAATTGGGAATCCGTGCCGCGTCCCGt                 >  2:472298/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     ccAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTGGGGa            >  1:429042/1‑250 (MQ=255)
                                                                                                                                                                                                                                                              gggCTTTTATCTT‑AATGCCCCGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTGGGGAATCAAAGCGt  <  1:105257/250‑1 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                          
GACACGACTCCAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCAC‑CACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGG‑A‑TTGTACTGGTG‑ATGCCAGAC‑ACCAGCCCGCGCGGCGAAAA‑GGTTGCCAACG‑ACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTT‑AATGCCACGCAACCGCCGTGGGC‑GACGCATTACC‑GGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTGGGGAATCAAAGCGT  >  NC_000913/2243962‑2244459

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