Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,596,106	G→C	9.2%	S164S (TCC→TCG)	arsB ←	arsenite/antimonite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,596,106	0	G	C	9.2%	219.2 / 9.8	76	S164S (TCC→TCG)	arsB	arsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (30/39);  new base C (3/4);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCA  >  minE/2596039‑2596176                                                                     |                                                                       tCGCGGCGATATCCACCGGCACCATTACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGa                                                                      >  1:3344881/1‑70 (MQ=255)  cgcgGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                                    >  1:762590/1‑71 (MQ=255)  cgcgGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                                    >  1:3327534/1‑71 (MQ=255)   gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                                   <  1:3357981/71‑1 (MQ=255)   gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                                   <  1:3281810/71‑1 (MQ=255)        gATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGatat                                                              <  1:212430/71‑1 (MQ=255)        gATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGatat                                                              <  1:1127373/71‑1 (MQ=255)         atatCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATAtt                                                             <  1:162433/71‑1 (MQ=255)          tatCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                                   <  1:2776925/64‑1 (MQ=255)          tatCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATAtt                                                             >  1:1732809/1‑70 (MQ=255)             ccACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                          >  1:3179704/1‑70 (MQ=255)               cccGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                          <  1:1863816/67‑1 (MQ=255)                ccGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAgg                                                      <  1:606113/71‑1 (MQ=255)                 cGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACg                                                                  >  1:2931477/1‑58 (MQ=255)                 cGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                               <  1:465757/61‑1 (MQ=255)                 cGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAg                                                       <  1:3276202/69‑1 (MQ=255)                   gCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                                   <  1:2632470/71‑1 (MQ=255)                     aCCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                                 <  1:2427126/71‑1 (MQ=255)                       cATAACGGAGGCGTACT‑GCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAgaaa                                               >  1:1756316/1‑69 (MQ=255)                        aTCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:855573/71‑1 (MQ=255)                         tCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGa                                                                 >  1:1560873/1‑51 (MQ=255)                         tCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              >  1:228468/1‑70 (MQ=255)                           aCCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                                 >  1:3771641/1‑65 (MQ=255)                           aCCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                                 >  1:2201873/1‑65 (MQ=255)                            ccGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                          <  1:1929019/71‑1 (MQ=255)                             cGAGGCGTATTCGCT‑AAAGCCGAGGCCAAAGAAATCAGCGGa                                                                      <  1:111261/42‑1 (MQ=255)                             cGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATaa                                         <  1:453881/71‑1 (MQ=255)                              gAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAg                                        >  1:1513322/1‑71 (MQ=255)                                    tATTCGCT‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGat                                                                <  1:644957/41‑1 (MQ=255)                                    tATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                          <  1:2234769/47‑1 (MQ=255)                                     aTTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAgg                                 >  1:1371377/1‑71 (MQ=255)                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGatat                                                              >  1:1864844/1‑41 (MQ=255)                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGa                                                                 <  1:1952124/38‑1 (MQ=255)                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATAtt                                                             <  1:1329411/42‑1 (MQ=255)                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                          <  1:3294984/45‑1 (MQ=255)                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGt                                                     <  1:2267484/50‑1 (MQ=255)                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGt                                                     <  1:496861/50‑1 (MQ=255)                                        cgcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCa                                                           >  1:633443/1‑42 (MQ=255)                                        cgcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCa                                                        >  1:3599048/1‑45 (MQ=255)                                        cgcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAgg                                                      <  1:389335/47‑1 (MQ=255)                                        cgcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAgg                                                      <  1:3072550/47‑1 (MQ=255)                                        cgcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTg                              >  1:274491/1‑71 (MQ=255)                                         gcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTgg                             >  1:1117304/1‑71 (MQ=255)                                         gcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTgg                             <  1:1896383/71‑1 (MQ=255)                                          cg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:2733741/53‑1 (MQ=255)                                          cg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:2396547/53‑1 (MQ=255)                                          cg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:2492156/53‑1 (MQ=255)                                           g‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATAtt                                                             <  1:876438/37‑1 (MQ=255)                                           g‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg                           <  1:2789780/71‑1 (MQ=255)                                             aaaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGc                                       >  1:744467/1‑58 (MQ=255)                                              aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                          >  1:2877939/1‑70 (MQ=255)                                              aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                          >  1:2715622/1‑70 (MQ=255)                                              aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                          >  1:635676/1‑70 (MQ=255)                                              aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGt                         <  1:203151/71‑1 (MQ=255)                                                gCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTAt                       <  1:2529413/71‑1 (MQ=255)                                                 ccGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATc                      >  1:3677093/1‑71 (MQ=255)                                                  cGAGGCCAAAGAAATCAGCGTAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCg                     <  1:3030136/71‑1 (MQ=255)                                                    aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTgg                                                  >  1:3900979/1‑40 (MQ=255)                                                    aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                   <  1:1731717/71‑1 (MQ=255)                                                      gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCaa                 >  1:2084086/1‑71 (MQ=255)                                                       ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                <  1:2061612/71‑1 (MQ=255)                                                       ccAAAGAAATAAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:1967548/41‑1 (MQ=255)                                                        cgaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTa                        <  1:2746529/60‑1 (MQ=255)                                                         gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACaa                                            >  1:239345/2‑41 (MQ=12)                                                         gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAAtg                                          <  1:881907/42‑2 (MQ=16)                                                         gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAtaa              <  1:1599426/70‑3 (MQ=255)                                                         aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATaa                                         <  1:2609176/44‑1 (MQ=255)                                                         aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATaa                                         <  1:2669200/44‑1 (MQ=255)                                                          aagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATaaa             <  1:2122001/71‑1 (MQ=255)                                                          aaGAAATCAGCGGAAACGATATTCACCAGGTTGgaga                                               >  1:1524877/1‑37 (MQ=255)                                                          aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa             >  1:3449102/1‑71 (MQ=255)                                                          aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa             >  1:333578/1‑71 (MQ=255)                                                           aGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATct                     >  1:675806/1‑61 (MQ=255)                                                           aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg                           >  1:3522283/1‑56 (MQ=255)                                                           aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt            >  1:3510425/1‑71 (MQ=255)                                                                    cGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgc    >  1:70191/1‑70 (MQ=255)                                                                     ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGTTATCGGCAATGAATCCGGCCGCCa  >  1:3568100/1‑71 (MQ=255)                                                                     |                                                                       TCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCA  >  minE/2596039‑2596176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑