Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,031,991 Δ1 bp coding (327/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,9910C.100.0% 161.9 / NA 39coding (327/1419 nt)dcmDNA cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (23/16);  total (23/16)

ACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTT  >  NC_000913/2031863‑2032122
                                                                                                                                |                                                                                                                                   
acGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAAt                                                                                                                          >  2:310066/1‑139 (MQ=255)
                tGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGa                                                                                                          <  1:753132/139‑1 (MQ=255)
                  aGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACg                                                                                                        >  2:86060/1‑139 (MQ=255)
                  aGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACg                                                                                                        >  2:307648/1‑139 (MQ=255)
                             gcgGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCg                                                                                                                               <  2:717108/105‑1 (MQ=255)
                             gcgGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCg                                                                                                                               >  1:717108/1‑105 (MQ=255)
                                     cTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCg                                                                                     >  2:392561/1‑139 (MQ=255)
                                       tCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGaa                                                                                   <  1:307648/139‑1 (MQ=255)
                                                     cGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCTCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGc                                                                     >  2:8009/1‑139 (MQ=255)
                                                      gTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCg                                                                    <  1:86060/139‑1 (MQ=255)
                                                                        tGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGgatg                                                  >  2:505806/1‑139 (MQ=255)
                                                                          gTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGgatgat                                                >  2:1037461/1‑139 (MQ=255)
                                                                             ggCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATa                                                                                 <  1:91387/103‑1 (MQ=255)
                                                                             ggCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATa                                                                                 >  2:91387/1‑103 (MQ=255)
                                                                             ggCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGAtgtg                                             >  2:658869/1‑139 (MQ=255)
                                                                             ggCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGAtgtg                                             >  1:267502/1‑139 (MQ=255)
                                                                               cTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGaa                                                                                   >  2:925524/1‑99 (MQ=255)
                                                                               cTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGaa                                                                                   <  1:925524/99‑1 (MQ=255)
                                                                                       gTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGttt                                   <  1:505806/139‑1 (MQ=255)
                                                                                         gCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTgg                                 >  1:718849/1‑139 (MQ=255)
                                                                                          cGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCg                                                                                        <  2:559087/83‑1 (MQ=255)
                                                                                          cGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCg                                                                                        >  1:559087/1‑83 (MQ=255)
                                                                                            tACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGaa                                                                  <  1:1023268/103‑1 (MQ=255)
                                                                                            tACCGCATGTTTGTTCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGaa                                                                  >  2:1023268/1‑103 (MQ=255)
                                                                                                       tgttCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGc                   >  1:928768/1‑139 (MQ=255)
                                                                                                         ttCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGa                                                                                                          >  2:1069579/1‑50 (MQ=38)
                                                                                                         ttCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGa                                                                                                          <  1:1069579/50‑1 (MQ=38)
                                                                                                          tCCATTCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGa                <  1:1037461/139‑1 (MQ=255)
                                                                                                              ttCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGa                                                                   <  1:917342/84‑1 (MQ=255)
                                                                                                              ttCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGa                                                                   >  2:917342/1‑84 (MQ=255)
                                                                                                              ttCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATggg            >  2:234527/1‑139 (MQ=255)
                                                                                                              ttCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATggg            >  2:1027209/1‑139 (MQ=255)
                                                                                                              ttCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATggg            >  1:447582/1‑139 (MQ=255)
                                                                                                               tCGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGc           >  2:506243/1‑139 (MQ=255)
                                                                                                                cGCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCg          <  1:43850/139‑1 (MQ=255)
                                                                                                                 gCTGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGa         >  1:311404/1‑139 (MQ=255)
                                                                                                                   tGGTAAACACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAAc       <  1:785702/139‑1 (MQ=255)
                                                                                                                        aaCACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCttt  <  1:69426/139‑1 (MQ=255)
                                                                                                                        aaCACGCATGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCttt  <  1:310066/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                   
ACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTT  >  NC_000913/2031863‑2032122

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