Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,061,202	A→G	N118S (AAC→AGC)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,202	0	A	G	100.0%	139.5 / NA	40	N118S (AAC→AGC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/22);  total (18/22)

TGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTAC  >  NC_000913/3061066‑3061333                                                                                                                                         |                                                                                                                                    tGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGcc                                                                                                                                   >  2:373174/1‑139 (MQ=255)      cATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCcgcgc                                                                                                                              >  2:87293/1‑139 (MQ=255)        tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                            >  2:151184/1‑139 (MQ=255)        tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                            >  1:851358/1‑139 (MQ=255)        tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                            >  2:871863/1‑139 (MQ=255)          cGTCAGTATGCTGGTTTTTCAACAGCAAAATAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTgg                                                                                                                          <  1:522749/139‑1 (MQ=255)                   gCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACg                                                                                                                 >  1:361638/1‑139 (MQ=255)                          tttCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGGaa                                                                                                          >  1:301426/1‑139 (MQ=255)                                cAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCggg                                                                                                    <  2:851358/139‑1 (MQ=255)                                 aGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGc                                                                                                   >  1:337007/1‑139 (MQ=255)                                 aGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGc                                                                                                   >  2:371431/1‑139 (MQ=255)                                            ccAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGGAAAGCGGGCGTCGCTATCGa                                                                                        <  2:301426/139‑1 (MQ=255)                                             cAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGac                                                                                       >  1:103267/1‑139 (MQ=255)                                             cAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGac                                                                                       >  1:558406/1‑139 (MQ=255)                                              aaCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGGTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                      <  2:717263/139‑1 (MQ=255)                                              aaCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                      <  1:871863/139‑1 (MQ=255)                                              aaCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                      <  2:193269/139‑1 (MQ=255)                                              aaCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                      <  2:1102575/139‑1 (MQ=255)                                                 gCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCg                                                                                   <  1:229463/139‑1 (MQ=255)                                                   tttttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                 <  2:361638/139‑1 (MQ=255)                                                     tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGa                                                                               >  1:131036/1‑139 (MQ=255)                                                        aTCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGa                                                                            <  1:926712/139‑1 (MQ=255)                                                             cGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTc                                                                                                <  1:708512/114‑1 (MQ=255)                                                             cGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTc                                                                                                >  2:708512/1‑114 (MQ=255)                                                               tAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                  <  2:797352/110‑1 (MQ=255)                                                               tAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                  >  1:797352/1‑110 (MQ=255)                                                                       ccgccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                             <  1:373174/139‑1 (MQ=255)                                                                           cgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGAc                                                         >  2:566584/1‑139 (MQ=255)                                                                                   aaGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATccc                                                 <  1:151184/139‑1 (MQ=255)                                                                                    aGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                                <  1:497891/139‑1 (MQ=255)                                                                                                    tttGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCgg                                >  2:955240/1‑139 (MQ=255)                                                                                                          cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          <  1:561145/139‑1 (MQ=255)                                                                                                            tGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGTTATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          <  1:1029512/137‑1 (MQ=255)                                                                                                            tGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          >  2:1029512/1‑137 (MQ=255)                                                                                                                         ggCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgcg           <  1:87293/139‑1 (MQ=255)                                                                                                                           cTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCa         <  2:337007/139‑1 (MQ=255)                                                                                                                           cTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCa         <  2:131036/139‑1 (MQ=255)                                                                                                                                  tCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTAc  <  1:371431/139‑1 (MQ=255)                                                                                                                                    cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCACTGGATAAAATGTCGGTTTCGATGACCATGAATGgc            <  2:698809/127‑1 (MQ=255)                                                                                                                                    cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCACTGGATAAAATGTCGGTTTCGATGACCATGAATGgc            >  1:698809/1‑127 (MQ=255)                                                                                                                                         |                                                                                                                                    TGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTAC  >  NC_000913/3061066‑3061333

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