Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,596,103	A→C	9.0%	A165A (GCT→GCG)	arsB ←	arsenite/antimonite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,596,103	0	A	C	9.0%	192.3 / 13.8	78	A165A (GCT→GCG)	arsB	arsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (31/40);  new base C (3/4);  total (34/44)
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.78e-01

ACAATCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAT  >  minE/2596035‑2596166                                                                      |                                                                acaaTCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGc                                                               >  1:781071/1‑71 (MQ=255)  caaTCTCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAg                                                                <  1:524461/69‑1 (MQ=255)   aaTCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGc                                                               >  1:3472154/1‑69 (MQ=255)     tCGCGGCGATATCCACCGGCACCATTACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGa                                                            >  1:3344881/1‑70 (MQ=255)      cgcgGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                          >  1:3327534/1‑71 (MQ=255)      cgcgGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                          >  1:762590/1‑71 (MQ=255)       gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                         <  1:3357981/71‑1 (MQ=255)       gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                         <  1:3281810/71‑1 (MQ=255)            gATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGatat                                                    <  1:1127373/71‑1 (MQ=255)            gATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGatat                                                    <  1:212430/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:496861/50‑1 (MQ=255)                                          ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGt                                           <  1:2267484/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:635676/1‑70 (MQ=255)                                                  aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                >  1:2715622/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:2669200/44‑1 (MQ=255)                                                             aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATaa                               <  1:2609176/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)                                                                      |                                                                ACAATCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAT  >  minE/2596035‑2596166

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