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.

GGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGC  >  NC_000913/2889200‑2889692
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
ggttcggttACCCCAGTAACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGGGGttt                                                                                                                                                                                                                                                        >  1:156798/1‑251 (MQ=255)
ggttcggtgGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTGCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGttt                                                                                                                                                                                                                                                        <  1:181437/251‑1 (MQ=255)
 gttcggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGtttt                                                                                                                                                                                                                                                       >  1:119650/1‑251 (MQ=255)
 gttcggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGtttt                                                                                                                                                                                                                                                       <  1:87599/251‑1 (MQ=255)
        tACCCCAGTCTCGCTGAGTTGTCACGACGGTTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAAttt                                                                                                                                                                                                                                                 <  1:287463/250‑1 (MQ=255)
          ccccAGTCACGCTGAGTTGTCACGACGGCTGCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTAcg                                                                                                                                                                                                                                              <  2:143555/251‑1 (MQ=255)
            ccAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTAcgcg                                                                                                                                                                                                                                            <  1:392923/251‑1 (MQ=255)
                 cACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCagg                                                                                                                                                                                                                                       <  1:316108/251‑1 (MQ=255)
                     cTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACCACCGTGGGTTTGAATTTTCGCGGCagggagg                                                                                                                                                                                                                                   >  1:162133/1‑251 (MQ=255)
                            gTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGTATTTACGCGGCAGGTAGGTCTGGcc                                                                                                                                                                                                                            >  1:151277/1‑251 (MQ=255)
                             tccctgtggttgcGGCCTCCGCCAGCGTATTCTCCTGCGGCATATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCCCCTCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCACGTAGGTCTGGcc                                                                                                                                                                                                                            <  2:147619/238‑1 (MQ=255)
                                       ttCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGGGGGGTTTGGTTTACGCGGCaggggggggttgccgcgggtcgggt                                                                                                                                                                                                                 >  2:76610/1‑229 (MQ=255)
                                          ggCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCCCTACCCGGGGTTTGGATTTACGCGGGaggtgggggggggcgagggtcgggtctt                                                                                                                                                                                                              >  2:238224/1‑227 (MQ=255)
                                           gCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCCGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGCTCACTACCCTGGGTTTGAATTTACGCGGCGGGGGGGTCTGGCCGAGGATCGGgtcttt                                                                                                                                                                                                             >  1:311790/1‑250 (MQ=255)
                                           gCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTAACTACCGTGGGTTTTAATTTTCGCGGCagggagggtgggccgggggtcgggtcttt                                                                                                                                                                                                             >  1:371993/1‑229 (MQ=255)
                                               ccGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTtcatca                                                                                                                                                                                                         <  1:74150/251‑1 (MQ=255)
                                                 gCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTTATCagt                                                                                                                                                                                                       >  1:249077/1‑251 (MQ=255)
                                                       gTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTGAATTTACGCGGGAGGTAGGGCTGGCCGAGGATCGGGTCTTTCTCagtaggggg                                                                                                                                                                                                 >  1:141244/1‑247 (MQ=255)
                                                           gCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACt                                                                                                                                                                                            <  1:561/251‑1 (MQ=255)
                                                             cccAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTGTGAATTTACGCGGGAGGTAGGTCTGGCCGAGGATCGGTTCTGCATCAGGAGTGGCGACttt                                                                                                                                                                                          <  2:24232/250‑1 (MQ=255)
                                                               cAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACttttt                                                                                                                                                                                        <  1:459809/251‑1 (MQ=255)
                                                               cAGCGGCAGATAGCCAAACTCACTCGCGGGGCGGGCGGAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATAGCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACttttt                                                                                                                                                                                        <  2:406994/251‑1 (MQ=255)
                                                                            ccAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGGGGTTTTTAATTTACGCGGCCGGTAGGGCTTGCCGGGGGTCGGTTCTTTCTCAGTAGTGGGGACttttttcttgttgcgcc                                                                                                                                                                            >  1:48466/1‑238 (MQ=255)
                                                                                   cactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTGGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGAtct                                                                                                                                                                     <  2:149181/249‑1 (MQ=255)
                                                                                        gcgGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGGGGCAGGTGGGTCTGGCCGGGGGACTGGTCTTCGTCACTACTGGCGACGTTTTTCTGGTCGAGCCAGATCTccacct                                                                                                                                                               >  2:79982/1‑246 (MQ=255)
                                                                                          ggTGCGGGCGTCGGTTTTCTTGTTGCCGTGTTCGCTGGAAAGCCAACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCAATCCCCGGGGTTTTAATTTTCGCCGCCGGCCGGTCTGGCCGAGGGTCGGTTATTCCTCCACAGTGGCGGTTTTTTCCTgggcgggcccgattctcccctcc                                                                                                                                                              >  2:313304/1‑229 (MQ=255)
                                                                                           gTGCGGGCGTAGGTTTTCTTGGTGCCGTGTTCGATGGAAAGCCCAACGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGCCTTTTTCCTGGTCGAGCCAGATCTCCGCATAcgc                                                                                                                                                           <  1:146700/251‑1 (MQ=255)
                                                                                             gCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATAcgcgc                                                                                                                                                         <  2:371993/251‑1 (MQ=255)
                                                                                                 gCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGt                                                                                                                                                     <  2:180690/251‑1 (MQ=255)
                                                                                                                    cGTGTACGATTGAAAGCCCCCCGCCCTCCAACAGGGTAAAGCCCACCAGCTTGCCGTTTTCTGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTGCCTGGTCGATCCAGATCTCCGCATACGCGCGGGTACGAGGCAAAAGATGCt                                                                                                                                   <  2:325724/249‑1 (MQ=255)
                                                                                                                     gtgtTCGATGGGAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCCATACCGCGGTTTTTAATTTACGCGGCCGGTAGGGCTGGCCGAGGATCGGGTCTTCCTCAGTAGTGGCGACTTTTTTCTGGGCGGGCCAGATCTCCCCATACGCCCGGGtgccaggccacacatgtcccg                                                                                                                                 >  2:174537/1‑231 (MQ=255)
                                                                                                                         tCGATGGAAAGCCCTCCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCTATCGCCACGAAGTTCATGTCGGTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTTGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGaaa                                                                                                                             <  1:502349/250‑1 (MQ=255)
                                                                                                                          cGATGGAAAGCCCACCGCCCACCAAAAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTTAATTTACGCGGGCGGTAGGTCTGGCCGAGGGGCGGTTCTTCATCAGTAGTGGCGACTTTTTTCTGGGCGAGCCAGATCTCCGCATACGCGCGGGGACGCGGCCACAGATGCTCAGAAAtc                                                                                                                           >  2:487268/1‑251 (MQ=255)
                                                                                                                           gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAtct                                                                                                                          <  2:151277/251‑1 (MQ=255)
                                                                                                                              ggAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAtcttc                                                                                                                        <  1:244039/250‑1 (MQ=255)
                                                                                                                               gAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTTCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATcttctt                                                                                                                      <  1:339207/251‑1 (MQ=255)
                                                                                                                                 aaGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCg                                                                                                                    <  1:161197/251‑1 (MQ=255)
                                                                                                                                   gcccaccgcccaccAACAGGTTCAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTGCAATTACGCGGCAAGGAGGGCTTGGCGAGGGTAGGGTCTTCATCAGTAGTGGCGACCTTTTCCTTGGCGGGCCAGATCTCCGCCTACGCGCGGGGACGAAGCGAAAGATGGCCAGAAAACTTCTTTGcc                                                                                                                  >  1:396845/1‑251 (MQ=255)
                                                                                                                                      cgccgcccaccAACAGGTTAAAGCCCTCCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGGGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGTCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCAc                                                                                                               <  1:101281/249‑1 (MQ=255)
                                                                                                                                           cccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTTAATTTACGCGGCAGGTAGGTCTGGCCGGGGATCGGGTCTTTCTCAGTAGGGGGGGATTTTTTCCGGGCGAGCCAAGACTCCCCCTAAGCGCGGGGACGAGGCACAAGATGCCCAGAAATCTTCTTCTCCCCCTCGt                                                                                                           >  2:365137/1‑250 (MQ=255)
                                                                                                                                                 aaCAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACCAACGGGGGTTTGGATTTTCGCGGCACGTAGGTCTGGCCGGGGCTCGGTTTTTTATCCGTAGTGGCGCCTTTTTCCCTTTCTCGCCCGATTTCCCCCTACCGGCGGGGACCAGGCAACAGAGGCGCATAAATCTTTTTCCCCCCCCCGTACGCTTc                                                                                                   >  2:457275/1‑251 (MQ=255)
                                                                                                                                                    aGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTGGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTTCTGGTCGAGCCAGATCTCCCCATACGCGCGGGGACGAGGCAACAGATGCTCCGAAATCTTCTTCCCCCCCTCGTACGCTCCcgc                                                                                                >  1:114953/1‑251 (MQ=255)
                                                                                                                                                               caccaGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTCCGCTTCCGCGTGCAGCTGca                                                                                     <  2:114953/251‑2 (MQ=255)
                                                                                                                                                                   aGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTc                                                                                 <  1:76610/251‑1 (MQ=255)
                                                                                                                                                                                  gCGATCGCCACGAAGTTCATGTCGTGGGCGTGCAGATCGATATCGTGCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCATGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCGGCGACTCGTAAGGGTTCGAGGt                                                                  <  2:341005/251‑1 (MQ=255)
                                                                                                                                                                                     aTCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGGGGGTTTGAATTTTCGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCCTCAGTAGTGGCGCCTTTTTCCTGGTCGAGCCAGATTTCCGCATACGCGCGGGTACGGGGGAACCAATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGGCTCGTACGGGTTCGAggggc                                                                >  1:238925/1‑247 (MQ=255)
                                                                                                                                                                                        gCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATGACTACCGTGGGTTTGGATTTACGGGGGCGGTGGGTGTGGACGCGGATCGGTTCTTTATTCGGAGAGGCGACTTTTTCCTGGACGAGCCAGATCTCCGCGTACGCACGGGTACGAGGGCAAAAATGCACAGAAATCTTCTTTGCCCCATCGTGAGATTCCGCGGGCAGCCGCGCATCGTCAGGGGTGGAGGTGCag                                                              >  1:84906/1‑249 (MQ=255)
                                                                                                                                                                                              aaGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGCAGGGAGGTCTGGCCGAGGAGCGGTTCTTCCTCAGTAGTGGCGGCTTTTTTCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAAACGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGAATCGTAAGGGGTCgggggggcgagtacggt                                                      >  1:484083/1‑235 (MQ=255)
                                                                                                                                                                                               aGTTCATGTCGTTGGCGTGCAGATCGATTTCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTGCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTa                                                     <  2:81902/251‑1 (MQ=255)
                                                                                                                                                                                                    aTGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTTCCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGtt                                                <  1:194916/251‑1 (MQ=255)
                                                                                                                                                                                                    aTGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCGTCATCAGTAGTGGCGACTTTTGCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGtt                                                <  1:492882/251‑1 (MQ=255)
                                                                                                                                                                                                      gTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTGAATTTACGCGGCAGGTAGGGCTCGCCGAGGATCGGGTCTTCATCAGTAGAGGCGACTTTTTCCTGGTCGAGCCAAATCTCCGCCTACGCGAGGGTACGCGGCAACAAATGGTCAGAAATCTTCTTCGCCCCATCGTACGCTCTCCCGTGCAGCTGCGCCACGTAAGGGTGCGAGGTGCAGAGTACGGTACGGTTcc                                              >  1:236730/1‑250 (MQ=255)
                                                                                                                                                                                                          ttGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTAAGCGGGAGGAGGGTCTGGCCGAGGATCGGGTTTTTATCAGTAGTGGCGCCTTTTTCCTGGTCGAGCCAAATCTCCGCTAAAGCGCGGGGACGAGGCAACAGATTCTCAGAAAACTTCTTTGCCCACTTGTCCGCTTCCGCGAGCATCTGCGACTTGTAAGGGGTCTAGGGGCAAAGTATGTTatgggtcctgtg                                          >  1:13435/1‑240 (MQ=255)
                                                                                                                                                                                                                ggCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGc                                    <  1:244423/250‑1 (MQ=255)
                                                                                                                                                                                                                 gCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGGAGTGGGGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGGACGAGGGAACAGATGCTCCGAAAACTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCt                                   >  1:394291/1‑251 (MQ=255)
                                                                                                                                                                                                                   aGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTGAAATTACGCGGCAGGGAGGGCTGGCCGGGGATCGGTTCTTCATCAGTAGGGGCGACTTTTTTCTGGTCGAGCCAGATCTCCGCATACGCGCGGGGACGGGGCAACAGATGCTCAGAAATCTTTTTCGCCCCCTCCTACGCTTCCGCGTGCAGCTGCGGCTCGTAAGGGGTTGGGGGGCCGAGTACGTTACGGTTCATGTCGTTAGCTGt                                 >  1:340004/1‑251 (MQ=255)
                                                                                                                                                                                                                      tcgttatCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGACGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGAGCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAGCTTCTTCGCCCGCGCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGATGCAGAGTACGTTACGGTCTATGTCGTCAGCTGttg                               <  2:224920/246‑3 (MQ=255)
                                                                                                                                                                                                                                     ggCGGGATCACTACCGTGGTTTTGAATTTACGCGGCATGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCGTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAgac                <  2:340004/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gggATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAtc                                                                                                                           >  2:145550/1‑141 (MQ=255)
                                                                                                                                                                                                                                           aTCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGAATTTTTCATGGTCGAGCCAGATCTCCGCATACGCGCGGGAACAAGGCAAAAACTGCTCAGAAATCTTCTTCGCCAACTCGTACGCTTCCGCGTGCGGCTGCGCCTCGTAAGGGTTCGAGTAGCAGAATAGGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGaccgaccg          >  1:502646/1‑250 (MQ=255)
                                                                                                                                                                                                                                            tCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGaccgaccg          >  1:389824/1‑249 (MQ=255)
                                                                                                                                                                                                                                              aCTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTg      <  2:156798/251‑1 (MQ=255)
                                                                                                                                                                                                                                                tACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgc     >  2:186910/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   cGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGGgcagc  >  2:105693/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
GGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGC  >  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: