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

TCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTG  >  NC_000913/3061073‑3061329                                                                                                                                  |                                                                                                                                tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                        >  2:191355/1‑139 (MQ=255) tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                        >  1:47341/1‑139 (MQ=255)         tATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCcgcg                                                                                                                           <  2:341976/128‑1 (MQ=255)         tATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCcgcg                                                                                                                           >  1:341976/1‑128 (MQ=255)                gTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTGTTTCCGTTGCGTTTGACCTTCCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCgg                                                                                                         >  1:1285081/1‑139 (MQ=255)                gTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCgg                                                                                                         >  1:1191194/1‑139 (MQ=255)                                  agTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTAt                                                                                       <  2:1191194/139‑1 (MQ=255)                                     ccAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCGAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGa                                                                                    <  2:857901/139‑1 (MQ=255)                                     ccAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGa                                                                                    <  2:747925/139‑1 (MQ=255)                                      cAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGac                                                                                   >  2:399173/1‑139 (MQ=255)                                         cGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACAcc                                                                                <  2:260944/139‑1 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                              <  1:1040231/120‑1 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                              >  2:1040231/1‑120 (MQ=255)                                               ttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGGCGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGaa                                                                          <  2:915597/139‑1 (MQ=255)                                                   cGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                              <  1:311916/115‑1 (MQ=255)                                                   cGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                              >  2:311916/1‑115 (MQ=255)                                                        tAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGaaa                                                                 >  1:386958/1‑139 (MQ=255)                                                                ccgccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                             >  2:291579/1‑119 (MQ=255)                                                                ccgccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                             <  1:291579/119‑1 (MQ=255)                                                                   ccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGa                                                      >  2:504581/1‑139 (MQ=255)                                                                             aGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                            <  2:736221/139‑1 (MQ=255)                                                                                 cTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCcgcgcg                                                                                                                         <  2:889090/58‑1 (MQ=255)                                                                                 cTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCcgcgcg                                                                                                                         >  1:889090/1‑58 (MQ=255)                                                                                 cTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTgg                                        >  1:949444/1‑139 (MQ=255)                                                                                 cTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTgg                                        >  1:276136/1‑139 (MQ=255)                                                                                         tGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGc                                           <  1:732944/128‑1 (MQ=255)                                                                                         tGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGc                                           >  2:732944/1‑128 (MQ=255)                                                                                          gCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGt                               <  1:504581/139‑1 (MQ=255)                                                                                                gACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGa                                                                                    >  1:1242507/1‑80 (MQ=255)                                                                                                gACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGa                                                                                    <  2:1242507/80‑1 (MQ=255)                                                                                                   cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                      <  2:386958/139‑1 (MQ=255)                                                                                                   cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                      <  1:664203/139‑1 (MQ=255)                                                                                                       ccacccacCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTc                                                        <  1:907058/101‑1 (MQ=255)                                                                                                       ccacccacCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTc                                                        >  2:907058/1‑101 (MQ=255)                                                                                                             accGTGGCTACGACTCCTATAGCCCGCGCGTGGCGGGCGACGTCGGCAAACCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGaa            <  2:1285081/139‑1 (MQ=255)                                                                                                               cGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATg          >  1:579060/1‑139 (MQ=255)                                                                                                                gTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                              >  2:769504/1‑54 (MQ=255)                                                                                                                gTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                              <  1:769504/54‑1 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCg                                                                                           >  1:118105/1‑56 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCg                                                                                           <  2:118105/56‑1 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                             >  1:838648/1‑70 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                             <  2:838648/70‑1 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgc        >  1:993193/1‑139 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgc        <  1:399173/139‑1 (MQ=255)                                                                                                                  ggCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGa                                                                        <  1:239747/74‑1 (MQ=255)                                                                                                                  ggCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGa                                                                        >  2:239747/1‑74 (MQ=255)                                                                                                                   gCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgcgc      <  1:440400/139‑1 (MQ=255)                                                                                                                       cGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGAt                                                >  1:1244713/1‑93 (MQ=255)                                                                                                                       cGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGAt                                                <  2:1244713/93‑1 (MQ=255)                                                                                                                       cGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTg  <  2:579060/139‑1 (MQ=255)                                                                                                                         aCTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGAc                                                     <  2:593990/86‑1 (MQ=255)                                                                                                                         aCTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGAc                                                     <  2:960330/86‑1 (MQ=255)                                                                                                                         aCTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGAc                                                     >  1:593990/1‑86 (MQ=255)                                                                                                                         aCTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGAc                                                     >  1:960330/1‑86 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                              >  1:528427/1‑41 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                              <  2:528427/41‑1 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                            >  1:154213/1‑91 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                            <  2:154213/91‑1 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATg                    >  1:1313745/1‑115 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATg                    <  2:1313745/115‑1 (MQ=255)                                                                                                                                  |                                                                                                                                TCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTG  >  NC_000913/3061073‑3061329

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