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

TCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTC  >  NC_000913/2030843‑2031331
                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                               
tctcTATCGCCGTATCCCGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGGGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgc                                                                                                                                                                                                                                                 <  2:312272/250‑1 (MQ=255)
  tctATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAACCCAATTAAgcgcg                                                                                                                                                                                                                                              >  2:550205/1‑251 (MQ=255)
         ccGTACCACGCGGGGCAATCGCTGGCATATTTTTGCTGCGAGAGGCCTTACCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCTCCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTGCTCCCGGCGCTTCAAAACCCATTAAgcgcgcattccc                                                                                                                                                                                                                                       <  2:136411/251‑7 (MQ=255)
                   cgTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTaa                                                                                                                                                                                                                           <  1:620350/251‑1 (MQ=255)
                    gtgtaattcacgcacATATTTTTGCTGCGAGTGGCCTTATAGTGAACGTCGGACATGTTATGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTGTTGGATCAATCAGTTTTGCAGCCGCGACACAGACCGGCCAGACCACAGAGTTACCGAACTGGCGATAGGCCTGCAG‑GTGCGAAACCGGGATACGAAATTTCGCTTCTCCCGGCGCATCAAACCCCCTTAAGCGC‑C‑CATTCCCGAgcctgtaac                                                                                                                                                                                                                          >  2:529250/15‑243 (MQ=255)
                      ggCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCg                                                                                                                                                                                                                        >  2:353277/1‑251 (MQ=255)
                             gcgcgcATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGa                                                                                                                                                                                                                 <  1:265938/251‑1 (MQ=255)
                                 gcATATTTGTGCTGCGAGGGGCCTTATCGTGAACGCCGGCCATGTTGTGCCGATTGCTGCCGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCGCCCGGCGCTTTAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGAt                                                                                                                                                                                                             <  2:683763/251‑1 (MQ=255)
                                   atagttgtgctgcGAGTGGCCTTAGCGTGAACGTCGGCCATGTTGGGCCTCTTGCGGACGCAACGCCTCCGACTGTTTGATTTTTGGCGCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGt                                                                                                                                                                                                           <  1:603125/244‑1 (MQ=255)
                                       ttttTCTGCGAGTGGCCTTATCTTGAACTGGGGCCCTGTTGTGCCTCTTGCTGACGCAACGCCACCTCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACGGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGc                                                                                                                                                                                                        <  2:700535/250‑1 (MQ=255)
                                        tttGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCATGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGCATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTg                                                                                                                                                                                                      >  1:644198/1‑251 (MQ=255)
                                         ttGCTGGGAGTGGCCTTATCGTGAACGTCGGCCAGGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTg                                                                                                                                                                                                      <  1:458174/250‑1 (MQ=255)
                                           gctgcGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGAtt                                                                                                                                                                                                   <  2:205649/251‑1 (MQ=255)
                                            ctgcgcgtggtattATCCTGAAAGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTc                                                                                                                                                                                                  >  2:468688/13‑251 (MQ=255)
                                                gAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGCGGCGTTAACCGACGTGGACGATGTTGCTGCTTCAGCg                                                                                                                                                                                              >  2:530336/1‑251 (MQ=255)
                                                 aGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGCCCACCGAGTTCCCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGCTGTTGCTGCTTCAACgg                                                                                                                                                                                             >  2:681529/1‑251 (MQ=255)
                                                   tGGCCTTATCGTGAACGTCGGCCATTTTGGGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCGCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGa                                                                                                                                                                                            <  2:607505/250‑1 (MQ=255)
                                                    ggCCTTATCGGGAACTTCGTCCATGTTGTGCCTCTTGCTGACGCACCGCCGCCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGACCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGAt                                                                                                                                                                                           <  2:25337/250‑1 (MQ=255)
                                                        ttATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGCTGTTGCTGATTCAGCGGATCGTc                                                                                                                                                                                       >  2:217681/1‑250 (MQ=255)
                                                              tGAACGGCGGCCATGGTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATCTTTGGCTCAAGCAGTCTTGCCACCGAGGAAAAGAGGGGAAAGAAAAAAGAGGGACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATCCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCCAAGTCt                                                                                                                                                                                >  2:258982/1‑251 (MQ=255)
                                                                aaCGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCttt                                                                                                                                                                              <  1:204474/251‑1 (MQ=255)
                                                                     cGGCCATGTTGTCCCTATTTCTGACGCAACGACCCCGCCTGTTTGATTTTTGGCTCAAGCAGTTGTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTGCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTCTCTCa                                                                                                                                                                          <  2:407177/250‑1 (MQ=255)
                                                                                      ttGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCa                                                                                                                                                         <  1:353277/250‑1 (MQ=255)
                                                                                        gCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGcc                                                                                                                                                      >  1:211324/1‑251 (MQ=255)
                                                                                           gACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCATCcgcg                                                                                                                                                   <  1:226594/251‑1 (MQ=255)
                                                                                               cAACGCCACAGCATGTTTGAATTTTGGCTCACGCAGTATTGCCCCCGCGGCAAAGCCCGGCACGACCACCGAGTTACCGAACTGGCGCTAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCACGAGGCGTTAACCGACGTGGACGCTGTTGCTGATTCCGCGGATCGTCAAAGTCTTTCTCACACGTGGCCATATCCCAGCcg                                                                                                                                                     >  2:415012/1‑245 (MQ=255)
                                                                                                aaCGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGcgatcga                                                                                                                                              <  1:127017/251‑1 (MQ=255)
                                                                                                         gCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAAttt                                                                                                                                     <  1:681529/251‑1 (MQ=255)
                                                                                                          ccTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGCTTAAAATTTc                                                                                                                                    >  1:8642/1‑251 (MQ=255)
                                                                                                                  aTTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCAt                                                                                                                            <  1:475107/251‑1 (MQ=255)
                                                                                                                     tttGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATcg                                                                                                                          <  2:210170/250‑2 (MQ=255)
                                                                                                                        ggCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGt                                                                                                                      <  1:530336/251‑1 (MQ=255)
                                                                                                                                cAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAAcgc                                                                                                              >  1:15449/1‑251 (MQ=255)
                                                                                                                                      tGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGCTTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCCTCTTTGTAATAACGCGCAGAc                                                                                                        >  1:85724/1‑251 (MQ=255)
                                                                                                                                          aCCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGAAAGGCCTG‑AGTGTCCGAAACAGGAATACGGCAATCCACTTATCCCGGCAATTCACAACCCATTAAGCGAGCACA‑‑CCAGCGGCGTTAACCGACACGGCCGATGCTGCCGATTCCGCGGATCGTAAAAGTCTTTCTCACCCGTGGACATATCCCAGCCGCGATCGACTAAAATTTCCGCGCCCTCTTTGTAATAACGCGCAGAAAGcc                                                                                                    >  1:649407/1‑250 (MQ=255)
                                                                                                                                              cgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTTCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTAcg                                                                                                <  1:432632/251‑1 (MQ=255)
                                                                                                                                              cgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAAGAACGCGCAGACAGCGTAcg                                                                                                <  1:550205/251‑1 (MQ=255)
                                                                                                                                              cgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCAGTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTAcg                                                                                                <  2:345659/251‑1 (MQ=255)
                                                                                                                                                       ccGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCt                                                                                       <  1:484909/251‑1 (MQ=255)
                                                                                                                                                         ggccagACCACCGAGTTGCCGAACTGGCGAGTGGCCTG‑AGTGTCCGAATCCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCAGTAA‑‑GCGCGCATTCCCGTGGCGTGAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGGCAGCGTACGCGTGACGCttt                                                                                     <  2:509402/246‑1 (MQ=255)
                                                                                                                                                              gACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGAGGGACGCTTTGCggg                                                                                >  2:302462/1‑250 (MQ=255)
                                                                                                                                                                 caccGAGTTACCGAACTGGCGATAGGCCTG‑AGGGACCGACACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGGTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCGCCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTGGCGGATTg                                                                             <  2:167051/251‑1 (MQ=255)
                                                                                                                                                                    cGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTc                                                                          <  1:468688/251‑1 (MQ=255)
                                                                                                                                                                          aCCGAACGGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATa                                                                     <  2:577524/250‑1 (MQ=255)
                                                                                                                                                                                  ggCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGACTTTCGCTTCTCCCGGAGCTTCAAAACCCCTTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGCTGTTGCTGATTCAGAGGATCGTCAAAGTCTTTCTCACCAGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATACCGCGCAGACAGCGTACTCGTGACGCTTTGCGGATTGTTCGGATAAACCATTcc                                                            >  2:48594/1‑251 (MQ=255)
                                                                                                                                                                                  ggCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTcc                                                            >  1:544590/1‑251 (MQ=255)
                                                                                                                                                                                   gCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCCCCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCa                                                           <  2:280317/251‑1 (MQ=255)
                                                                                                                                                                                   gCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCa                                                           <  2:237845/251‑1 (MQ=255)
                                                                                                                                                                                        aGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑GCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAc                                                       <  1:242769/250‑1 (MQ=255)
                                                                                                                                                                                             tG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAA‑‑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)
                                                                                                                                                                                                               aatacgtcaTTTCGCTTCTCACGGCGCTTCACAACTCAT‑AAAC‑CGCCCATTACCAAAGCGTTAACCACCGTGAGCGATTTTGCTGATTCAGCGGATCGTGAAAGTCTTTCTCACGCGCGTCCATATCCCAGCCGCGATCGATTCAAACTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTACGGATAAACCATTCAACAACCGAAGCCGTTACCGCGAGCCTg                                 >  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)
                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                               
TCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTG‑AGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTC  >  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: ‑