Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,031,097 (GC)4→3 coding (1220‑1221/1419 nt) dcm ← DNA‑cytosine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,031,0910G.96.7% 256.7 / 1.4 61coding (1227/1419 nt)dcmDNA‑cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (23/35);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*NC_0009132,031,0920C.95.1% 257.1 / 4.1 61coding (1226/1419 nt)dcmDNA‑cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base . (23/35);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 7.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTC  >  NC_000913/2030843‑2031331
                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                               
tctcTATCGCCGTATCCCGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGGGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgc                                                                                                                                                                                                                                                 <  2:312272/250‑1 (MQ=255)
  tctATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAACCCAATTAAgcgcg                                                                                                                                                                                                                                              >  2:550205/1‑251 (MQ=255)
         ccGTACCACGCGGGGCAATCGCTGGCATATTTTTGCTGCGAGAGGCCTTACCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCTCCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTGCTCCCGGCGCTTCAAAACCCATTAAgcgcgcattccc                                                                                                                                                                                                                                       <  2:136411/251‑7 (MQ=255)
                   cgTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTaa                                                                                                                                                                                                                           <  1:620350/251‑1 (MQ=255)
                    gtgtaattcacgcacATATTTTTGCTGCGAGTGGCCTTATAGTGAACGTCGGACATGTTATGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTGTTGGATCAATCAGTTTTGCAGCCGCGACACAGACCGGCCAGACCACAGAGTTACCGAACTGGCGATAGGCCTGCAGGTGCGAAACCGGGATACGAAATTTCGCTTCTCCCGGCGCATCAAACCCCCTTAAGCGCCCATTCCCGAgcctgtaac                                                                                                                                                                                                                          >  2:529250/15‑243 (MQ=255)
                      ggCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCg                                                                                                                                                                                                                        >  2:353277/1‑251 (MQ=255)
                             gcgcgcATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGa                                                                                                                                                                                                                 <  1:265938/251‑1 (MQ=255)
                                 gcATATTTGTGCTGCGAGGGGCCTTATCGTGAACGCCGGCCATGTTGTGCCGATTGCTGCCGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCGCCCGGCGCTTTAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGAt                                                                                                                                                                                                             <  2:683763/251‑1 (MQ=255)
                                   atagttgtgctgcGAGTGGCCTTAGCGTGAACGTCGGCCATGTTGGGCCTCTTGCGGACGCAACGCCTCCGACTGTTTGATTTTTGGCGCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGt                                                                                                                                                                                                           <  1:603125/244‑1 (MQ=255)
                                       ttttTCTGCGAGTGGCCTTATCTTGAACTGGGGCCCTGTTGTGCCTCTTGCTGACGCAACGCCACCTCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACGGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGc                                                                                                                                                                                                        <  2:700535/250‑1 (MQ=255)
                                        tttGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCATGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGCATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTg                                                                                                                                                                                                      >  1:644198/1‑251 (MQ=255)
                                         ttGCTGGGAGTGGCCTTATCGTGAACGTCGGCCAGGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTg                                                                                                                                                                                                      <  1:458174/250‑1 (MQ=255)
                                           gctgcGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGAtt                                                                                                                                                                                                   <  2:205649/251‑1 (MQ=255)
                                            ctgcgcgtggtattATCCTGAAAGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTc                                                                                                                                                                                                  >  2:468688/13‑251 (MQ=255)
                                                gAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGCGGCGTTAACCGACGTGGACGATGTTGCTGCTTCAGCg                                                                                                                                                                                              >  2:530336/1‑251 (MQ=255)
                                                 aGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGCCCACCGAGTTCCCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGCTGTTGCTGCTTCAACgg                                                                                                                                                                                             >  2:681529/1‑251 (MQ=255)
                                                   tGGCCTTATCGTGAACGTCGGCCATTTTGGGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCGCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGa                                                                                                                                                                                            <  2:607505/250‑1 (MQ=255)
                                                    ggCCTTATCGGGAACTTCGTCCATGTTGTGCCTCTTGCTGACGCACCGCCGCCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGACCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGAt                                                                                                                                                                                           <  2:25337/250‑1 (MQ=255)
                                                        ttATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGCTGTTGCTGATTCAGCGGATCGTc                                                                                                                                                                                       >  2:217681/1‑250 (MQ=255)
                                                              tGAACGGCGGCCATGGTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATCTTTGGCTCAAGCAGTCTTGCCACCGAGGAAAAGAGGGGAAAGAAAAAAGAGGGACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATCCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCCAAGTCt                                                                                                                                                                                >  2:258982/1‑251 (MQ=255)
                                                                aaCGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCttt                                                                                                                                                                              <  1:204474/251‑1 (MQ=255)
                                                                     cGGCCATGTTGTCCCTATTTCTGACGCAACGACCCCGCCTGTTTGATTTTTGGCTCAAGCAGTTGTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTGCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTCTCTCa                                                                                                                                                                          <  2:407177/250‑1 (MQ=255)
                                                                                      ttGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCa                                                                                                                                                         <  1:353277/250‑1 (MQ=255)
                                                                                        gCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGcc                                                                                                                                                      >  1:211324/1‑251 (MQ=255)
                                                                                           gACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCATCcgcg                                                                                                                                                   <  1:226594/251‑1 (MQ=255)
                                                                                               cAACGCCACAGCATGTTTGAATTTTGGCTCACGCAGTATTGCCCCCGCGGCAAAGCCCGGCACGACCACCGAGTTACCGAACTGGCGCTAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCACGAGGCGTTAACCGACGTGGACGCTGTTGCTGATTCCGCGGATCGTCAAAGTCTTTCTCACACGTGGCCATATCCCAGCcg                                                                                                                                                     >  2:415012/1‑245 (MQ=255)
                                                                                                aaCGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGcgatcga                                                                                                                                              <  1:127017/251‑1 (MQ=255)
                                                                                                         gCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAAttt                                                                                                                                     <  1:681529/251‑1 (MQ=255)
                                                                                                          ccTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGCTTAAAATTTc                                                                                                                                    >  1:8642/1‑251 (MQ=255)
                                                                                                                  aTTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCAt                                                                                                                            <  1:475107/251‑1 (MQ=255)
                                                                                                                     tttGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATcg                                                                                                                          <  2:210170/250‑2 (MQ=255)
                                                                                                                        ggCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGt                                                                                                                      <  1:530336/251‑1 (MQ=255)
                                                                                                                                cAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAAcgc                                                                                                              >  1:15449/1‑251 (MQ=255)
                                                                                                                                      tGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGCTTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCCTCTTTGTAATAACGCGCAGAc                                                                                                        >  1:85724/1‑251 (MQ=255)
                                                                                                                                          aCCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGAAAGGCCTGAGTGTCCGAAACAGGAATACGGCAATCCACTTATCCCGGCAATTCACAACCCATTAAGCGAGCACA‑‑CCAGCGGCGTTAACCGACACGGCCGATGCTGCCGATTCCGCGGATCGTAAAAGTCTTTCTCACCCGTGGACATATCCCAGCCGCGATCGACTAAAATTTCCGCGCCCTCTTTGTAATAACGCGCAGAAAGcc                                                                                                    >  1:649407/1‑250 (MQ=255)
                                                                                                                                              cgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTTCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTAcg                                                                                                <  1:432632/251‑1 (MQ=255)
                                                                                                                                              cgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAAGAACGCGCAGACAGCGTAcg                                                                                                <  1:550205/251‑1 (MQ=255)
                                                                                                                                              cgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCAGTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTAcg                                                                                                <  2:345659/251‑1 (MQ=255)
                                                                                                                                                       ccGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCt                                                                                       <  1:484909/251‑1 (MQ=255)
                                                                                                                                                         ggccagACCACCGAGTTGCCGAACTGGCGAGTGGCCTGAGTGTCCGAATCCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCAGTAA‑‑GCGCGCATTCCCGTGGCGTGAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGGCAGCGTACGCGTGACGCttt                                                                                     <  2:509402/246‑1 (MQ=255)
                                                                                                                                                              gACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGAGGGACGCTTTGCggg                                                                                >  2:302462/1‑250 (MQ=255)
                                                                                                                                                                 caccGAGTTACCGAACTGGCGATAGGCCTGAGGGACCGACACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGGTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCGCCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTGGCGGATTg                                                                             <  2:167051/251‑1 (MQ=255)
                                                                                                                                                                    cGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTc                                                                          <  1:468688/251‑1 (MQ=255)
                                                                                                                                                                          aCCGAACGGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATa                                                                     <  2:577524/250‑1 (MQ=255)
                                                                                                                                                                                  ggCGATAGGCCTGAGTGTCCGAAACCGGAATACGGACTTTCGCTTCTCCCGGAGCTTCAAAACCCCTTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGCTGTTGCTGATTCAGAGGATCGTCAAAGTCTTTCTCACCAGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATACCGCGCAGACAGCGTACTCGTGACGCTTTGCGGATTGTTCGGATAAACCATTcc                                                            >  2:48594/1‑251 (MQ=255)
                                                                                                                                                                                  ggCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTcc                                                            >  1:544590/1‑251 (MQ=255)
                                                                                                                                                                                   gCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCCCCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCa                                                           <  2:280317/251‑1 (MQ=255)
                                                                                                                                                                                   gCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCa                                                           <  2:237845/251‑1 (MQ=255)
                                                                                                                                                                                        aGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAc                                                       <  1:242769/250‑1 (MQ=255)
                                                                                                                                                                                             tGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAACGTACGAGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGCAGc                                                 >  2:403311/1‑251 (MQ=255)
                                                                                                                                                                                                aGTGTCCGCAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGATGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCTATTAAAATTTCCGCGCCTTCTTTTTAATAACGCGCAGACAGCGTTCGCGTGACGCTTGGCGGATTGTTCGGATAAACCAGTTCATAACCGAAGCCg                                               <  2:460141/251‑1 (MQ=255)
                                                                                                                                                                                                 gtgtCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAAGCATTCCATAACCGAAGCCGt                                              >  1:383749/1‑251 (MQ=255)
                                                                                                                                                                                                     ccGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTAcc                                          <  1:6493/251‑1 (MQ=255)
                                                                                                                                                                                                       gAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGACAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACcgc                                        >  1:678585/1‑251 (MQ=255)
                                                                                                                                                                                                        aaaCCGGAATACGGAATTTCTCTTCCCCCGGCGCTTCAAACCCCATTAA‑‑GCGCGAATTCCCGAGTCGTTATCCGACGTGGACGAGGTTGCTTATGCAGCGGATCGTCAAAGTCTTTCTCGCCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTTTAATAACGCGCAGACAGCGTACGCGTGACGCTTGGCGTAGTGTTCGGAGAAGTCATTGCAGAACCGAATCCGTTACcgcg                                       <  2:608334/251‑1 (MQ=255)
                                                                                                                                                                                                           ccGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGCTTGTTCGGATAAACCATTCCATAACCGAAACCGTTACCgcgcgc                                    >  2:6098/1‑251 (MQ=255)
                                                                                                                                                                                                               aatacgtcaTTTCGCTTCTCACGGCGCTTCACAACTCATAAACCGCCCATTACCAAAGCGTTAACCACCGTGAGCGATTTTGCTGATTCAGCGGATCGTGAAAGTCTTTCTCACGCGCGTCCATATCCCAGCCGCGATCGATTCAAACTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTACGGATAAACCATTCAACAACCGAAGCCGTTACCGCGAGCCTg                                 >  2:92424/9‑250 (MQ=255)
                                                                                                                                                                                                                          tCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCg                      <  1:162052/250‑1 (MQ=255)
                                                                                                                                                                                                                              ttCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCatat                 >  2:541939/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tcCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCg               >  2:16165/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cccGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATc                <  1:404752/248‑1 (MQ=255)
                                                                                                                                                                                                                                  cccGGCGCTTCAAAACCCATTAA‑‑GCGCGCATGCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGAt             <  2:24521/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGAtcg           >  1:33640/1‑248 (MQ=255)
                                                                                                                                                                                                                                             aaaaCCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTc  >  1:66923/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                               
TCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTC  >  NC_000913/2030843‑2031331

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