Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,5670CT39.7% 60.5 / 47.9 59intergenic (+174/+541)eco/mqoserine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (22/13);  new base T (12/11);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAT  >  NC_000913/2304324‑2304816
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
aGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCGCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGccc                                                                                                                                                                                                                                                       <  1:102852/251‑1 (MQ=17)
   caaggtggATGTGAAGTCCCGCGTCGTGAAGGCGGAAGACAAAATTGGCAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCCGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGcccaag                                                                                                                                                                                                                                                    <  2:485334/244‑4 (MQ=11)
       gTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACCTCTAACGCCCAAGCCgg                                                                                                                                                                                                                                                >  1:423910/1‑251 (MQ=17)
        tAGATGCGAAGTACCGCGTCTGGAAGGCGGCAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGGTACCCCGCATCCGACATCTAACGcccaaaccggg                                                                                                                                                                                                                                               >  2:676881/1‑243 (MQ=14)
              ggaattacggcagtgggAAGGCGGAGGAGAAAATTGAAAAG‑GCGGTAGTTCGCTAAACTGCCGTAAAGTGCGGCCCCCCGTAGGTCAGACAGGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTg                                                                                                                                                                                                                                         <  2:913772/236‑1 (MQ=11)
               gAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTGGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                                                                                                                        <  1:289909/251‑1 (MQ=14)
                aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCCCATCCGACATCTAACGCCCCAGCCGGTTGCCtcat                                                                                                                                                                                                                                       >  1:222546/1‑248 (MQ=14)
                 aGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCCTCCGACATCTAAAGCCCAAACCGGGTGccagctt                                                                                                                                                                                                                                      >  1:439365/1‑246 (MQ=14)
                     ccGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACACCGCATCCGCCATCTAACGCCCAAGCCGGGTGCCTGATGCGAc                                                                                                                                                                                                                                  >  1:286099/1‑251 (MQ=14)
                      ggcgTCTGGAAGGCGGAAAAGAAAATTG‑CAACCGCGTTATTTCGCTAAACTGCCGTGAAGGGCGGCCCCCCGGAGGTCAGAAAGGGCGGCCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                                                                                                                                                                 <  2:794707/250‑1 (MQ=11)
                         gTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg                                                                                                                                                                                                                              <  1:59361/251‑1 (MQ=14)
                           cTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCCTCCGACCTCTAACGCCCCACCCGGTTGCCCGATGCGACGCTGgc                                                                                                                                                                                                                            >  2:651653/1‑251 (MQ=14)
                             ggAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGACCGATGCGACGCCGgcgc                                                                                                                                                                                                                          >  2:102807/1‑251 (MQ=14)
                                aGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGACGGTCACGCCGCATCCGACATCCCACGCCCGAGCCGGTTGCCTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGGTACCCCGCATCCGGCATCTAACACCCAAGCCCGGTGCCTGCTTCGACGCAGgccggtc                                                                                                                                                                                                                       >  2:636937/1‑246 (MQ=11)
                                 ggCGGAAGAGAAAATTGACACC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                                                                                                                                                      <  1:1220214/251‑1 (MQ=11)
                                    ggAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGGTGCGACGCTGGCGCGTCtt                                                                                                                                                                                                                     >  2:410071/1‑249 (MQ=11)
                                    ggAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCtt                                                                                                                                                                                                                     <  1:410071/249‑1 (MQ=11)
                                      aaGAGAAAATTGACAAC‑GCAGTAGTTCGCTAAACGGCCGTGAAGTGCGGCCCCCCGTAGGCCAGACAGGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                                                                                                                                                 <  1:1064165/251‑1 (MQ=11)
                                         agaAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCCACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAgac                                                                                                                                                                                                              >  1:1030064/1‑249 (MQ=11)
                                           aaaaTTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCACGTCTTATCAGGCCt                                                                                                                                                                                                            <  1:203955/251‑1 (MQ=11)
                                              aTTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCCTCCGACCTCTAACg                                                                                                                                                                                                                                                          >  1:488436/1‑202 (MQ=11)
                                              aTTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAACCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                                                                                                                                                      >  2:76724/1‑238 (MQ=11)
                                              aTTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                                                                                                                                                      <  1:76724/238‑1 (MQ=11)
                                                                                                                                                                          cggtagttcgctaaactgccgtgAAGAGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAACCGGGTGCCTGATGCGGCCCTGGCGCGGCTTAATAGGCCTACACCGCtggga                                                                              >  1:953255/21‑247 (MQ=11)
                                                                                                                                                                              agttcgctaaactgccgtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa                                                                                                  <  1:1138903/210‑1 (MQ=255)
                                                                                                                                                                                 tcgctaacctgccgtgAAGTGCGGCCCCCCGTGGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTc                                                                      <  2:759812/237‑1 (MQ=255)
                                                                                                                                                                                 gcgctaatctgccgtgAAGTGCGCCACCCCGGAGGTC‑GACCAAGTCGGTCACGCCGCACCTGACCTCCACCGCCCGAGCCGGTTGCGGGAGGCGACGCTGGCGCGTCTGATCAGGGCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTcc                                                                     <  2:286099/238‑1 (MQ=11)
                                                                                                                                                                                  cgctaaactgccgtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACCTCT‑AACGCCCAAGCCGGCTGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTG‑TGAAGTGCTcc                                                                     >  1:1030686/13‑250 (MQ=255)
                                                                                                                                                                                   gctaaactgccgtgAAGTGCGGCCCCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtg                                                                                 <  1:1072084/227‑1 (MQ=255)
                                                                                                                                                                                      aaactgccgtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑cccc                                                               <  1:51688/243‑1 (MQ=255)
                                                                                                                                                                                       aactgccgtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACATCT‑AACGCCCAACCCGGTTGCCTGATGCCACGCTGGCGCGTCTTATCAGGCCCTCACCGCTG‑TGAAgaggcccaccccc                                                               >  2:128721/8‑239 (MQ=14)
                                                                                                                                                                                       aactgccgtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑ccccc                                                              >  1:299610/8‑250 (MQ=255)
                                                                                                                                                                                          tgccgtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTACCTGATGCGACGCTGGCGCGTCTTATCAAGCCTACACCGCGG‑TGAAGTGCTCCA‑CCCCGt                                                             >  1:637755/5‑249 (MQ=255)
                                                                                                                                                                                            ccgtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTTTCAGGCCTACACCGCTG‑TGAAGGGCTCCC‑CCCCGTAGGt                                                         >  2:874627/3‑251 (MQ=255)
                                                                                                                                                                                              gtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACCTCT‑AACGCCCAAGCCg                                                                                                                               >  2:385461/1‑181 (MQ=255)
                                                                                                                                                                                              gtgAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCg                                                                                                                               <  1:385461/181‑1 (MQ=255)
                                                                                                                                                                                                  attGCGGCACCCCGTAGGCCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGt                                                                                                                             <  1:1102962/177‑1 (MQ=255)
                                                                                                                                                                                                  aGTGCGGCACCCCGTAGGTCAGCC‑ACGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTCTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATACGGCATCTTAAC‑CCCCAGCCCGTTGCCTGTTTCGACGCCGGCGCGTCATATCAGGCCTCAACCCCTG‑TGAAGTGATCCA‑ccctggatggtggact                                                   >  2:1215334/1‑238 (MQ=14)
                                                                                                                                                                                                  aGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGGTACGCCGCATCCGGCATCT‑TACGCCCCAACCGGt                                                                                                                             >  2:1102962/1‑179 (MQ=14)
                                                                                                                                                                                                  aGTGCGGCACCCCGTAGGCCAGAC‑AAGGCGGCCACGCAGCATCCATCATCCAACGCCCGTGCCGGTTGCCTGAGGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATa                                                   <  2:274029/251‑1 (MQ=14)
                                                                                                                                                                                                       ggctccccGTAGGTCAGCC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCg                                              <  1:638718/247‑1 (MQ=255)
                                                                                                                                                                                                       ggCACCCCGTCGGTAAGAC‑AAGGCGGTCACGCCGCATCCGACATCAAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCGTACGAAATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCATTTCAGGCCTGCACCCATG‑TTAAGGGCTCCAACCCCGTGGGGTGG‑TAANGCg                                              >  2:1155007/1‑251 (MQ=11)
                                                                                                                                                                                                       ggCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCg                                              >  1:667142/1‑251 (MQ=255)
                                                                                                                                                                                                          gatcggtgggTCAGCC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAGGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGtt                                           <  2:216166/243‑1 (MQ=255)
                                                                                                                                                                                                          aCCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCCGTTGCCTGGTGCGACGCTGGCGCGTCTTATAAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGGAGGTTGGATAAGGCgggt                                           >  2:521175/1‑249 (MQ=255)
                                                                                                                                                                                                          aCCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGGTACGCCGCGTGCGACAGCC‑AACGCCCCAGGCGGTTGCCTGATGCGACGCCGGCGGGTCTTATCAGGCCGACAACCCTG‑GGAAGTGGTCCA‑CCCCCTAGGtgaggtaaggggggg                                           >  2:484969/1‑237 (MQ=11)
                                                                                                                                                                                                             ccGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCCACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTcc                                                                                                                                                                                        >  2:569943/1‑110 (MQ=17)
                                                                                                                                                                                                             ccGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTcc                                                                                                                                                                                        <  1:569943/110‑1 (MQ=32)
                                                                                                                                                                                                               gTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATACGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCgggtagatg                                      >  1:216166/1‑244 (MQ=255)
                                                                                                                                                                                                                tAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGGTGGGTAAGGCgggaggtcgg                                     >  2:1192949/1‑243 (MQ=255)
                                                                                                                                                                                                                tAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGt                                            <  1:1192949/244‑1 (MQ=38)
                                                                                                                                                                                                                  ggTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCCGATGCGACGCTGGCGCGTCTTATCAGGCCCACACCCCTG‑GGAAGTGCTCCA‑CCCCGTAGGGTGGACAAGGCgggttagccccc                                   >  1:525660/1‑241 (MQ=11)
                                                                                                                                                                                                                    tCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAACCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑GGAAATGCTCCA‑CCCCGTCGGGTGGATAAAGCGGTTaagacgcagc                                 >  2:28461/1‑242 (MQ=11)
                                                                                                                                                                                                                     cAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACATCT‑AACGCCCAAGCCGGGTGCCTGATGCGGAGCTGGGGCGTCTTATCAAGCCTACAACGCTGGTGA‑GTGCTCCC‑CCCCTGAGGTTGGATAAGGAGGTTTCACCGCaacc                                >  2:1083239/1‑248 (MQ=11)
                                                                                                                                                                                                                       gAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                              <  1:737032/251‑1 (MQ=35)
                                                                                                                                                                                                                        aC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTAAGCCGCCTCCGAc                             >  1:1114744/1‑251 (MQ=25)
                                                                                                                                                                                                                             ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTgg                                                      >  1:1022514/1‑222 (MQ=33)
                                                                                                                                                                                                                                        cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGGTGGATAAGGGGGGGTCGCCGCAACCGGCACCCAACCTCCGAGc              >  2:914375/1‑251 (MQ=11)
                                                                                                                                                                                                                                         gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGGTTCGCCGCATCCGGCATCCCACGTCCGAGcc             >  1:700427/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ccctcctccactcAACGCCCGAGCCGGTTGCTCGATGCGACGCTGGGGCGTCTTATCAGGACTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCACCCGACATCT‑AACTCCCAAGCCGGTTGCCTTATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGcc             <  2:953255/239‑1 (MQ=11)
                                                                                                                                                                                                                                                    ccAACGCCCGGGCCGGTTGCCGGATGCTACGCTGGCGCGTCTTATCAGTCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg    <  2:700427/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAt  <  1:914375/251‑1 (MQ=37)
                                                                                                                                                                                                                                                     cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGGTGGATAAGGCGGGTAAGCCCCATCCCACATCCAAACTCCCAGCCGGGTGcccgg   >  2:210737/1‑246 (MQ=11)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
AGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAT  >  NC_000913/2304324‑2304816

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