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%	69.1 / NA	20	intergenic (+459/‑170)	arsC/yhiS	arsenate reductase/putative uncharacterized protein YhiS
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (8/12);  total (8/12)

CAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCC‑TTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAT  >  NC_000913/3651044‑3651191                                                                               |                                                                       cAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCC‑ttt                                                                     >  1:709293/1‑81 (MQ=255)    agaACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCC‑tttttt                                                                  >  2:709297/1‑81 (MQ=255)        cATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCC‑tttttttta                                                               >  2:709305/1‑78 (MQ=255)                   tttAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAatat                                                     <  2:709272/80‑1 (MQ=255)                      aaaaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTAcc                                                 >  2:709295/1‑81 (MQ=255)                       aaaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCg                                                >  1:709301/1‑81 (MQ=255)                       aaaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCg                                                <  2:709273/81‑1 (MQ=255)                        aaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCgg                                               >  2:709299/1‑81 (MQ=255)                          tGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGt                                             <  2:709274/81‑1 (MQ=255)                          tGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGt                                             <  1:709275/81‑1 (MQ=255)                           gAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTg                                            >  1:709292/1‑81 (MQ=255)                            aGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTgg                                           <  1:709276/81‑1 (MQ=255)                                cATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGgaga                                        >  1:709303/1‑80 (MQ=255)                                     tGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATaaaa                                   <  1:709277/80‑1 (MQ=255)                                                  cGATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGt                     <  1:709278/81‑1 (MQ=255)                                                   gATTTCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTa                    <  2:709279/81‑1 (MQ=255)                                                     tttCCTTTACGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAAt                  <  2:709280/81‑1 (MQ=255)                                                              cGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAgg         <  1:709281/81‑1 (MQ=255)                                                              cGATACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAg          >  1:709298/1‑80 (MQ=255)                                                                aTACTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGaa       <  2:709282/81‑1 (MQ=255)                                                                   cTATCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGa     >  2:709300/1‑80 (MQ=255)                                                                     aTCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGata   >  1:709304/1‑80 (MQ=255)                                                                      tCACTACCCTTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGatat  <  2:709283/80‑1 (MQ=255)                                                                               |                                                                       CAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCC‑TTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAT  >  NC_000913/3651044‑3651191

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