Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,062,332	G→A	A495T (GCC→ACC)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,062,332	0	G	A	100.0%	78.1 / NA	23	A495T (GCC→ACC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/0);  total (23/0)

CGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTGGGTGAAATTTCCGATGCGCTGGAAGTCG  >  NC_000913/3062200‑3062461                                                                                                                                     |                                                                                                                                  cGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGt                                                                                                                            >  1:1017499/1‑140 (MQ=255) cGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGt                                                                                                                            >  1:967812/1‑140 (MQ=255)     tCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAAcc                                                                                                                        >  2:330676/1‑140 (MQ=255)     tCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAAcc                                                                                                                        >  1:276401/1‑140 (MQ=255)     tCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAAcc                                                                                                                        >  1:859542/1‑140 (MQ=255)               tgtCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGa                                                                                                              >  1:217711/1‑140 (MQ=255)               tgtCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGa                                                                                                              >  1:297223/1‑140 (MQ=255)                 tCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAAc                                                                                                            >  2:157523/1‑140 (MQ=255)                   aacaaGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGc                                                                                                          >  2:736119/1‑140 (MQ=255)                     caaGTACAAACTGGATCACGAAGACGAAACCGATGTATTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGccc                                                                                                        >  1:703324/1‑140 (MQ=255)                            aaaCGGGACCAAAAAGACAAACCCGTTGTTTTGGAGACCGAAAGGGTGGGGGTGCGAAAGGAAAAATTTGCTCCGCGGAAACGCATCCGCCCCCCCGGTGATGATCCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCa                                                                                                 <  2:823362/140‑1 (MQ=255)                                  gATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcgccgc                                                                                           >  1:257403/1‑140 (MQ=255)                                        gAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCa                                                                                     >  1:499894/1‑140 (MQ=255)                                                                aTCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGAAACGCATTCGCGCCACCCGTGAGGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGc                                                             >  1:926604/1‑140 (MQ=255)                                                                    aCAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGtt                                                         >  2:632244/1‑140 (MQ=255)                                                                        cGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATg                                                     >  2:541993/1‑140 (MQ=255)                                                                        cGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATg                                                     >  2:63114/1‑140 (MQ=255)                                                                         gTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGc                                                    >  2:103556/1‑140 (MQ=255)                                                                          tGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGcc                                                   >  1:263295/1‑140 (MQ=255)                                                                             gggTGCGAAACGACAAATTTGCTCCGCGGAACCGCTTCCGCCCCCCCGGTGATGATCCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCt                                                <  1:309919/139‑1 (MQ=255)                                                                                    tAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTc                                         >  2:406502/1‑140 (MQ=255)                                                                                          gCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCa                                   >  2:11662/1‑140 (MQ=255)                                                                                          gCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCa                                   >  1:352415/1‑140 (MQ=255)                                                                                                  cTTCGCTGGAACGCATTCGCGCCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTGGGt                           >  2:468/1‑140 (MQ=255)                                                                                                                     cgcCACCCGTGATGATACCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTGGGTGAAATTTCCGATGCGCTgg        >  2:322226/1‑140 (MQ=255)                                                                                                                           cgggaatgatcCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTGGGTGAAATTTCCGATGCGCTGGAAGTCg  <  1:736119/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  CGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTGGGTGAAATTTCCGATGCGCTGGAAGTCG  >  NC_000913/3062200‑3062461

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