Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,650,642	G→A	A136T (GCG→ACG)	arsC →	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,650,642	0	G	A	100.0%	40.1 / NA	16	A136T (GCG→ACG)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/7);  total (9/7)

TGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCACCCGATTCTGATTAATCGCCCGATTGTGGTGACGCCGCTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTAGGCCTGACAGGCGTAGCGCATCAGGCGATTTTGCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGC  >  NC_000913/3650412‑3650883                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  tgcGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCACCCGATTCTGATTAATCGCCCGATTGTGGTGACGCCGCTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAACggg                                                                                                                                                                                                                                               >  1:188974/1‑235 (MQ=255)                                                                               tGCTTCAGCACCCGATTCTGATTAATCGCCCGATTGTGGAGCCGCCGCTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGAGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGAGGAAACGGGTCAGCGCCTGAAATAACGCGACGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTGTcc                                                                                                                                                 >  1:245957/1‑251 (MQ=255)                                                                                 cTTCAGCACCCGATTCTGATTAATCGCCCGATTGTGGTGACGCCGCTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATACTGCCAGATGCGCAAAAAGGCGCATTATCCCAGGAAGATGGCGAGAAAGTGGTTGATGAAACGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCg                                                                                                                                                >  1:285573/1‑250 (MQ=255)                                                                                              ttCTGATTAATCGCCCGATTGTGGTGACGCCGCTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAACGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTggn                                                                                                                                                            <  1:150735/225‑2 (MQ=255)                                                                                                                        aCGCCGCTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGTAAGATGGCGAGAAAGGTGTTGATGAAACGGGGAAGCGCCTGAAATAAAGCGTCGATATCCCCCCCACGTGTTGTTAGAAAAGTGCGCTGTATTTATGTCGGATTCGACGCTGGTGCGGCTTATCCGAGCTACAAAGTTATGAAAATTCGATGAATAGTATGTTAtt                                                                                                        <  1:474721/251‑1 (MQ=255)                                                                                                                             gcTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGAGCATTCTCAAAGGAAGATGGCGAGAAAGTGGTTGATGAAACGGGTCAGCGCATGAAATAAAGCGGCGat                                                                                                                                                                                                                      >  1:344630/1‑136 (MQ=255)                                                                                                                             gcTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATGCNGCCNGATGNGCNAAAAGGCGCATTCTCNAAGGAAGATGGCGAGAAAGTGGTTGATGAAACGGGTAAGCGCCTNNANNNANGCGGCGat                                                                                                                                                                                                                      <  2:344630/136‑1 (MQ=255)                                                                                                                                               ccgccCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCGCCAAGGAAGATGGCGAGAAAGTGGTGGATGAAACGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGgcg                                                                                                                                                           <  1:498166/177‑1 (MQ=255)                                                                                                                                                              ggtgCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGAGAGAAAGTGGTTGATGAAACGAGTAAGCGCCTGAAATAAAGCGGCAATATCCCCCCCACAGGTTGTTGGAAAAGTGCGCTTTATTTATGTCAGATGCGACGCTGGCGCGACTTATCCGACGTACAAAGTTATGAAAATTCGATGACTAGAATGTTATTTGTAGGCCTGACAGTCGTAGAGCATAAGGCGATTTTGc                                                                  >  1:362041/1‑251 (MQ=255)                                                                                                                                                                                      gcAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAACGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGGTTTATTTATGTCGGATGCGACGCTGGCGCGGCTTATCCGACCGACAAAGTTATGAAAAATCGATGAATAGTATGTTATTTGTAGGCCTGACAGGCGTAGCGCATCAGGc                                                                          <  1:468756/219‑1 (MQ=255)                                                                                                                                                                                               ngcATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAACGGGTAAGCGCCTGAAATAACGCGGCGATATCCTCCCCACAGGTTGTTAGAAAAGTGCGCtttat                                                                                                                                                                                  >  1:93139/2‑106 (MQ=255)                                                                                                                                                                                               cgcATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAACGGGTAAGCGCCTGAAATACAGCAGCGCTATCCCCCCCACAGGTTGTTAGAAAAGCGCGCTTTATTTATGTCGGATGCGACGCTGGCGAGTCTTATACGACCTACAGAGTTATGAAAATTCGATGAATAGTATGTTCTTTGTAGGCCTGAAAGGCGTAGCGCATCAGGCGATTTTGCTATTTACACAGTACTCCCTGCGTGAGATTCCa                                  >  1:295077/1‑250 (MQ=255)                                                                                                                                                                                                                  ggCGANAAAGTNGNNGATGAAACGGGTAAGCGCCNGAAATAAAGCGGCGANATCCNCCCcan                                                                                                                                                                                                           >  2:202947/1‑61 (MQ=255)                                                                                                                                                                                                                  ggCGAGAAAGTGGTTGATGAAACGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTAGGCCTGACAGGCGTAGCGCATCAGGCGATTTTGCTATTTaca                                                          <  1:202947/207‑1 (MQ=255)                                                                                                                                                                                                                   ncGAGAAAGTGGTTGATGAAACGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACCAAGTTATGAAAATTCGATGAATAGTa                                                                                                               >  1:90689/2‑153 (MQ=255)                                                                                                                                                                                                                                tgatgaAACGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTAGGCCTGACAGGCGTAGCGCATCAGGCGATTTTGCTATTTACACAGTACTCCCTGCGTGATATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAgc  <  1:476601/249‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  TGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCACCCGATTCTGATTAATCGCCCGATTGTGGTGACGCCGCTGGGAACTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTAGGCCTGACAGGCGTAGCGCATCAGGCGATTTTGCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGC  >  NC_000913/3650412‑3650883

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