Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,658,730	Δ1 bp	coding (662/2427 nt)	ynfE →	putative selenate reductase YnfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,658,727	0	A	.	100.0%	75.9 / NA	19	coding (659/2427 nt)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (9/10);  total (9/10)

CTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTAT  >  NC_000913/1658650‑1658800                                                                              |                                                                          cTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGCCAACAGCACCACGGATATTGaa                                                                          <  1:304091/79‑1 (MQ=255)   aCAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGCCAACAGCACCACGGATATTGaaacag                                                                      >  2:304114/1‑78 (MQ=255)             gcAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTc                                                           <  1:304092/81‑1 (MQ=255)             gcAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTc                                                           >  2:304118/1‑81 (MQ=255)              cAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCg                                                          <  2:304093/81‑1 (MQ=255)                   ttCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATg                                                     >  2:304112/1‑81 (MQ=255)                     ccTGTGCCATGCCCTACACCTACGGCAGTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGtt                                                   >  1:304115/1‑81 (MQ=255)                                       ccTACGGCAGTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAg                                  >  2:304121/1‑80 (MQ=255)                                       ccTACGGCAGTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAg                                  >  2:304132/1‑80 (MQ=255)                                           cGGCAGTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAAcc                             <  1:304094/81‑1 (MQ=255)                                                gTAATGATGCCAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCAt                        <  1:304095/81‑1 (MQ=255)                                                       tgccAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCggt                 <  1:304096/78‑1 (MQ=255)                                                        gccAACAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCggtg                >  2:304123/3‑81 (MQ=255)                                                            aCAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGtggtgg            >  2:304120/1‑81 (MQ=255)                                                             cAGCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGtggtgg            >  2:304127/1‑80 (MQ=255)                                                               gCACCACGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCAt         <  1:304097/81‑1 (MQ=255)                                                                caccacGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCAt         >  2:304124/1‑80 (MQ=255)                                                                   cacGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACt     <  1:304098/81‑1 (MQ=255)                                                                    acGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACtt    <  1:304099/81‑1 (MQ=255)                                                                     cGGATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTa   <  1:304100/81‑1 (MQ=255)                                                                      ggATATTG‑AAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTAt  <  2:304101/81‑1 (MQ=255)                                                                              |                                                                          CTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTAT  >  NC_000913/1658650‑1658800

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