Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	636,528	T→C	E14E (GAA→GAG)	ybdN ←	putative PAPS reductase/DUF3440 domain‑containing protein YbdN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	636,528	0	T	C	100.0%	55.5 / NA	17	E14E (GAA→GAG)	ybdN	putative PAPS reductase/DUF3440 domain‑containing protein YbdN
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/11);  total (6/11)

CCGCCAGAAAATGAGACACAAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCTTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCTGAAGACGCTCTCTGAACGTTA  >  NC_000913/636460‑636592                                                                     |                                                                 ccgccAGAAAATGAGACACAAACGCGGGGTAGGGTGTTGAGTGTTCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTg                                                      <  2:138402/81‑1 (MQ=255)  cgccAGAAAATGAGACACAAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGa                                                     <  2:138403/81‑1 (MQ=255)   gccAGAAAATGAGACACAAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAg                                                    <  1:138404/81‑1 (MQ=255)    ccAGAAAATGAGACACAAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGc                                                   <  1:138405/81‑1 (MQ=255)      aGAAAATGAGACACAAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCgg                                                 <  2:138406/81‑1 (MQ=255)            tGAGACACAAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAgg                                           <  1:138407/81‑1 (MQ=255)              agaCACAAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGaa                                         <  2:138408/81‑1 (MQ=255)                 cacaAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAAttt                                      >  2:138428/1‑81 (MQ=255)                    aaaCGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTta                                    >  1:138429/1‑80 (MQ=255)                              tAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACAt                         <  1:138409/81‑1 (MQ=255)                                gggTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCt                       <  2:138410/81‑1 (MQ=255)                                gggTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCt                       >  2:138426/1‑81 (MQ=255)                                    gtTGAGTGTCCAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCTg                      >  2:138416/1‑78 (MQ=255)                                              cAGGTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCTGAAGACGctctct          <  1:138411/80‑1 (MQ=255)                                                 gTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCTGAAGACGCTCTCTGAAc      <  1:138412/81‑1 (MQ=255)                                                 gTAATACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCTGAAGACGCTCTCTGAAc      >  2:138431/1‑81 (MQ=255)                                                     tACGCTCTCGTGCAGCCTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCTGAAGACGCTCTCTGAACGTTa  >  2:138430/1‑81 (MQ=255)                                                                     |                                                                 CCGCCAGAAAATGAGACACAAACGCGGGGTAGGGTGTTGAGTGTCCAGGTAATACGCTCTCGTGCAGCTTCCAGAATATTGAGCGGAAGAGGAATTTTATAAATAGACATCTGAAGACGCTCTCTGAACGTTA  >  NC_000913/636460‑636592

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