Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,5370AG54.5% ‑3.6 / 78.4 67intergenic (+144/+571)eco/mqoserine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (19/11);  new base G (18/18);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCC  >  NC_000913/2304298‑2304786
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                             
aaGCTGCCGATCGTGGTTTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGt                                                                                                                                                                                                                                                   <  2:388069/251‑1 (MQ=255)
aaGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCCGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCggg                                                                                                                                                                                                                                                   >  1:875474/1‑250 (MQ=255)
   gtgctgatcggtgggtATACGCCAGAAAAGGTAGATGTGAAGTCCCGCGGTGGGAAGGCGGAAGAGAAAATTGCCAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGGGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGttac                                                                                                                                                                                                                                                <  2:1007686/237‑4 (MQ=14)
   cTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGATAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAggaggatac                                                                                                                                                                                                                                                >  1:485334/1‑244 (MQ=21)
      ccGATCGTGGTGTATACGCCAGACAATTGAGATGTGAAGTACCGCGTCTGTAAGACGGAAGAAAAAATTGACACC‑GCGGTAGTTCGCTAAACTGCCGTGAAGGGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTAcgcc                                                                                                                                                                                                                                             <  2:1238213/251‑1 (MQ=17)
      ccGATCGTGGTGTATACGCCAGACAATGGAGAGGTGAAGTACCGCGTCGGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGATAAACTGCCGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTAcgcc                                                                                                                                                                                                                                             <  2:1158518/251‑1 (MQ=18)
        gATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTAcgctgc                                                                                                                                                                                                                                           >  2:1220214/1‑248 (MQ=255)
            gtggtgTATCCGCAAGACAATGTAGATGTGAAGTACCGCGGCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATcc                                                                                                                                                                                                                                       <  2:920700/251‑1 (MQ=17)
            gtggtgTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATcc                                                                                                                                                                                                                                       >  2:289909/1‑251 (MQ=255)
                   ataCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCgggttcgccgcagcccacatct                                                                                                                                                                                                                                >  2:638718/1‑231 (MQ=18)
                   ataCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATct                                                                                                                                                                                                                                <  1:1166439/251‑2 (MQ=255)
                     aCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACATcta                                                                                                                                                                                                                               >  1:759812/1‑248 (MQ=255)
                       gCCAGACAATTTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACACC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCTCCCCGTAGGTCGGAT‑AAGGCGGTTACGCCGCATCCGACagctaacg                                                                                                                                                                                                                            <  1:235692/251‑8 (MQ=17)
                          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)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                             
AAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAC‑GCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGAC‑AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCT‑AACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTG‑TGAAGTGCTCCA‑CCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCC  >  NC_000913/2304298‑2304786

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