Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,062,980	C→A	Q711K (CAG→AAG)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,062,980	0	C	A	100.0%	161.4 / NA	47	Q711K (CAG→AAG)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/24);  total (23/24)

ATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGC  >  NC_000913/3062846‑3063116                                                                                                                                       |                                                                                                                                         atatCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGcat                                                                                                                                     <  1:1178336/140‑1 (MQ=255) atatCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGcat                                                                                                                                     <  2:744950/140‑1 (MQ=255) atatCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGcat                                                                                                                                     <  1:279210/140‑1 (MQ=255) atatCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGcat                                                                                                                                     <  2:1208965/140‑1 (MQ=255) atatCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGcat                                                                                                                                     <  1:878431/140‑1 (MQ=255)        cGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGAtt                                                                                                                              >  1:663951/1‑140 (MQ=255)         gTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTa                                                                                                                             >  1:119633/1‑140 (MQ=255)           ggTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAAt                                                                                                                           >  2:1038745/1‑140 (MQ=255)            gTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATg                                                                                                                          >  1:1019889/1‑140 (MQ=255)                cgGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGc                                                                                                                      <  2:954746/140‑1 (MQ=255)                    tGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACg                                                                                                                  <  2:1000115/140‑1 (MQ=255)                      gCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCt                                                                                                                >  2:613163/1‑140 (MQ=255)                           aTTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAg                                                                                                           >  2:1142650/1‑140 (MQ=255)                                 ccgcAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATAAAGCCACGCTGGCAGAAAGTa                                                                                                     >  2:157104/1‑140 (MQ=255)                                         ttACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTcgccgc                                                                                             <  2:26205/140‑1 (MQ=255)                                             gCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTAc                                                                                         <  1:1038745/140‑1 (MQ=255)                                               cTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGt                                                                                       <  1:716958/140‑1 (MQ=255)                                               cTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGt                                                                                       <  2:1019889/140‑1 (MQ=255)                                                 tCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCa                                                                                     >  2:316717/1‑140 (MQ=255)                                                 tCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCa                                                                                     >  1:232963/1‑140 (MQ=255)                                                 tCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCa                                                                                     >  2:86891/1‑140 (MQ=255)                                                 tCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCa                                                                                     >  2:667964/1‑140 (MQ=255)                                                 tCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGAATTCGCCGCTTACGTCa                                                                                     >  1:131250/1‑140 (MQ=255)                                                 tCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACA‑‑GTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAgg                                                                                   >  1:976423/1‑140 (MQ=255)                                                        agagCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAg                                                                              <  1:239626/140‑1 (MQ=255)                                                            cgcgGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCCGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTg                                                                          >  2:1119180/1‑140 (MQ=255)                                                             gcggcgTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGc                                                                         <  2:1111373/140‑1 (MQ=255)                                                                         cgATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCa                                                             <  2:287731/140‑1 (MQ=255)                                                                                    ccAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGc                                                  <  1:950398/140‑1 (MQ=255)                                                                                            acCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGaa                                          >  2:878895/1‑140 (MQ=255)                                                                                             ccTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGaaa                                         >  1:119937/1‑140 (MQ=255)                                                                                                  gCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGt                                    >  2:857682/1‑140 (MQ=255)                                                                                                      gACAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCAcc                                >  2:866972/1‑140 (MQ=255)                                                                                                        cAGTGTGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCcg                              >  1:493584/1‑140 (MQ=255)                                                                                                          gtgtGCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCcgcg                            <  2:803317/140‑1 (MQ=255)                                                                                                           tgtgCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGt                           >  1:808207/1‑140 (MQ=255)                                                                                                             tgCGCGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGtca                         >  2:379556/1‑140 (MQ=255)                                                                                                              gcgcGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGtcat                        <  2:119937/140‑1 (MQ=255)                                                                                                              gcgcGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGtcat                        <  2:663951/140‑1 (MQ=255)                                                                                                              gcgcGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGtcat                        <  1:667964/140‑1 (MQ=255)                                                                                                              gcgcGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGtcat                        <  1:683888/140‑1 (MQ=255)                                                                                                              gcgcGACGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGtcat                        <  1:857682/140‑1 (MQ=255)                                                                                                                   aCGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGGCGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGc                   <  1:613163/140‑1 (MQ=255)                                                                                                                   aCGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGc                   <  2:232963/140‑1 (MQ=255)                                                                                                                    cGTACTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCa                  <  1:497162/140‑1 (MQ=255)                                                                                                                        cTGAATCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTaa              >  2:792261/1‑140 (MQ=255)                                                                                                                           aaTCTGATAAGCAAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTa           >  2:1005851/1‑140 (MQ=255)                                                                                                                                    agcaagcaTCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgctgc  <  1:792261/140‑1 (MQ=255)                                                                                                                                    agcaagcaTCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgctgc  <  1:86891/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         ATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGC  >  NC_000913/3062846‑3063116

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