Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,353,609	T→C	T148A (ACC→GCC)	yhcC ←	radical SAM family oxidoreductase YhcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,353,609	0	T	C	100.0%	65.7 / NA	21	T148A (ACC→GCC)	yhcC	radical SAM family oxidoreductase YhcC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/9);  total (12/9)

GCCGACAATTAAGTGGGAACACACCTTCAGCCCGCGCTGACGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGTTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCTGGCTC  >  NC_000913/3353490‑3353736                                                                                                                        |                                                                                                                                gcCGACAATTAAGTGGGAACACACCTTCAGCCCGCGCTGACGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGc                                                                                                                    >  2:460467/1‑133 (MQ=255)                       acCTTCAGCCCGCGCTGACGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTGGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCt                                                                                              <  2:435710/133‑1 (MQ=255)                        ccTTCAGCCCGCGCTGACGTGCCAGCTGGGTTGTACGCTGATAGCAGCCAAAAGCATGGCCGCGGTTGATGCGATGCAGTGTTCAGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTg                                                                                             <  2:275831/133‑1 (MQ=255)                        ccTTCAGCCCGCGCTGACGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTg                                                                                             <  2:283040/133‑1 (MQ=255)                                 cgcgCTGACGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTtata                                                                                    >  1:62581/1‑133 (MQ=255)                                 cgcgCTGACGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGCAGCCCCAGCTCCAGTCCCACTTCGTAACCCTGGTCCTtata                                                                                    >  2:224520/1‑133 (MQ=255)                                   cgcTGACGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATAtt                                                                                  >  2:47070/1‑133 (MQ=255)                                      tGACGTGCCAGCTGGTTTATACCCTGATAAAAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAGCCCTGGTCCTTATATTCGc                                                                               <  1:200587/133‑1 (MQ=255)                                         cGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCaaa                                                                            <  2:344253/133‑1 (MQ=255)                                            gCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCa                                                                         >  2:457787/1‑133 (MQ=255)                                                              tGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCg                                                       >  1:44921/1‑133 (MQ=255)                                                                aTAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGc                                                     >  1:307567/1‑133 (MQ=255)                                                                          aaaTCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATGTTCGCAAAGCAGTTCCGGCACCGCACCCGGGACGCAGTCcg                                           >  2:154710/1‑133 (MQ=255)                                                                            aTCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCcggg                                         <  1:47070/133‑1 (MQ=255)                                                                                tGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCcgggcggg                                     <  1:457787/133‑1 (MQ=255)                                                                                    cgcgGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTAcc                                 <  2:307567/133‑1 (MQ=255)                                                                                            atgcgatgcAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCcacacaca                         >  2:58773/1‑133 (MQ=255)                                                                                                 atgcAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACACAGTCCGGGCGGGTACCCACATACAAACCg                    >  1:192039/1‑133 (MQ=255)                                                                                                   gcAGTGTTTTGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGAc                  >  1:328634/1‑133 (MQ=255)                                                                                                            tGTCGTGGGCGGCTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGc         >  1:344732/1‑133 (MQ=255)                                                                                                            tGTCGTGGACGGCTTGTATCCCCAGCTCCAACCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGc         <  1:224520/133‑1 (MQ=255)                                                                                                                   ggcggcTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCAACACACAAACCTACTATATTGACCTGGCTc  >  1:239899/1‑133 (MQ=255)                                                                                                                        |                                                                                                                                GCCGACAATTAAGTGGGAACACACCTTCAGCCCGCGCTGACGTGCCAGCTGGGTTGTACGCTGATAACAGGCAAAATCATGGCCGCGGTTGATGCGATGCAGTGTTTTGTCGTGGGCGGTTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCTGGCTC  >  NC_000913/3353490‑3353736

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