Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,660,717	G→A	A54A (GCG→GCA)	ynfF →	putative selenate reductase YnfF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,660,717	0	G	A	100.0%	100.2 / NA	28	A54A (GCG→GCA)	ynfF	putative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/12);  total (16/12)

TGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACT  >  NC_000913/1660639‑1660793                                                                               |                                                                             tGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGt                                                                            >  1:308989/1‑81 (MQ=255)   cgcTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCt                                                                          <  1:308990/81‑1 (MQ=255)     cTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTg                                                                         >  1:309020/1‑80 (MQ=255)        gCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGt                                                                     >  2:309016/1‑81 (MQ=255)          aaGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTc                                                                   >  1:309027/1‑81 (MQ=255)          aaGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTc                                                                   >  1:309017/1‑81 (MQ=255)            gCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCt                                                                 <  2:308991/81‑1 (MQ=255)              aGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGc                                                               <  2:843614/81‑1 (MQ=255)                         tctGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTg                                                    >  1:309015/1‑81 (MQ=255)                             ccATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCggg                                                <  2:308992/81‑1 (MQ=255)                               aTTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAg                                              >  2:309021/1‑81 (MQ=255)                                     ccAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGct                                         <  2:308993/80‑1 (MQ=255)                                         aTGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGctgtctg                                    <  1:308994/81‑1 (MQ=255)                                            gTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGtt                                  <  1:308995/80‑1 (MQ=255)                                               cGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGt                              >  1:309034/1‑81 (MQ=255)                                                tggcggcAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGGGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGtt                             >  2:309033/2‑81 (MQ=255)                                                    ggcAGAGGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCa                         >  1:309025/1‑81 (MQ=255)                                                       agagGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCAtgt                      >  2:309024/1‑81 (MQ=255)                                                       agagGCTCCGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCAtgt                      <  1:308996/81‑1 (MQ=255)                                                               cGGTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGAt               >  1:309023/1‑80 (MQ=255)                                                                 gTAGAAGAGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGac            <  1:308997/81‑1 (MQ=255)                                                                    gaagaGAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACAcc         >  2:309018/1‑81 (MQ=255)                                                                      agagAAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCgt       >  2:308998/1‑81 (MQ=255)                                                                       gagaAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCgtg      <  1:308999/81‑1 (MQ=255)                                                                       gagaAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCgtg      <  1:309000/81‑1 (MQ=255)                                                                        agaAAGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCgtgt     <  2:309001/81‑1 (MQ=255)                                                                          aaaGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTAc   >  1:309026/1‑81 (MQ=255)                                                                           aaGCAGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGCGTACt  >  2:309030/1‑81 (MQ=255)                                                                               |                                                                             TGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACT  >  NC_000913/1660639‑1660793

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