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%	151.6 / NA	43	N118S (AAC→AGC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/23);  total (20/23)

TCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTAC  >  NC_000913/3061073‑3061333                                                                                                                                  |                                                                                                                                    tCCGTCAGTATGCTGGTTTTTCAACAGCCAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                            >  2:498331/1‑139 (MQ=255) tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCGACCCACCGTGGCTACGACTCTGATAGCCCGCGCGt                                                                                                                            >  1:714363/1‑139 (MQ=255) tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                            >  1:312793/1‑139 (MQ=255) tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                            >  1:84289/1‑139 (MQ=255)    gTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTggt                                                                                                                         >  2:359947/1‑138 (MQ=255)            gCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACg                                                                                                                 >  1:520944/1‑139 (MQ=255)                     tCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAg                                                                                                        <  2:312793/139‑1 (MQ=255)                          aGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGc                                                                                                   >  2:402461/1‑139 (MQ=255)                           gCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                  <  2:179998/139‑1 (MQ=255)                             aaaaGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTc                                                                                                <  2:551767/139‑1 (MQ=255)                                     ccAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGa                                                                                        <  2:289360/139‑1 (MQ=255)                                         cGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACAcc                                                                                    >  2:571121/1‑139 (MQ=255)                                            tttttATCGCCATAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                 <  2:520944/139‑1 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                 <  2:757014/137‑1 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                 >  1:757014/1‑137 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGa                                                                               <  1:622873/139‑1 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGa                                                                               >  2:416585/1‑139 (MQ=255)                                                      cGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGa                                                                       >  1:278136/1‑139 (MQ=255)                                                      cGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCGACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGa                                                                       <  1:371356/139‑1 (MQ=255)                                                                ccgccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                             <  1:571121/139‑1 (MQ=255)                                                                      ggCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCg                                                                                               >  2:232520/1‑99 (MQ=255)                                                                      ggCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCg                                                                                               <  1:232520/99‑1 (MQ=255)                                                                                         tGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                 >  1:666936/1‑94 (MQ=255)                                                                                         tGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                 <  2:666936/94‑1 (MQ=255)                                                                                          gCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGt                                   >  2:610336/1‑139 (MQ=255)                                                                                            gTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGacac                                                                                     <  2:601407/87‑1 (MQ=255)                                                                                            gTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGacac                                                                                     >  1:601407/1‑87 (MQ=255)                                                                                                   cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGTGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          <  1:498331/139‑1 (MQ=255)                                                                                                   cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGTGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          <  1:359947/139‑1 (MQ=255)                                                                                                   cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          <  2:84289/139‑1 (MQ=255)                                                                                                   cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          <  2:177306/139‑1 (MQ=255)                                                                                                     tGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                  <  1:670690/65‑1 (MQ=255)                                                                                                     tGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                  >  2:670690/1‑65 (MQ=255)                                                                                                     tGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATg                        >  2:771629/1‑139 (MQ=255)                                                                                                               cGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATg                                                                        >  2:471614/1‑81 (MQ=255)                                                                                                               cGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATg                                                                        <  1:471614/81‑1 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          <  2:528527/125‑1 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                          >  1:528527/1‑125 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgc            <  1:416585/139‑1 (MQ=255)                                                                                                                  ggCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgcg           <  1:610336/139‑1 (MQ=255)                                                                                                                    cTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgcg           >  2:304461/1‑137 (MQ=255)                                                                                                                    cTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgcg           <  1:304461/137‑1 (MQ=255)                                                                                                                           tCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTAc  <  1:701353/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                    TCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTAC  >  NC_000913/3061073‑3061333

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