Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,278	A→G	7.9%	I6I (ATT→ATC)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,278	0	A	G	7.9%	249.6 / 11.1	88	I6I (ATT→ATC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand):  ref base A (41/40);  new base G (5/2);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCG  >  minE/2595210‑2595348                                                                     |                                                                       ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa                                                                       >  1:252935/1‑70 (MQ=255)  gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                                     <  1:990271/71‑1 (MQ=255)  gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                                     <  1:744594/71‑1 (MQ=255)  gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                                     <  1:1268239/71‑1 (MQ=255)    tCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                                                    >  1:115969/1‑70 (MQ=255)    tCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                                                    >  1:3776337/1‑70 (MQ=255)     cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                                                   <  1:1919791/70‑1 (MQ=255)     cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                                                   <  1:1623907/70‑1 (MQ=255)        tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                                               >  1:1913872/1‑71 (MQ=255)        tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                                               >  1:1754070/1‑71 (MQ=255)        tgCCGCTGTTGCAGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                                                 <  1:648216/69‑1 (MQ=255)          ccGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGAAATGtt                                                             <  1:533673/71‑1 (MQ=255)           cGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTg                                                            >  1:1906664/1‑71 (MQ=255)               gTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCa                                                        >  1:2620261/1‑71 (MQ=255)               gTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCa                                                        >  1:3673576/1‑71 (MQ=255)                ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGt                                                              <  1:321464/64‑1 (MQ=255)                  gCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGTCGGGTTGTGATAAATGGTAATGTTGCTCatat                                                     >  1:3007022/1‑71 (MQ=255)                  gCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGtt                                                             >  1:2215643/1‑63 (MQ=255)                  gCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCatat                                                     >  1:730369/1‑71 (MQ=255)                  gCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCatat                                                     >  1:716957/1‑71 (MQ=255)                   cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCatat                                                     <  1:393707/70‑1 (MQ=255)                   cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATc                                                    <  1:1020104/71‑1 (MQ=255)                    ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                                                 >  1:769447/1‑57 (MQ=255)                     gATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                                                  >  1:1753051/1‑71 (MQ=255)                      atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCa                                                        <  1:2435469/64‑1 (MQ=255)                      agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAg                                                  >  1:3643465/3‑70 (MQ=255)                       tcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt                                                          <  1:2984954/61‑1 (MQ=255)                       tcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa                                                <  1:3167255/71‑1 (MQ=255)                         atcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAAt                                                                      <  1:3366685/47‑1 (MQ=255)                         atcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                                                   <  1:1638371/50‑1 (MQ=255)                          tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                                                        <  1:1790447/44‑1 (MQ=255)                           ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                                                   >  1:3276253/1‑48 (MQ=255)                           ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCa                                                   >  1:3938327/1‑64 (MQ=255)                           ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                                                  >  1:1224926/1‑65 (MQ=255)                           ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                                                  >  1:1876939/1‑65 (MQ=255)                           ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtct                                             >  1:3640407/1‑70 (MQ=255)                           ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtct                                             >  1:3268418/1‑70 (MQ=255)                           ctcCAGCGTATTACGCGACGTGCCGAAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                                                  >  1:2807732/1‑65 (MQ=255)                            tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCatat                                                     >  1:2004395/1‑61 (MQ=255)                            tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                                           >  1:3187360/1‑71 (MQ=255)                              cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtct                                             <  1:3338274/67‑1 (MQ=255)                              cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt                                          >  1:213189/1‑70 (MQ=255)                               aGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGtt                                                             >  1:1513749/1‑50 (MQ=255)                                 cGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACa                                      <  1:2099097/71‑1 (MQ=255)                                  gTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa                                      >  1:1597207/1‑70 (MQ=255)                                   tATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa                                      >  1:2672352/1‑69 (MQ=255)                                    aTTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                                                 <  1:38501/41‑1 (MQ=255)                                    aTTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                                                 <  1:546729/41‑1 (MQ=255)                                       aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa                                     >  1:2202706/1‑66 (MQ=255)                                       aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg                                 >  1:1639357/1‑70 (MQ=255)                                       aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg                                 <  1:2030771/70‑1 (MQ=255)                                        cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa                               <  1:2880073/71‑1 (MQ=255)                                        cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa                               <  1:2384280/71‑1 (MQ=255)                                        cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                                <  1:2113878/70‑1 (MQ=255)                                        cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                                <  1:2634848/70‑1 (MQ=255)                                         gcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                              <  1:2336324/71‑1 (MQ=255)                                          cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                                >  1:745275/1‑68 (MQ=255)                                           gACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                            <  1:77634/71‑1 (MQ=255)                                           gACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                            <  1:885400/71‑1 (MQ=255)                                            aggTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtc                                              >  1:2156008/3‑52 (MQ=255)                                            aCGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagag                           <  1:3456575/71‑1 (MQ=255)                                            aCGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagag                           <  1:3804033/71‑1 (MQ=255)                                              gTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAg                                                  >  1:2692657/1‑46 (MQ=255)                                               tGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                            <  1:427661/67‑1 (MQ=255)                                                 ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata                                                      >  1:4053232/1‑39 (MQ=255)                                                 ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga                          >  1:3346957/1‑67 (MQ=255)                                                 ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga                          <  1:3826643/67‑1 (MQ=255)                                                 ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACg                      <  1:2828871/71‑1 (MQ=255)                                                 ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGATCata                                                      >  1:2431832/1‑39 (MQ=255)                                                  cGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag                             >  1:2236511/1‑63 (MQ=255)                                                  cGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag                             >  1:1790879/1‑63 (MQ=255)                                                   acAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtct                                             <  1:4081679/45‑1 (MQ=255)                                                    cAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACagagagc                        >  1:1754559/1‑65 (MQ=255)                                                    cAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt                                                 <  1:3542693/41‑1 (MQ=255)                                                    cAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                              >  1:1723584/1‑60 (MQ=255)                                                    cAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                              >  1:2392040/1‑60 (MQ=255)                                                        ccGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGcc               <  1:1336311/71‑1 (MQ=255)                                                         cGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAt                                               <  1:1970770/38‑1 (MQ=255)                                                         cGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAgcgc                >  1:3884315/1‑69 (MQ=255)                                                          gggTTGTGATAAATGGTAATGTTGCTCATATCAGt                                                 <  1:2982313/35‑1 (MQ=255)                                                              tgtgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag                             >  1:1906276/1‑51 (MQ=255)                                                              tgtgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAgagaga                        <  1:1168801/56‑1 (MQ=255)                                                                tgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa                                      <  1:3507236/40‑1 (MQ=255)                                                                tgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGcc               <  1:2358810/63‑1 (MQ=255)                                                                tgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcag        <  1:1479919/70‑1 (MQ=255)                                                                 gATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa                                     >  1:829721/1‑40 (MQ=255)                                                                  aTAGATGGTAATGTTGCTCATATCAGTATCTAATTACaaa                                    >  1:3230023/1‑40 (MQ=38)                                                                  aTAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCa              >  1:3451772/1‑62 (MQ=255)                                                                  aTAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagc       >  1:31842/1‑69 (MQ=255)                                                                   tAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                                <  1:135342/43‑1 (MQ=255)                                                                     aaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCg  >  1:2351324/1‑71 (MQ=255)                                                                     aaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCg  <  1:3550584/71‑1 (MQ=255)                                                                     aaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCg  >  1:1405204/1‑71 (MQ=255)                                                                     |                                                                       GGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCG  >  minE/2595210‑2595348

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