Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,651,128	(T)7→8	intergenic (+466/‑163)	arsC → / → yhiS	arsenate reductase/putative uncharacterized protein YhiS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,651,121	1	.	T	100.0%	63.2 / NA	19	intergenic (+459/‑170)	arsC/yhiS	arsenate reductase/putative uncharacterized protein YhiS
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (10/9);  total (10/9)

CTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCC‑TTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGA  >  NC_000913/3651059‑3651188                                                                |                                                                    cTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAat                                                    <  2:699283/81‑1 (MQ=255)        aaaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTAcc                                              <  2:699284/80‑1 (MQ=255)        aaaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCg                                             >  2:699307/1‑81 (MQ=255)               cgcATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGgag                                      <  1:699285/81‑1 (MQ=255)                gcATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGgaga                                     <  1:699286/81‑1 (MQ=255)                  aTTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGAt                                    >  2:699312/1‑80 (MQ=255)                     ttGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATaaa                                 <  2:699287/80‑1 (MQ=255)                           gttgttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGaaa                          >  1:699313/1‑81 (MQ=255)                           gttgttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGaa                           <  1:699288/80‑1 (MQ=255)                               ttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCa                       <  2:699289/80‑1 (MQ=255)                               ttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCa                       >  2:699306/1‑80 (MQ=255)                                tAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCa                       >  1:699317/1‑79 (MQ=255)                                tAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCa                       >  2:699310/1‑79 (MQ=255)                                 aaCGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAg                    >  2:699314/1‑81 (MQ=255)                                        tCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATta             >  2:699305/1‑81 (MQ=255)                                            ttACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTatat          >  2:699315/1‑80 (MQ=255)                                                 aTACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGaa    <  2:699290/81‑1 (MQ=255)                                                  tACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAg   <  1:699291/81‑1 (MQ=255)                                                    cTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGa  >  2:699308/1‑80 (MQ=255)                                                                |                                                                    CTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCC‑TTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGA  >  NC_000913/3651059‑3651188

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