Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,031,309	(T)6→7	coding (1009/1419 nt)	dcm ←	DNA‑cytosine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,031,303	1	.	T	100.0%	155.7 / NA	44	H339N (CAT→AAT)	dcm	DNA‑cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (19/25);  total (19/25)

CCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATG‑TTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGCTGATATCACGCAGGG  >  NC_000913/2031165‑2031434                                                                                                                                            |                                                                                                                                    ccGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATG‑t                                                                                                                                    >  1:1072192/1‑140 (MQ=255)           atCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATG‑tttttttcgca                                                                                                                          >  2:33173/1‑135 (MQ=255)           atCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATG‑tttttttcgca                                                                                                                          <  2:410228/140‑6 (MQ=255)           atCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATG‑tttttttcgca                                                                                                                          >  2:65410/1‑135 (MQ=255)            tCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCat                                                                                                                          <  2:627585/140‑1 (MQ=255)               cAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTACATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATc                                                                                                                       >  1:1012776/1‑140 (MQ=255)                aGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCATACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCg                                                                                                                      >  1:343416/1‑140 (MQ=255)                aGCCGCGATAGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCg                                                                                                                      >  1:945946/1‑140 (MQ=255)                 gCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGa                                                                                                                     >  2:901922/1‑140 (MQ=255)                 gCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGa                                                                                                                     <  1:33173/140‑1 (MQ=255)                 gCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGa                                                                                                                     >  1:593809/1‑140 (MQ=255)                 gCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATG‑tttttttcg                                                                                                                            <  1:245880/132‑4 (MQ=255)                 gCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATG‑tttttttcg                                                                                                                            >  2:245880/1‑129 (MQ=255)                   cgcgATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATa                                                                                                                   <  2:593809/140‑1 (MQ=255)                     cgatcgatTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATaga                                                                                                                 <  2:1072192/140‑1 (MQ=255)                     cgatcgatTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATaga                                                                                                                 <  1:901922/140‑1 (MQ=255)                       atcgatTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAgg                                                                                                               <  1:45449/140‑1 (MQ=255)                       atcgatTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAgg                                                                                                               >  1:1133762/1‑140 (MQ=255)                        tcgatTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCATACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGt                                                                                                              <  2:343416/140‑1 (MQ=255)                                     ccGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCAc                                                                                                 <  2:945946/140‑1 (MQ=255)                                          ccATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCg                                                                                            <  2:88422/140‑1 (MQ=255)                                          ccATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCg                                                                                            <  2:328588/140‑1 (MQ=255)                                           cATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGt                                                                                           <  2:255810/140‑1 (MQ=255)                                               tttGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAgg                                                                                       <  1:780518/140‑1 (MQ=255)                                                     ataaCGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGatatat                                                                                 <  2:522894/140‑1 (MQ=255)                                                        aCGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTc                                                                              >  1:641163/1‑140 (MQ=255)                                                         cgcgCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTc                                                                                          >  1:549997/1‑127 (MQ=255)                                                         cgcgCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTc                                                                                          <  2:549997/127‑1 (MQ=255)                                                         cgcgCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCg                                                                             >  2:543031/1‑140 (MQ=255)                                                           cgcAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGcc                                                                           >  1:838174/1‑140 (MQ=255)                                                                      tACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCg                                                                 <  1:859365/139‑1 (MQ=255)                                                                      tACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCg                                                                 >  2:859365/1‑139 (MQ=255)                                                                          cgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTcc                                                            <  1:782424/140‑1 (MQ=255)                                                                            tGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCaa                                                          >  2:1137468/1‑140 (MQ=255)                                                                                                      aCCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTgcg                                <  1:65410/140‑1 (MQ=255)                                                                                                            ccATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAGCAGCTGCGCCAGCGTCACTCGCTGCGCAGGGa                          >  2:1017217/1‑140 (MQ=255)                                                                                                            ccATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGa                          <  2:216424/140‑1 (MQ=255)                                                                                                                aaCCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACa                      >  1:417349/1‑140 (MQ=255)                                                                                                                aaCCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACa                      <  2:1087679/140‑1 (MQ=255)                                                                                                                 aCCGAAGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACAt                     <  2:1133762/140‑1 (MQ=255)                                                                                                                     aaGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTgcg                                <  1:775911/125‑1 (MQ=255)                                                                                                                     aaGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTgcg                                >  2:775911/1‑125 (MQ=255)                                                                                                                     aaGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGc                 <  2:417349/140‑1 (MQ=255)                                                                                                                     aaGCCGTTACCGCGCGCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGc                 >  1:566127/1‑140 (MQ=255)                                                                                                                               cgcgcgCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGCTGATATCACg       <  2:480630/140‑1 (MQ=255)                                                                                                                               cgcgcgCCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCCACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGCTGATATCACg       <  2:1012776/140‑1 (MQ=255)                                                                                                                                   cgcCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGCTGATATCACg       <  1:446500/136‑1 (MQ=255)                                                                                                                                   cgcCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGCTGATATCACg       >  2:446500/1‑136 (MQ=255)                                                                                                                                   cgcCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGCTGATATCACGCAgg   >  2:924276/1‑140 (MQ=255)                                                                                                                                    gcCTGATGTTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGCTGATATCACGCAggg  >  2:852546/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                    CCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATG‑TTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTCGCTGCGCAGGGAAACATTCGCTGATATCACGCAGGG  >  NC_000913/2031165‑2031434

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