Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,064,741	T→C	F250S (TTC→TCC)	scpB →	methylmalonyl‑CoA decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,064,741	0	T	C	100.0%	135.7 / NA	39	F250S (TTC→TCC)	scpB	methylmalonyl‑CoA decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/16);  total (23/16)

CCGATGAATTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAG  >  NC_000913/3064663‑3064813                                                                               |                                                                         ccGATGAATTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTccc                                                                        >  2:562690/1‑81 (MQ=255)  cGATGAATTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTccc                                                                        <  1:562646/80‑1 (MQ=255)   gATGAATTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTc                                                                      >  2:562686/1‑81 (MQ=255)   gATGAATTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTc                                                                      >  2:562694/1‑81 (MQ=255)    aTGAATTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCg                                                                     <  1:562647/81‑1 (MQ=255)       aaTTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGaa                                                                   >  2:562692/1‑80 (MQ=255)       aaTTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGaa                                                                   <  2:562648/80‑1 (MQ=255)        aTTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGaaaa                                                                 <  1:562650/81‑1 (MQ=255)        aTTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGaaaa                                                                 <  1:562649/81‑1 (MQ=255)        aTTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGaaa                                                                  >  1:562681/1‑80 (MQ=255)         tttGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGaaaaa                                                                >  2:562709/1‑81 (MQ=255)           tGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAAc                                                               <  1:562651/80‑1 (MQ=255)           tGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACg                                                              >  2:562700/1‑81 (MQ=255)                   ttCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCt                                                       >  1:562707/1‑80 (MQ=255)                      aGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAAtt                                                   >  1:562695/1‑81 (MQ=255)                      aGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAAtt                                                   >  1:562698/1‑81 (MQ=255)                         ggATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCg                                                <  2:562652/81‑1 (MQ=255)                           aTGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGtt                                              <  1:562653/81‑1 (MQ=255)                                  gcgcgGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCAt                                        >  2:562684/1‑80 (MQ=255)                                    gcgGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTaa                                     >  2:562703/1‑81 (MQ=255)                                     cgGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAAt                                    >  2:562705/1‑81 (MQ=255)                                     cgGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAAt                                    >  2:562710/1‑81 (MQ=255)                                      ggTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATc                                   >  1:562687/1‑81 (MQ=255)                                      ggTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATc                                   >  2:562715/1‑81 (MQ=255)                                        tgtATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATccc                                 >  1:562704/1‑81 (MQ=255)                                          tATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCt                                <  1:562655/80‑1 (MQ=255)                                          tATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCt                                <  1:562654/80‑1 (MQ=255)                                              aCAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGcgaa                           >  2:562689/1‑81 (MQ=255)                                               cAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGcgaac                          <  2:562656/81‑1 (MQ=255)                                               cAGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGcgaac                          >  2:562706/1‑81 (MQ=255)                                                aGCGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGcgaac                          >  1:562711/1‑80 (MQ=255)                                                  cGAAGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGcgaacgaa                       <  2:562657/81‑1 (MQ=255)                                                    aaGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAgg                     <  1:562658/81‑1 (MQ=255)                                                    aaGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAg                      >  2:562713/1‑80 (MQ=255)                                                     aGATTACCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAg                      >  2:562725/1‑79 (MQ=255)                                                          aCCAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTaa                >  2:562691/1‑80 (MQ=255)                                                            cAGGAAGGGATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTaaaaa             <  1:562659/81‑1 (MQ=255)                                                                 aGGGATGAACGCTTCACTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGaa        <  1:562660/81‑1 (MQ=255)                                                                   ggATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAAc      <  1:562661/81‑1 (MQ=255)                                                                    gATGAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACt     <  1:562662/81‑1 (MQ=255)                                                                       gAACGCTTCCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAg  >  1:562717/1‑81 (MQ=255)                                                                               |                                                                         CCGATGAATTTGAACGTATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAG  >  NC_000913/3064663‑3064813

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