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%	164.5 / NA	51	*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 (32/19);  total (32/19)

TAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACG  >  NC_000913/3352994‑3353255                                                                                                                                |                                                                                                                                       taataCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTgg                                                                                                                             >  1:274824/1‑139 (MQ=255)       tGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAg                                                                                                                       >  1:203465/1‑139 (MQ=255)         cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  2:128308/1‑139 (MQ=255)         cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  2:104973/1‑139 (MQ=255)         cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  2:585213/1‑139 (MQ=255)           ccAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCAc                                                                                                                   <  1:356817/139‑1 (MQ=255)              gcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           >  1:136531/1‑128 (MQ=255)              gcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           <  2:136531/128‑1 (MQ=255)              gcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                >  2:38635/1‑139 (MQ=255)                   cGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTcc                                                                                                           <  1:155370/139‑1 (MQ=255)                                  aGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTg                                                                                            <  1:104973/139‑1 (MQ=255)                                    ggAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTa                                                                                          >  1:551921/1‑139 (MQ=255)                                        ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGcc                                                                                      >  1:88396/1‑139 (MQ=255)                                         tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                     <  2:368679/139‑1 (MQ=255)                                         tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                     >  1:450433/1‑139 (MQ=255)                                           cTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                   <  2:88396/139‑1 (MQ=255)                                                  aCAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattc                                                                            <  2:497499/139‑1 (MQ=255)                                                     aCCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTg                                                                                            >  1:396584/1‑120 (MQ=255)                                                     aCCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTg                                                                                            <  2:396584/120‑1 (MQ=255)                                                                    aCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCg                                                          >  1:432726/1‑139 (MQ=255)                                                                     cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 >  1:65510/1‑115 (MQ=255)                                                                     cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 <  2:65510/115‑1 (MQ=255)                                                                         aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGAc                                                        <  1:172009/136‑1 (MQ=255)                                                                         aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGAc                                                        >  2:172009/1‑136 (MQ=255)                                                                         aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCAt                                                     >  2:250651/1‑139 (MQ=255)                                                                          ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                   >  2:687709/1‑108 (MQ=255)                                                                          ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                   <  1:687709/108‑1 (MQ=255)                                                                                cAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTc                                              >  1:245916/1‑139 (MQ=255)                                                                                    tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTg                                                                                            >  1:365091/1‑89 (MQ=255)                                                                                    tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTg                                                                                            <  2:365091/89‑1 (MQ=255)                                                                                    tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATc                                          >  2:582959/1‑139 (MQ=255)                                                                                       gTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCATGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                       <  2:551921/139‑1 (MQ=255)                                                                                           taAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTc                                   >  2:362686/1‑139 (MQ=255)                                                                                                 gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcacca                             >  2:477686/1‑139 (MQ=255)                                                                                                 gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcacca                             >  2:41401/1‑139 (MQ=255)                                                                                                 gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcacca                             <  1:126044/139‑1 (MQ=255)                                                                                                 gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcac                               <  1:62577/137‑1 (MQ=255)                                                                                                 gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcac                               >  2:62577/1‑137 (MQ=255)                                                                                                  tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAg                                                                      <  1:584685/97‑1 (MQ=255)                                                                                                  tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAg                                                                      >  2:584685/1‑97 (MQ=255)                                                                                                  tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                       >  2:134180/1‑128 (MQ=255)                                                                                                  tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                       <  1:134180/128‑1 (MQ=255)                                                                                                  tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcaccac                            >  2:226103/1‑139 (MQ=255)                                                                                                  tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcaccac                            >  2:205188/1‑139 (MQ=255)                                                                                                  tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcaccac                            >  2:534338/1‑139 (MQ=255)                                                                                                  tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCAACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcaccac                            >  2:575751/1‑139 (MQ=255)                                                                                                     gagaAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGc                         <  2:135429/139‑1 (MQ=255)                                                                                                            acGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCAAGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTGGCCCCACAGCGGagcaa                  >  1:473547/2‑139 (MQ=255)                                                                                                                  ttCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGAc                                                        >  2:356343/1‑95 (MQ=255)                                                                                                                  ttCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGAc                                                        <  1:356343/95‑1 (MQ=255)                                                                                                                            cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACg  >  1:151731/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                       TAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACG  >  NC_000913/3352994‑3353255

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