Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,227,425	G→A	V1200I (GTT→ATT)	metH →	cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,227,425	0	G	A	100.0%	84.3 / NA	25	V1200I (GTT→ATT)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/1);  total (24/1)

GCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAG  >  NC_000913/4227289‑4227559                                                                                                                                         |                                                                                                                                       gctgGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTg                                                                                                                                     >  1:247134/1‑140 (MQ=255) gctgGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTg                                                                                                                                     <  1:626533/140‑1 (MQ=255) gctgGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTg                                                                                                                                     >  1:61467/1‑140 (MQ=255) gctgGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTg                                                                                                                                     >  1:296917/1‑140 (MQ=255)  ctgGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGa                                                                                                                                    >  1:464839/1‑140 (MQ=255)                    acTGGCATAAAACTCACAGAATCTTTCGCCATGTGGCCGGGTGCATGGGTTGGGGGTTGGTACTTCAGCCCCCCGGACAGAAAGTCCTACGCTGTAGCACAAATTCAGCGCGATAAGATTGAAGATTATGCCCGCCGTaa                                                                                                                  <  2:570450/140‑1 (MQ=255)                         cATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTa                                                                                                             >  2:696089/1‑140 (MQ=255)                                cTCACAGATTCTTCCGCCGTGTGGCCCGGTGCATTGGTTTGGGTTTGGTCCTCCCGCCCCCGGGACAGCAATTACTCCGTTGTAGAACAAATTCGGCGCGATAAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGt                                                                                                      <  2:509759/140‑1 (MQ=255)                                             tCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAg                                                                                         >  2:261767/1‑140 (MQ=255)                                               gCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAgcg                                                                                       >  1:721215/1‑140 (MQ=255)                                                   tgtgGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCgctgg                                                                                   >  2:442785/1‑140 (MQ=255)                                                             tGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGa                                                                         >  1:291210/1‑140 (MQ=255)                                                             tGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGa                                                                         >  1:283553/1‑140 (MQ=255)                                                             ggCATGGTTTTGGGGTTGGTCCTCCACCCCCCCGGACAGAAAGTACTACGCTGTAGCACAAATCCAGCGCGATAAGATTGAAGATTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGa                                                                         <  2:543819/139‑1 (MQ=255)                                                                      tggggtttGGTACTCCACCCCCCCGGACAGAAAGTCCTCCGCTGTAGCAAAAATTCAGCGCGATCAGATTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGt                                                                <  2:61467/134‑1 (MQ=255)                                                                       tCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTa                                                               >  1:53744/1‑140 (MQ=255)                                                                              ggTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGAcgcg                                                        >  1:662982/1‑140 (MQ=255)                                                                                   ttCACCCCCCCGGACAGCAAGTACTACGCTGTAGCAAAAATCCAGCGCAATAAGTTTGAAGTTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTg                                                   <  2:662982/140‑1 (MQ=255)                                                                                         cACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcac                                             >  1:285917/1‑140 (MQ=255)                                                                                          aCCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcaca                                            >  2:513146/1‑140 (MQ=255)                                                                                                    cAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCAc                                  >  1:143713/1‑140 (MQ=255)                                                                                                      aGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtt                                >  2:425852/1‑140 (MQ=255)                                                                                                      aGTACTACGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtt                                >  1:365455/1‑140 (MQ=255)                                                                                                          gacagCTGTAGCAAAATTTCGGCGCGATCAGATTGAAGTTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                            <  2:728439/136‑1 (MQ=255)                                                                                                          ctccgttGTAACAAAAATTCAGCGCAATCAGATAGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                            <  2:476587/134‑1 (MQ=255)                                                                                                          ctacGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                            <  1:696089/140‑1 (MQ=255)                                                                                                          ctacGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                            <  1:442785/140‑1 (MQ=255)                                                                                                          ctacGCTGTACCAAAAATCCAGCGCGAAAAGATTAAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                            <  2:53744/140‑1 (MQ=255)                                                                                                             cGCTGTAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTa                         >  1:204270/1‑140 (MQ=255)                                                                                                                  tAGCACAAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaaca                    >  1:439486/1‑140 (MQ=255)                                                                                                                     cacaAATTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                 >  1:511820/1‑140 (MQ=255)                                                                                                                        aaaTTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg              <  2:143713/140‑1 (MQ=255)                                                                                                                         aaTCCAGCGCGATAAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAggg             <  2:464839/140‑1 (MQ=255)                                                                                                                          aTTCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaacctgtc            >  1:728439/1‑135 (MQ=255)                                                                                                                            tCAGCGCGATCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc          >  1:375903/1‑140 (MQ=255)                                                                                                                                    ttCAGATTGAAGTTTATACCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg  <  2:285917/139‑1 (MQ=255)                                                                                                                                    aTCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg  >  2:373467/1‑140 (MQ=255)                                                                                                                                    aTCAGATTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTATTCCTTACAACAAACAGGGCGCTCAATGAg  >  2:420611/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       GCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAG  >  NC_000913/4227289‑4227559

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