Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,353,378	C→G	A225P (GCG→CCG)	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,378	0	C	G	100.0%	145.1 / NA	40	A225P (GCG→CCG)	yhcC	radical SAM family oxidoreductase YhcC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (20/20);  total (20/20)

GCGAATCATTTCTCCGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGCCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGTTTCAACCACTCGCT  >  NC_000913/3353301‑3353449                                                                              |                                                                        gCGAATCATTTCTCCGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCa                                                                       >  1:626957/1‑80 (MQ=255)  cGAATCATTTCTCCGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCAt                                                                      >  2:626940/1‑80 (MQ=255)   gAATCATTTCTCCGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATaa                                                                    >  1:626985/1‑81 (MQ=255)    aaTCATTTCTCCGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAAt                                                                   <  2:626916/81‑1 (MQ=255)     aTCATTTCTCCGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAAtg                                                                  <  2:626917/81‑1 (MQ=255)           tctcCGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGcc                                                             <  1:626918/80‑1 (MQ=255)              ccGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCttt                                                          >  1:626973/1‑80 (MQ=255)              ccGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCttt                                                          >  1:626967/1‑80 (MQ=255)               cGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTca                                                        >  1:626964/1‑81 (MQ=255)               cGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTc                                                         <  1:626919/80‑1 (MQ=255)               aggCAGTGAGCGTGTAATACACCAGTTCAATACCGTTCAAACGCCACGCTTCCCAGGCTTTCGGCATAATGCTGCCttt                                                          <  1:626921/78‑1 (MQ=255)                 gCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTcaca                                                      <  2:626920/81‑1 (MQ=255)                   aGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAat                                                    <  2:626922/81‑1 (MQ=255)                    gTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAata                                                   <  2:626923/81‑1 (MQ=255)                       aGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATg                                                 >  1:626971/1‑80 (MQ=255)                       aGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATg                                                 <  1:626924/80‑1 (MQ=255)                        gCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATg                                                 >  2:626974/1‑79 (MQ=255)                         cGTGTAATCCTCCAGTTCGATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAg                                              <  2:626925/81‑1 (MQ=255)                          gtgtAATCCTCCAGTTCAATACCGTTCAAACACCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGc                                             <  2:626926/81‑1 (MQ=255)                            gtAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCg                                            >  1:626970/1‑80 (MQ=255)                               aTCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATGTGCAGCGGATg                                        <  1:626927/81‑1 (MQ=255)                                  ctccAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAg                                     >  2:626979/1‑81 (MQ=255)                                     cAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCtt                                  <  1:626928/81‑1 (MQ=255)                                     cAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCtt                                  >  2:626980/1‑81 (MQ=255)                                           aaTACCGTTCAGACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGGTGCAGCTTTATGcc                            <  1:626929/81‑1 (MQ=255)                                              aCCGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTc                         >  1:626978/1‑81 (MQ=255)                                               ccGTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTc                         <  1:626930/80‑1 (MQ=255)                                                 gTTCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCAc                      <  2:626931/81‑1 (MQ=255)                                                   tCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGc                    >  2:626968/1‑81 (MQ=255)                                                   tCAAACGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGc                    >  2:626976/1‑81 (MQ=255)                                                     aaaCGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCg                  >  1:626996/1‑81 (MQ=255)                                                      aaCGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCgg                 >  2:626982/1‑81 (MQ=255)                                                        cGCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGtt               <  2:626932/81‑1 (MQ=255)                                                         gCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGttt              >  1:626989/1‑81 (MQ=255)                                                         gCCCCGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGttt              >  1:626988/1‑81 (MQ=255)                                                           cccGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGTTTc             >  2:626987/1‑80 (MQ=255)                                                            ccGCTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGTTTCaa           <  2:626933/81‑1 (MQ=255)                                                               cTTCCCAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGTTTCAACCa        <  2:626934/81‑1 (MQ=255)                                                                    cAGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGTTTCAACCACTCGc   <  2:626935/81‑1 (MQ=255)                                                                     aGGCTTTCGGCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGTTTCAACCACTCGCt  >  2:626990/1‑81 (MQ=255)                                                                              |                                                                        GCGAATCATTTCTCCGGCAGTGAGCGTGTAATCCTCCAGTTCAATACCGTTCAAACGCCCCGCTTCCCAGGCTTTCGCCATAATGCTGCCTTTCACAATATGCAGCGGATGCAGCTTTATGCCGTCCACGCCGGTTTCAACCACTCGCT  >  NC_000913/3353301‑3353449

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