Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,226,924	A→G	T1033A (ACC→GCC)	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,924	0	A	G	100.0%	146.8 / NA	41	T1033A (ACC→GCC)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/21);  total (20/21)

GAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAG  >  NC_000913/4226795‑4227059                                                                                                                                  |                                                                                                                                        gagaAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAAc                                                                                                                               >  1:264977/1‑140 (MQ=255)  agaAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACa                                                                                                                              >  2:530661/1‑140 (MQ=255)  agaAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACa                                                                                                                              >  2:568603/1‑140 (MQ=255)      aaCGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCt                                                                                                                          <  1:694901/140‑1 (MQ=255)      aaCGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCt                                                                                                                          <  1:618532/140‑1 (MQ=255)         gCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCg                                                                                                                       <  2:264977/140‑1 (MQ=255)              aTCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAac                                                                                                                  <  2:678943/140‑1 (MQ=255)                       gcgtggTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGtctc                                                                                                         <  1:530661/140‑1 (MQ=255)                                       cccGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTgcgc                                                                                         >  1:806598/1‑140 (MQ=255)                                               aCCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCtt                                                                                 >  1:555385/1‑140 (MQ=255)                                               aCCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCtt                                                                                 <  1:390281/140‑1 (MQ=255)                                               aCCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCtt                                                                                 >  1:865589/1‑140 (MQ=255)                                                    gtgGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTgg                                                                            <  1:507246/140‑1 (MQ=255)                                                    gtgGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTgg                                                                            >  1:607803/1‑140 (MQ=255)                                                     tgGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGt                                                                           >  1:748856/1‑140 (MQ=255)                                                             gACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTaaa                                                                        >  1:903696/1‑135 (MQ=255)                                                             gACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTaaa                                                                        <  2:903696/135‑1 (MQ=255)                                                             gACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGa                                                                   >  2:488278/1‑140 (MQ=255)                                                                          aCCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAACCAGATTACATCGGCGCa                                                      <  2:748856/140‑1 (MQ=255)                                                                            cGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCAtt                                                    >  2:113159/1‑140 (MQ=255)                                                                                    aaCGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTAGCTAACTACTGTCACACTGACTTCGTTGCGGCGAAGCTTTCTGGTAAAGCAGATTACACCGGCGCATTTGCCGTGa                                            <  2:375429/140‑1 (MQ=255)                                                                                            cccATGTGATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCggg                                    <  1:488278/140‑1 (MQ=255)                                                                                                 gtgATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGc                                   <  2:581334/136‑1 (MQ=255)                                                                                                 gtgATCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGc                                   >  1:581334/1‑136 (MQ=255)                                                                                                    aTCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAaga                            >  1:449221/1‑140 (MQ=255)                                                                                                     tCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAagag                           <  1:113159/140‑1 (MQ=255)                                                                                                     tCAACGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAagag                           <  2:607803/140‑1 (MQ=255)                                                                                                       aaCGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCg                                      <  1:129242/127‑1 (MQ=255)                                                                                                       aaCGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCg                                      >  2:129242/1‑127 (MQ=255)                                                                                                         cGTCAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACg                       >  2:196905/1‑140 (MQ=255)                                                                                                            cAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCAc                    >  1:55729/1‑140 (MQ=255)                                                                                                            cAGCCACCATCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCAc                    <  1:196905/140‑1 (MQ=255)                                                                                                               ccaccaTCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGc                                   <  1:97842/122‑1 (MQ=255)                                                                                                               ccaccaTCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGc                                   >  2:97842/1‑122 (MQ=255)                                                                                                               ccaccaTCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCa                     >  1:86522/1‑136 (MQ=255)                                                                                                               ccaccaTCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCa                     <  2:86522/136‑1 (MQ=255)                                                                                                               ccaccaTCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCAc                    <  2:21693/137‑1 (MQ=255)                                                                                                               ccaccaTCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCAc                    >  1:21693/1‑137 (MQ=255)                                                                                                                   caTCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGa             <  1:705202/140‑1 (MQ=255)                                                                                                                   caTCTGCGTCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGa             >  2:809087/1‑140 (MQ=255)                                                                                                                           tCAACAGGCCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTg     <  1:809087/140‑1 (MQ=255)                                                                                                                               cAGGCCGGAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCTGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAg  <  2:873660/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        GAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAG  >  NC_000913/4226795‑4227059

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