Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,353,121	C→T	*310* (TAG→TAA)	yhcC ←	putative Fe‑S oxidoreductase, Radical SAM superfamily protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,353,121	0	C	T	100.0%	135.6 / NA	44	*310* (TAG→TAA)	yhcC	putative Fe‑S oxidoreductase, Radical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/21);  total (23/21)

TCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACG  >  NC_000913/3352988‑3353255                                                                                                                                      |                                                                                                                                       tCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACATCCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    <  2:411660/138‑1 (MQ=255) tCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    >  1:411660/1‑138 (MQ=255)          taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTgg                                                                                                                             <  1:114650/136‑1 (MQ=255)          taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTgg                                                                                                                             >  2:114650/1‑136 (MQ=255)          taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                          >  1:348080/1‑139 (MQ=255)          taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                          >  1:428465/1‑139 (MQ=255)               cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  2:438506/1‑139 (MQ=255)               cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  2:416014/1‑139 (MQ=255)               cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  2:127479/1‑139 (MQ=255)                  cagcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACg                                                                                                                  >  2:835133/1‑139 (MQ=255)                    gcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                >  2:397571/1‑139 (MQ=255)                         cGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTcc                                                                                                           <  1:438506/139‑1 (MQ=255)                              ttGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAgcg                                                                                                      <  1:835133/139‑1 (MQ=255)                                      ttAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATccc                                                                                              <  1:397571/139‑1 (MQ=255)                                         gggAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGt                                                                                           <  1:416014/139‑1 (MQ=255)                                             aTTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGc                                                                                       <  2:348080/139‑1 (MQ=255)                                              ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                <  2:57087/113‑1 (MQ=255)                                              ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                >  1:57087/1‑113 (MQ=255)                                               tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                <  1:429317/112‑1 (MQ=255)                                               tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                >  2:429317/1‑112 (MQ=255)                                               tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                     <  2:98645/139‑1 (MQ=255)                                                ccTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCAt                                                                                    >  1:28380/1‑139 (MQ=255)                                                 cTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCGTCAACGACTTACTCCGTTGGAGTTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                   >  2:380118/1‑139 (MQ=255)                                                      cGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcat                                                                              <  2:28380/139‑1 (MQ=255)                                                           aCCTGCCCCAATGATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGa                                                                         >  1:500966/1‑139 (MQ=255)                                                            ccTGCCCGAATAATACCCAATTCTGACAGTTACGTAATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTc                                                                                                                                   <  2:718061/80‑1 (MQ=255)                                                            ccTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTc                                                                                                                                   >  1:718061/1‑80 (MQ=255)                                                                 cccAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGa                                                                   <  1:127479/139‑1 (MQ=255)                                                                   cAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATc                                                                 <  1:321163/139‑1 (MQ=255)                                                                     ataataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCa                                                               >  2:526861/1‑139 (MQ=255)                                                                      taataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATccc                                                                                              >  1:463670/1‑107 (MQ=255)                                                                      taataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATccc                                                                                              <  2:463670/107‑1 (MQ=255)                                                                               aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    <  2:20108/60‑1 (MQ=255)                                                                               aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    >  1:20108/1‑60 (MQ=255)                                                                                                       gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGc                                                                                       <  1:608245/81‑1 (MQ=255)                                                                                                       gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGc                                                                                       >  2:608245/1‑81 (MQ=255)                                                                                                        tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATa                                                                       >  1:211987/1‑96 (MQ=255)                                                                                                        tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATa                                                                       <  2:211987/96‑1 (MQ=255)                                                                                                             gaAATGCGGAATTCCCCTCAACGACTTACTCCGTTGAAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCgg                       <  1:380118/139‑1 (MQ=255)                                                                                                                                  cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 >  1:229352/1‑60 (MQ=255)                                                                                                                                  cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 <  2:557317/60‑1 (MQ=255)                                                                                                                                  cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 <  2:229352/60‑1 (MQ=255)                                                                                                                                  cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 >  1:557317/1‑60 (MQ=255)                                                                                                                                  cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACg  >  1:448748/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                       TCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACG  >  NC_000913/3352988‑3353255

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