Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,4480TG28.3% 97.6 / 12.9 60T218P (ACC→CCC) cysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (14/29);  new base G (15/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGC  >  NC_000913/2889200‑2889692
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
ggttcggttACCCCAGTAACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACC‑GGGGt‑tt                                                                                                                                                                                                                                                        >  1:156798/1‑251 (MQ=255)
ggttcggtgGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTGCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGt‑tt                                                                                                                                                                                                                                                        <  1:181437/251‑1 (MQ=255)
 gttcggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGt‑ttt                                                                                                                                                                                                                                                       >  1:119650/1‑251 (MQ=255)
 gttcggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGt‑ttt                                                                                                                                                                                                                                                       <  1:87599/251‑1 (MQ=255)
        tACCCCAGTCTCGCTGAGTTGTCACGACGGTTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAAttt                                                                                                                                                                                                                                                 <  1:287463/250‑1 (MQ=255)
          ccccAGTCACGCTGAGTTGTCACGACGGCTGCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTAcg                                                                                                                                                                                                                                              <  2:143555/251‑1 (MQ=255)
            ccAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTAcgcg                                                                                                                                                                                                                                            <  1:392923/251‑1 (MQ=255)
                 cACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCagg                                                                                                                                                                                                                                       <  1:316108/251‑1 (MQ=255)
                     cTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACCACC‑GTGGG‑TTTGAATTTTCGCGGCagggagg                                                                                                                                                                                                                                   >  1:162133/1‑251 (MQ=255)
                            gTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGTATTTACGCGGCAGGTAGGTCTGGcc                                                                                                                                                                                                                            >  1:151277/1‑251 (MQ=255)
                             tccctgtggttgcGGCCTCCGCCAGCGTATTCTCCTGCGGCATATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCCCCTCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCACGTAGGTCTGGcc                                                                                                                                                                                                                            <  2:147619/238‑1 (MQ=255)
                                       ttCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACC‑GGGGGGTTTGG‑TTTACGCGGCaggggggggttgccgcgggtcgggt                                                                                                                                                                                                                 >  2:76610/1‑229 (MQ=255)
                                          ggCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCCCTACC‑CGGGG‑TTTGGATTTACGCGGGaggtgggggggggcgagggtcgggtctt                                                                                                                                                                                                              >  2:238224/1‑227 (MQ=255)
                                           gCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCCGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGCTCACTACC‑CTGGG‑TTTGAATTTACGCGGCGGGGGGGTCTGGCCGAGGATCGGgtcttt                                                                                                                                                                                                             >  1:311790/1‑250 (MQ=255)
                                           gCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTAACTACC‑GTGGG‑TTTTAATTTTCGCGGCagggagggtgggccgggggtcgggtcttt                                                                                                                                                                                                             >  1:371993/1‑229 (MQ=255)
                                               ccGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTtcatca                                                                                                                                                                                                         <  1:74150/251‑1 (MQ=255)
                                                 gCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTTATCagt                                                                                                                                                                                                       >  1:249077/1‑251 (MQ=255)
                                                       gTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACC‑GTGGG‑TTTGAATTTACGCGGGAGGTAGGGCTGGCCGAGGATCGGGTCTTTCTCagtaggggg                                                                                                                                                                                                 >  1:141244/1‑247 (MQ=255)
                                                           gCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACt                                                                                                                                                                                            <  1:561/251‑1 (MQ=255)
                                                             cccAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGG‑TGTGAATTTACGCGGGAGGTAGGTCTGGCCGAGGATCGGTTCTGCATCAGGAGTGGCG‑ACttt                                                                                                                                                                                          <  2:24232/250‑1 (MQ=255)
                                                               cAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACttttt                                                                                                                                                                                        <  1:459809/251‑1 (MQ=255)
                                                               cAGCGGCAGATAGCCAAACTCACTCGCGGGGCGGGCGGAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATAGCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACttttt                                                                                                                                                                                        <  2:406994/251‑1 (MQ=255)
                                                                            ccAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACC‑GGGGT‑TTTTAATTTACGCGGCCGGTAGGGCTTGCCGGGGGTCGGTTCTTTCTCAGTAGTGGGG‑ACttttttcttgttgcgcc                                                                                                                                                                            >  1:48466/1‑238 (MQ=255)
                                                                                   cactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTGGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GAtct                                                                                                                                                                     <  2:149181/249‑1 (MQ=255)
                                                                                        gcgGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GGGGT‑TTTGAATTTACGGGGCAGGTGGGTCTGGCCGGGGGACTGGTCTTCGTCACTACTGGCG‑ACGTTTTTCTGGTCGAGCCA‑GATCTccacct                                                                                                                                                               >  2:79982/1‑246 (MQ=255)
                                                                                          ggTGCGGGCGTCGGTTTTCTTGTTGCCGTGTTCGCTGGAAAGCCAACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCAATCCCCGGGGT‑TTT‑AATTTTCGCCGCCGGCCGGTCTGGCCGAGGGTCGGTTATTCCTCCACAGTGGCG‑GTTTTTTCCTgggcgggcccgattctcccctcc                                                                                                                                                              >  2:313304/1‑229 (MQ=255)
                                                                                           gTGCGGGCGTAGGTTTTCTTGGTGCCGTGTTCGATGGAAAGCCCAACGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GGGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑CCTTTTTCCTGGTCGAGCCA‑GATCTCCGCATAcgc                                                                                                                                                           <  1:146700/251‑1 (MQ=255)
                                                                                             gCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATAcgcgc                                                                                                                                                         <  2:371993/251‑1 (MQ=255)
                                                                                                 gCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGt                                                                                                                                                     <  2:180690/251‑1 (MQ=255)
                                                                                                                    cGTGTACGATTGAAAGCCCCCCGCCCTCCAACAGGGTAAAGCCCACCAGCTTGCCGTTTTCTGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTGCCTGGTCGATCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAAAAGATGCt                                                                                                                                   <  2:325724/249‑1 (MQ=255)
                                                                                                                     gtgtTCGATGGGAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCCATACC‑GCGGT‑TTTTAATTTACGCGGCCGGTAGGGCTGGCCGAGGATCGGGTCTTCCTCAGTAGTGGCG‑ACTTTTTTCTGGGCGGGCCA‑GATCTCCCCATACGCCCGGGtgccaggccacacatgtcccg                                                                                                                                 >  2:174537/1‑231 (MQ=255)
                                                                                                                         tCGATGGAAAGCCCTCCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCTATCGCCACGAAGTTCATGTCGGTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTTGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGaaa                                                                                                                             <  1:502349/250‑1 (MQ=255)
                                                                                                                          cGATGGAAAGCCCACCGCCCACCAAAAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGG‑TTTTAATTTACGCGGGCGGTAGGTCTGGCCGAGGGGCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTTCTGGGCGAGCCA‑GATCTCCGCATACGCGCGGGGACGCGG‑CCACAGATGCTCAGAAAtc                                                                                                                           >  2:487268/1‑251 (MQ=255)
                                                                                                                           gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAAtct                                                                                                                          <  2:151277/251‑1 (MQ=255)
                                                                                                                              ggAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAAtcttc                                                                                                                        <  1:244039/250‑1 (MQ=255)
                                                                                                                               gAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTTCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATcttctt                                                                                                                      <  1:339207/251‑1 (MQ=255)
                                                                                                                                 aaGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGG‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCg                                                                                                                    <  1:161197/251‑1 (MQ=255)
                                                                                                                                   gcccaccgcccaccAACAGGTTCAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGG‑TTTGCAATTACGCGGCAAGGAGGGCTTGGCGAGGGTAGGGTCTTCATCAGTAGTGGCG‑ACCTTTTCCTTGGCGGGCCA‑GATCTCCGCCTACGCGCGGGGACGAAG‑CGAAAGATGGCCAGAAAACTTCTTTGcc                                                                                                                  >  1:396845/1‑251 (MQ=255)
                                                                                                                                      cgccgcccaccAACAGGTTAAAGCCCTCCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGGGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGTCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCAc                                                                                                               <  1:101281/249‑1 (MQ=255)
                                                                                                                                           cccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACC‑GTGGG‑TTTTAATTTACGCGGCAGGTAGGTCTGGCCGGGGATCGGGTCTTTCTCAGTAGGGGGGGATTTTTTCC‑GGGCGAGCCAAGA‑CTCCCCCTAAGCGCGGGGACGAGG‑CACAAGATGCCCAGAAATCTTCTTCTCCCCCTCGt                                                                                                           >  2:365137/1‑250 (MQ=255)
                                                                                                                                                 aaCAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACCAAC‑GGGGG‑TTTGGATTTTCGCGGCACGTAGGTCTGGCCGGGGCTCGGTTTTTTATCCGTAGTGGCG‑CCTTTTTCCCTTTCTCGCCC‑GATTTCCCCCTACCGGCGGGGACCAGG‑CAACAGAGGCGCATAAATCTTTTTCCCCCCCCCGTACGCT‑Tc                                                                                                   >  2:457275/1‑251 (MQ=255)
                                                                                                                                                    aGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTGGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTTCTGGTCGAGCCA‑GATCTCCCCATACGCGCGGGGACGAGG‑CAACAGATGCTCCGAAATCTTCTTCCCCCCCTCGTACGCT‑CCcgc                                                                                                >  1:114953/1‑251 (MQ=255)
                                                                                                                                                               caccaGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTCCGCT‑TCCGCGTGCAGCTGca                                                                                     <  2:114953/251‑2 (MQ=255)
                                                                                                                                                                   aGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTc                                                                                 <  1:76610/251‑1 (MQ=255)
                                                                                                                                                                                  gCGATCGCCACGAAGTTCATGTCGTGGGCGTGCAGATCGATATCGTGCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCATGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCGGCGACTCGTAAGGGTTCGAGGt                                                                  <  2:341005/251‑1 (MQ=255)
                                                                                                                                                                                     aTCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACC‑GGGGG‑TTTGAATTTTCGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCCTCAGTAGTGGCG‑CCTTTTTCCTGGTCGAGCCA‑GATTTCCGCATACGCGCGGGTACGGGG‑GAACCAATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGGCTCGTACGGGTTCGAggggc                                                                >  1:238925/1‑247 (MQ=255)
                                                                                                                                                                                        gCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATGACTACC‑GTGGG‑TTTGGATTTACGGGGGCGGTGGGTGTGGACGCGGATCGGTTCTTTATTCGGAGAGGCG‑ACTTTTTCCTGGACGAGCCA‑GATCTCCGCGTACGCACGGGTACGAGGGCAA‑AAATGCACAGAAATCTTCTTTGCCCCATCGTGAGAT‑TCCGCGGGCAGCCGCGCATCGTCAGGGGTGGAGGTGCag                                                              >  1:84906/1‑249 (MQ=255)
                                                                                                                                                                                              aaGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GGGGT‑TTTGAATTTACGCGGCAGGGAGGTCTGGCCGAGGAGCGGTTCTTCCTCAGTAGTGGCG‑GCTTTTTTCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAAACGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGAATCGTAAGGGGTCgggggggcgagtacggt                                                      >  1:484083/1‑235 (MQ=255)
                                                                                                                                                                                               aGTTCATGTCGTTGGCGTGCAGATCGATTTCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑GCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTa                                                     <  2:81902/251‑1 (MQ=255)
                                                                                                                                                                                                    aTGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTTCC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGtt                                                <  1:194916/251‑1 (MQ=255)
                                                                                                                                                                                                    aTGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCGTCATCAGTAGTGGCG‑ACTTTTGCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGtt                                                <  1:492882/251‑1 (MQ=255)
                                                                                                                                                                                                      gTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGG‑TTTGAATTTACGCGGCAGGTAGGGCTCGCCGAGGATCGGGTCTTCATCAGTAGAGGCG‑ACTTTTTCCTGGTCGAGCCA‑AATCTCCGCCTACGCGAGGGTACGCGG‑CAACAAATGGTCAGAAATCTTCTTCGCCCCATCGTACGCTCTCC‑CGTGCAGCTGCGCCACGTAAGGGTGCGAGGTGCAGAGTACGGTACGGTTcc                                              >  1:236730/1‑250 (MQ=255)
                                                                                                                                                                                                          ttGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTAAGCGGGAGGAGGGTCTGGCCGAGGATCGGGTTTTTATCAGTAGTGGCG‑CCTTTTTCCTGGTCGAGCCA‑AATCTCCGCTAAAGCGCGGGGACGAGG‑CAACAGATTCTCAGAAAACTTCTTTGCCCACTTGTCCGCT‑TCCGCGAGCATCTGCGACTTGTAAGGGGTCTAGGGGCAAAGTATGTTatgggtcctgtg                                          >  1:13435/1‑240 (MQ=255)
                                                                                                                                                                                                                ggCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGc                                    <  1:244423/250‑1 (MQ=255)
                                                                                                                                                                                                                 gCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGG‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGGAGTGGGG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGGACGAGG‑GAACAGATGCTCCGAAAACTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCt                                   >  1:394291/1‑251 (MQ=255)
                                                                                                                                                                                                                   aGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGG‑TTTGAAATTACGCGGCAGGGAGGGCTGGCCGGGGATCGGTTCTTCATCAGTAGGGGCG‑ACTTTTTTCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGGACGGGG‑CAACAGATGCTCAGAAATCTTTTTCGCCCCCTCCTACGCT‑TCCGCGTGCAGCTGCGGCTCGTAAGGGGTTGGGGGGCCGAGTACGTTACGGTTCATGTCGTTAGCTGt                                 >  1:340004/1‑251 (MQ=255)
                                                                                                                                                                                                                      tcgttatCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGACGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GAGCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAAGCTTCTTCGCCCGCGCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGATGCAGAGTACGTTACGGTCTATGTCGTCAGCTGttg                               <  2:224920/246‑3 (MQ=255)
                                                                                                                                                                                                                                     ggCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCATGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCG‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAgac                <  2:340004/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gggATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAAtc                                                                                                                           >  2:145550/1‑141 (MQ=255)
                                                                                                                                                                                                                                           aTCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑AATTTTTCATGGTCGAGCCA‑GATCTCCGCATACGCGCGGGAACAAGG‑CAAAAACTGCTCAGAAATCTTCTTCGCCAACTCGTACGCT‑TCCGCGTGCGGCTGCGCCTCGTAAGGGTTCGAGTAGCAGAATAGGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGaccgaccg          >  1:502646/1‑250 (MQ=255)
                                                                                                                                                                                                                                            tCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGaccgaccg          >  1:389824/1‑249 (MQ=255)
                                                                                                                                                                                                                                              aCTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTg      <  2:156798/251‑1 (MQ=255)
                                                                                                                                                                                                                                                tACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgc     >  2:186910/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   c‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGGgcagc  >  2:105693/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
GGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACC‑GTGGT‑TTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCA‑GATCTCCGCATACGCGCGGGTACGAGG‑CAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCT‑TCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGC  >  NC_000913/2889200‑2889692

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