Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,495,571	G→A	S95N (AGC→AAC)	ydcJ →	DUF1338 domain‑containing protein YdcJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,495,571	0	G	A	100.0%	64.6 / NA	19	S95N (AGC→AAC)	ydcJ	DUF1338 domain‑containing protein YdcJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/13);  total (6/13)

GCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGATGATGCTTC  >  NC_000913/1495495‑1495640                                                                             |                                                                      gCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTAt                                                                   <  1:253949/81‑1 (MQ=255) gCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTAt                                                                   <  2:253948/81‑1 (MQ=255)       cgcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGAt                                                             >  1:253977/1‑81 (MQ=255)          gTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGAtctc                                                          >  1:253973/1‑81 (MQ=255)           ttGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGAtctct                                                         >  2:253971/1‑81 (MQ=255)              ggACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCg                                                       <  2:253950/80‑1 (MQ=255)                          agCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGc                                          >  2:253981/1‑81 (MQ=255)                          agCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGc                                          <  2:253951/81‑1 (MQ=255)                               gCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGt                                      <  2:253952/80‑1 (MQ=255)                                cTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTAc                                    <  2:253953/81‑1 (MQ=255)                                         ggcggATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACa                            <  1:253954/80‑1 (MQ=255)                                         ggcggATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAg                           <  1:253955/81‑1 (MQ=255)                                          gcggATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGc                          <  1:253956/81‑1 (MQ=255)                                           cggATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCa                         >  1:253987/1‑81 (MQ=255)                                            ggATGTTTGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAt                        >  2:253975/1‑81 (MQ=255)                                                 tttGCCATTATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCg                    <  2:253957/80‑1 (MQ=255)                                                        ttATGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtg            <  2:253958/81‑1 (MQ=255)                                                           tGGGGATGTACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtgatg         <  1:253959/81‑1 (MQ=255)                                                                   tACCCGGTGAACTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGATGATGCTtc  <  1:253960/80‑1 (MQ=255)                                                                             |                                                                      GCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGATGATGCTTC  >  NC_000913/1495495‑1495640

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