Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,032,339	A→G	intergenic (‑22/+45)	dcm ← / ← yedJ	DNA‑cytosine methyltransferase/putative HD superfamily phosphohydrolase YedJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,032,339	0	A	G	100.0%	84.9 / NA	25	intergenic (‑22/+45)	dcm/yedJ	DNA‑cytosine methyltransferase/putative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/14);  total (11/14)

CGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTT  >  NC_000913/2032213‑2032461                                                                                                                               |                                                                                                                           cgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCg                                                                                                                      >  1:389870/1‑133 (MQ=255) cgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCg                                                                                                                      >  1:504637/1‑133 (MQ=255)  gcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGa                                                                                                                     <  1:408976/133‑1 (MQ=255)  gcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGa                                                                                                                     <  2:110263/133‑1 (MQ=255)   ccaccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCCCCGACCTTCTGAGCTGGCCGAc                                                                                                                    >  1:220009/1‑133 (MQ=255)                aaCATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAAt                                                                                                       <  1:110977/133‑1 (MQ=255)                     catAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAg                                                                                                  >  2:373253/1‑133 (MQ=255)                     catAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCGGAGGTGGACGACAGGTTAACAAAATTACAg                                                                                                  >  1:167349/1‑133 (MQ=255)                             tGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                          >  1:22698/1‑133 (MQ=255)                                  cagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCa                                                                                     <  1:373253/133‑1 (MQ=255)                                                   ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                                                    >  1:355422/1‑133 (MQ=255)                                                   ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                                                    <  2:403686/133‑1 (MQ=255)                                                        gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                                               >  1:445240/1‑133 (MQ=255)                                                          ggCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAc                                                             <  1:423945/133‑1 (MQ=255)                                                          ggCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAc                                                             <  2:389870/133‑1 (MQ=255)                                                           gCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACg                                                            <  2:22698/133‑1 (MQ=255)                                                           gCATTCCCGGTGCTGCATGAATCGTTTAATGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACg                                                            <  2:144736/133‑1 (MQ=255)                                                             aTTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                                               >  2:412845/1‑128 (MQ=255)                                                             aTTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                                               <  1:412845/128‑1 (MQ=255)                                                              ttCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCAt                                                         <  2:336571/133‑1 (MQ=255)                                                              ttCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCAt                                                         <  2:355422/133‑1 (MQ=255)                                                                               aTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg                                        >  1:113861/1‑133 (MQ=255)                                                                                tCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa                                       <  2:315502/133‑1 (MQ=255)                                                                                    ttACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACAcc                                   >  2:389475/1‑133 (MQ=255)                                                                                                        aTAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCTTCTATCACCTTGTGATCGACCCCTTCATTCTCCCCCGCCAGtt               <  2:445240/133‑1 (MQ=255)                                                                                                            aTTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGc           >  1:492594/1‑133 (MQ=255)                                                                                                                     gCCTTCTGAGGTGGCATACGGGTTAACATAATTTCGGCCATTTAGGCCCACGCGCCCAGCCACGGAATGAAACCCCATCTATCACTGCGTGATCGACTCCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCtt  <  2:438662/133‑1 (MQ=255)                                                                                                                               |                                                                                                                           CGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTT  >  NC_000913/2032213‑2032461

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