Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,328,943	T→C	N57S (AAC→AGC)	yciK ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,328,943	0	T	C	100.0%	66.3 / NA	20	N57S (AAC→AGC)	yciK	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/14);  total (6/14)

CCAGTTGTTGGCAATTTTCGGACGTGCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGTTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATCACTGTCGCACCATAGCGTGCA  >  NC_000913/1328866‑1329017                                                                              |                                                                           ccAGTTGTTGGCAATTTTCGGACGTGCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTa                                                                         >  2:257321/1‑81 (MQ=255)    gttgttGGCAATTTTCGGACGTGCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTAtgt                                                                      <  1:257283/81‑1 (MQ=255)    gttgttGGCAATTTTCGGACGTGCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTAtgt                                                                      >  1:257326/1‑81 (MQ=255)          ggCAATTTTCGGACGTGCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGtggctg                                                                 <  2:257284/80‑1 (MQ=255)                    ggACGTGCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGAc                                                       <  2:257285/80‑1 (MQ=255)                     gACGTGCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGt                                                     <  1:257286/81‑1 (MQ=255)                          gCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAAttt                                                >  2:257319/1‑81 (MQ=255)                           cAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAAtttt                                               <  2:257287/81‑1 (MQ=255)                           cAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAAtttt                                               <  2:257288/81‑1 (MQ=255)                                    aaaTCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTa                                      >  2:257320/1‑81 (MQ=255)                                     aaTCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTAc                                     <  1:257289/81‑1 (MQ=255)                                         gagaATAAACCACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGcc                                 >  1:257316/1‑81 (MQ=255)                                                  ccACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAAt                        <  2:257290/81‑1 (MQ=255)                                                   cACTGTGGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATc                       <  2:257291/81‑1 (MQ=255)                                                        tgGCTGACGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATCACt                    <  1:257292/79‑1 (MQ=255)                                                              aCGCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATCACTGTCGCAc             <  1:257293/80‑1 (MQ=255)                                                                gCCCAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATCACTGTCGCACCa           >  1:257333/1‑80 (MQ=255)                                                                   cAGTTTCTTCGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATCACTGTCGCACCATAGc       <  2:257294/81‑1 (MQ=255)                                                                        tcttcGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATCACTGTCGCACCATAGCGTGc   <  2:257295/80‑1 (MQ=255)                                                                         cttcGCTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATCACTGTCGCACCATAGCGTGCa  <  2:257296/80‑1 (MQ=255)                                                                              |                                                                           CCAGTTGTTGGCAATTTTCGGACGTGCAGGTCAGCAAATCGAGAATAAACCACTGTGGCTGACGCCCAGTTTCTTCGTTTATGTGGCTGGCTACCTGACGTAATTTTTCTTCATTACGGCCCAACAGAATCACTGTCGCACCATAGCGTGCA  >  NC_000913/1328866‑1329017

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: ‑