Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,649,690	A→T	A252A (GCA→GCT)	arsB →	arsenite/antimonite:H(+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,649,690	0	A	T	100.0%	120.2 / NA	33	A252A (GCA→GCT)	arsB	arsenite/antimonite:H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (16/17);  total (16/17)

GCTGGGTTGTTTTACTGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCAGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGTAAAGTCCTGCGCGGTGC  >  NC_000913/3649614‑3649767                                                                             |                                                                              gctggGTTGTTTTACTGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGctgctg                                                                           <  2:704838/81‑1 (MQ=255)    gggTTGTTTTACTGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTgg                                                                        >  2:704871/1‑81 (MQ=255)    gggTTGTTTTACTGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTgg                                                                        >  1:704873/1‑81 (MQ=255)       ttgttTTACTGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGgcg                                                                     >  1:704886/1‑81 (MQ=255)         gttTTACTGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGgcgcg                                                                   >  2:704888/1‑81 (MQ=255)         gttTTACTGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGgcgcg                                                                   >  1:704887/1‑81 (MQ=255)                tGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGatat                                                            >  1:704883/1‑81 (MQ=255)                  cTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGatat                                                            <  1:704839/79‑1 (MQ=255)                     ctgctgGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATattatt                                                        <  1:704840/80‑1 (MQ=255)                        ctgGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGAGCGCGCTGATATTATTTgtc                                                    <  2:704841/81‑1 (MQ=255)                           gtggGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTgtcgtc                                                 >  1:704885/1‑81 (MQ=255)                           gtggGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTgtcgtc                                                 <  2:704843/81‑1 (MQ=255)                           gtggGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTgtcgtc                                                 >  2:704878/1‑81 (MQ=255)                           gtggGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTgtcgtc                                                 <  1:704844/81‑1 (MQ=255)                           gtggGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTgtcgt                                                  <  2:704842/80‑1 (MQ=255)                                  ttttCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAAc                                          >  1:704881/1‑81 (MQ=255)                                   tttCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAAc                                          <  1:704845/80‑1 (MQ=255)                                       gTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGt                                     >  1:704896/1‑81 (MQ=255)                                         ccTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCa                                   >  1:704880/1‑81 (MQ=255)                                         ccTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCa                                   >  2:704895/1‑81 (MQ=255)                                          cTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCAt                                  <  2:704846/81‑1 (MQ=255)                                             gAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCg                               <  2:704847/81‑1 (MQ=255)                                             gAACCGCTCGGCATTCCGGTGAGCGCCATTGCTGATGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCg                               <  2:704848/81‑1 (MQ=255)                                               aCCGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGAt                             >  2:704899/1‑81 (MQ=255)                                                ccGCTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGAt                             >  1:704897/1‑80 (MQ=255)                                                   cTCGGCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAAt                         <  1:704849/81‑1 (MQ=255)                                                      ggCATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACg                      <  2:704850/81‑1 (MQ=255)                                                        cATTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACggg                    <  1:704851/81‑1 (MQ=255)                                                         aTTCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGt                   <  1:704852/81‑1 (MQ=255)                                                           tCCGGTGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGTa                  >  2:704867/1‑80 (MQ=255)                                                                tGAGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGTAAAGTcc            >  1:704898/1‑81 (MQ=255)                                                                  aGCGCCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGTAAAGTCCt           <  2:704853/80‑1 (MQ=255)                                                                     gcCATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGTAAAGTCCTgcgc       <  2:704854/81‑1 (MQ=255)                                                                      ccATTGCTGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGTAAAGTCCTGcgcg      <  1:704855/81‑1 (MQ=255)                                                                          tgctgcTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGTAAAGTCCTGCGCGGTGc  >  1:704891/1‑81 (MQ=255)                                                                             |                                                                              GCTGGGTTGTTTTACTGCTTCTGCTGGTGGGATTTTTCGTCCTGGAACCGCTCGGCATTCCGGTGAGCGCCATTGCAGCTGTGGGCGCGCTGATATTATTTGTCGTCGCTAAACGCGGTCATGCGATTAATACGGGTAAAGTCCTGCGCGGTGC  >  NC_000913/3649614‑3649767

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