Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,226,907	A→G	H1027R (CAC→CGC)	metH →	cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,226,907	0	A	G	100.0%	153.3 / NA	44	H1027R (CAC→CGC)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/20);  total (24/20)

GACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGA  >  NC_000913/4226771‑4227036                                                                                                                                         |                                                                                                                                  gACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGccgcca                                                                                                                                <  2:270515/140‑1 (MQ=255)           aTAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCaa                                                                                                                      >  2:463727/1‑140 (MQ=255)              aaTTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAg                                                                                                                   <  2:793931/140‑1 (MQ=255)                    gcgcCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGaa                                                                                                             >  2:188754/1‑140 (MQ=255)                    gcgcCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGaa                                                                                                             >  2:552535/1‑140 (MQ=255)                          agaAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACa                                                                                                       >  2:29260/1‑140 (MQ=255)                          agaAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACa                                                                                                       >  1:656461/1‑140 (MQ=255)                                      aTCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAac                                                                                           >  1:195312/1‑140 (MQ=255)                                          gcgtggcgtgGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAactact                                                                                       <  1:552535/140‑1 (MQ=255)                                           cgtggcgtggTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCt                                                                                                   >  2:690827/1‑127 (MQ=255)                                           cgtggcgtggTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCt                                                                                                   <  1:690827/127‑1 (MQ=255)                                            gtggcgtggTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGt                                                                                     <  1:337266/140‑1 (MQ=255)                                               gcgtggTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCt                                                                                                   <  2:475524/123‑1 (MQ=255)                                               gcgtggTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCt                                                                                                   >  1:475524/1‑123 (MQ=255)                                                         ccTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCg                                                                        >  2:245251/1‑140 (MQ=255)                                                              tCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGtc                                                                                    >  1:89210/1‑123 (MQ=255)                                                              tCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGtc                                                                                    <  2:89210/123‑1 (MQ=255)                                                                       aCCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCtt                                                          <  1:188754/140‑1 (MQ=255)                                                                       aCCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCtt                                                          >  2:199541/1‑140 (MQ=255)                                                                       aCCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCtt                                                          <  1:221651/140‑1 (MQ=255)                                                                       aCCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCtt                                                          <  1:245251/140‑1 (MQ=255)                                                                        ccGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCttt                                                         <  1:199541/140‑1 (MQ=255)                                                                             tgGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGt                                                    >  1:865680/1‑140 (MQ=255)                                                                              gggCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTa                                                   <  1:110443/140‑1 (MQ=255)                                                                                 cGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAg                                                <  2:512384/140‑1 (MQ=255)                                                                                       cATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGAtt                                          >  2:653374/1‑140 (MQ=255)                                                                                          tGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACa                                       >  1:80218/1‑140 (MQ=255)                                                                                          tGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACa                                       >  1:271609/1‑140 (MQ=255)                                                                                          tGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACa                                       >  1:690136/1‑140 (MQ=255)                                                                                                   ccGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCAt                              <  2:271609/140‑1 (MQ=255)                                                                                                     gTGACGAAACGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCAttt                            <  2:756523/140‑1 (MQ=255)                                                                                                           aaaCGCGTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCATCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTg                      >  2:381367/1‑140 (MQ=255)                                                                                                                cgTACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTACCTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTgg                 <  2:690136/140‑1 (MQ=255)                                                                                                                  tACCCATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCg               >  2:118003/1‑140 (MQ=255)                                                                                                                    cccATGTGATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCATCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCCGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCggg             <  1:381367/140‑1 (MQ=255)                                                                                                                         gtgATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGa        >  1:330771/1‑140 (MQ=255)                                                                                                                         gtgATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGa        >  1:487746/1‑140 (MQ=255)                                                                                                                         gtgATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGa        >  1:816286/1‑140 (MQ=255)                                                                                                                           gATCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAag      >  2:570368/1‑140 (MQ=255)                                                                                                                             tCAACGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAagag    <  1:570368/140‑1 (MQ=255)                                                                                                                               aaCGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGa        >  2:216993/1‑134 (MQ=255)                                                                                                                               aaCGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGa        <  1:216993/134‑1 (MQ=255)                                                                                                                               aaCGTCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGa  <  2:656461/140‑1 (MQ=255)                                                                                                                                   tCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAagag    >  2:522733/1‑134 (MQ=255)                                                                                                                                   tCAGCCGCCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAagag    <  1:522733/134‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  GACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGA  >  NC_000913/4226771‑4227036

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: ‑