Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,062,645	G→A	R599H (CGC→CAC)	scpA →	methylmalonyl‑CoA mutase

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

CAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTG  >  NC_000913/3062512‑3062781                                                                                                                                      |                                                                                                                                         cAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGgcgc                                                                                                                                    >  1:106605/1‑140 (MQ=255) cAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGgcgc                                                                                                                                    >  2:359597/1‑140 (MQ=255)         tcatcaGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGa                                                                                                                            <  1:509506/140‑1 (MQ=255)                    agaAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGcc                                                                                                                 <  1:169671/140‑1 (MQ=255)                              cTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGAtc                                                                                                       >  1:340918/1‑140 (MQ=255)                                     ttCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGttt                                                                                                <  1:442821/140‑1 (MQ=255)                                     ttCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGttt                                                                                                <  1:92306/140‑1 (MQ=255)                                             cATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAg                                                                                        >  1:495675/1‑140 (MQ=255)                                                tgttgCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGAtt                                                                                     <  2:495675/140‑1 (MQ=255)                                                     cgcAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGc                                                                                <  1:2092/140‑1 (MQ=255)                                                           cGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATg                                                                          <  1:328925/140‑1 (MQ=255)                                                           cGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATg                                                                          <  1:190440/140‑1 (MQ=255)                                                           cGGAGCAGTTCCTTGCCGACAATGGTCGTCGACCGCGCATTCTGATCGCTAAGATGGGCCAGGAGGGACACGATCATGGCGCGAAAGTGATCGCGAGCGTCTATTCCGGTCTCGGGTTCGACGTAGATTTAAGCCCGATg                                                                          <  2:424994/140‑1 (MQ=255)                                                               gCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCg                                                                             >  2:66354/1‑133 (MQ=255)                                                               gCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCg                                                                             <  1:66354/133‑1 (MQ=255)                                                                cAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTtctc                                                                     <  2:67324/140‑1 (MQ=255)                                                                cAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTtctc                                                                     <  2:310317/140‑1 (MQ=255)                                                                cAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTtctc                                                                     <  2:255147/140‑1 (MQ=255)                                                                cAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTtctc                                                                     <  1:277606/140‑1 (MQ=255)                                                                cAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTtctc                                                                     <  1:359597/140‑1 (MQ=255)                                                                  gTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTa                                                                   <  2:459744/140‑1 (MQ=255)                                                                        tGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATggt                                                                        >  1:278362/1‑127 (MQ=255)                                                                        tGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATg                                                                          >  1:63598/1‑127 (MQ=255)                                                                        tGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATg                                                                          <  2:63598/127‑1 (MQ=255)                                                                        tGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGtt                                                                        <  2:278362/129‑1 (MQ=255)                                                                             aCAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAagag                                                        >  1:486915/1‑140 (MQ=255)                                                                               aaTGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGAt                                                      >  1:458955/1‑140 (MQ=255)                                                                                         gCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTg                                            >  2:316348/1‑140 (MQ=255)                                                                                          cccGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTgg                                           <  2:340918/140‑1 (MQ=255)                                                                                                   tCTGATCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGaaa                                  <  2:298930/140‑1 (MQ=255)                                                                                                       aTCGCTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAacga                              <  2:486915/140‑1 (MQ=255)                                                                                                           cTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGtt                          <  2:458955/140‑1 (MQ=255)                                                                                                           cTAAGATGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGtt                          >  1:119424/1‑140 (MQ=255)                                                                                                                 tGGGCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTa                    <  2:119424/140‑1 (MQ=255)                                                                                                                   ggCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGt                  >  2:283264/1‑140 (MQ=255)                                                                                                                    gCCAGGATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTg                 >  1:6166/1‑140 (MQ=255)                                                                                                                         gATGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGgcgc            >  1:358182/1‑140 (MQ=255)                                                                                                                          aTGGACACGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCa           <  2:292601/140‑1 (MQ=255)                                                                                                                              acacGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCt       >  1:277369/1‑140 (MQ=255)                                                                                                                                acGATCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCa     <  2:6166/140‑1 (MQ=255)                                                                                                                                   aTCACGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTg  >  1:519151/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                         CAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTG  >  NC_000913/3062512‑3062781

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