Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,215	T→G	12.9%	T27T (ACA→ACC) ‡	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,215	0	T	G	12.9%	182.2 / 26.9	85	T27T (ACA→ACC) ‡	arsC	arsenate reductase
Reads supporting (aligned to +/- strand):  ref base T (45/29);  new base G (5/6);  total (50/35)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGT  >  minE/2595147‑2595283                                                                     |                                                                     tCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCtgt                                                                    <  1:785359/71‑1 (MQ=255) tCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCtg                                                                     <  1:1097938/70‑1 (MQ=255) tCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCt                                                                      <  1:625538/69‑1 (MQ=255)   ggCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCtgtg                                                                   <  1:2749632/70‑1 (MQ=255)      aaTGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGc                                                               >  1:3033206/1‑71 (MQ=255)      aaTGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGc                                                               >  1:2136200/1‑71 (MQ=255)       aTGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGc                                                                  >  1:1616424/1‑67 (MQ=255)       aTGAGTTTGACAAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCt                                                              <  1:3505962/71‑1 (MQ=255)        tGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTg                                                             >  1:2862579/1‑71 (MQ=255)          aGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGtt                                                           <  1:928467/71‑1 (MQ=255)          aGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGt                                                            <  1:1997456/70‑1 (MQ=255)          aGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGt                                                            <  1:1969198/70‑1 (MQ=255)               gACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGa                                                      >  1:453227/1‑71 (MQ=255)                 ccAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGat                                                     >  1:1443439/1‑70 (MQ=255)                 ccAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGa                                                      >  1:3577499/1‑69 (MQ=255)                     ttCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtct                                                >  1:2661998/1‑71 (MQ=255)                       cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCtg                                                                     >  1:1700646/1‑48 (MQ=255)                           gcgcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCt                                                                      >  1:2429815/1‑43 (MQ=255)                           gcgcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCt                                                                      >  1:2988032/1‑43 (MQ=255)                           gcgcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                          >  1:1073230/1‑71 (MQ=255)                                 tGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAc                                    >  1:3429468/1‑71 (MQ=255)                                  ggCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCt                                                              >  1:249500/1‑44 (MQ=255)                                   gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCa                                             >  1:1874774/1‑60 (MQ=255)                                   gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                  <  1:3148513/71‑1 (MQ=255)                                   gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                  <  1:2682631/71‑1 (MQ=255)                                   gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGATGTTGCGGATCATCTCCAGCGTATTAcgc                                  <  1:572767/71‑1 (MQ=255)                                    cGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgcg                                 <  1:3269313/71‑1 (MQ=255)                                       aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                     >  1:833701/1‑64 (MQ=255)                                       aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                              >  1:3149622/1‑71 (MQ=255)                                       aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                              <  1:879926/71‑1 (MQ=255)                                       aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                              >  1:3305580/1‑71 (MQ=255)                                         tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                            >  1:2816590/1‑71 (MQ=255)                                         tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                            >  1:3519560/1‑71 (MQ=255)                                         tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                            >  1:2440572/1‑71 (MQ=255)                                          ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                           >  1:1903062/1‑71 (MQ=255)                                          ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                           >  1:1793698/1‑71 (MQ=255)                                           tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                           <  1:210215/67‑1 (MQ=255)                                           tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                           >  1:2230615/4‑70 (MQ=255)                                           tCCAGATAATGGATAATAGTCGGTTCTGTGCCGTTGTTGCGGATCAtctc                                               <  1:124172/50‑1 (MQ=255)                                           tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                           <  1:1160763/70‑1 (MQ=255)                                           tCCAGATAATGGATAATAGTCGGTTCTGTGACGCTGTTGCGGATCAtctc                                               <  1:2477670/50‑1 (MQ=255)                                            ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtctc                                               >  1:3559069/1‑49 (MQ=255)                                            ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                            >  1:2373334/1‑52 (MQ=255)                                            ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                          <  1:3564941/70‑1 (MQ=255)                                            ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                          <  1:1578844/70‑1 (MQ=255)                                             gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGccac                        <  1:1365175/70‑3 (MQ=255)                                             cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGt                                                            >  1:1840798/1‑35 (MQ=255)                                             cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                        <  1:3867493/71‑1 (MQ=255)                                             cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                        <  1:545457/71‑1 (MQ=255)                                              aGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCcaca                       <  1:3064865/71‑1 (MQ=255)                                              aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                     >  1:3419162/1‑57 (MQ=255)                                              aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                     >  1:3889924/1‑57 (MQ=255)                                              aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                         >  1:2546880/1‑69 (MQ=255)                                              aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                         >  1:325666/1‑69 (MQ=255)                                                aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAgg                     >  1:31737/1‑71 (MQ=255)                                                 tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                    >  1:2706900/1‑71 (MQ=255)                                                 tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                    <  1:2711013/71‑1 (MQ=255)                                                  aaTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCa                       >  1:2572014/1‑67 (MQ=255)                                                   aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGGGACGTGCCGCAGGcc                   <  1:1000847/70‑1 (MQ=255)                                                   aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                         <  1:3754353/64‑1 (MQ=255)                                                    tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAt                                       <  1:3019104/49‑1 (MQ=255)                                                    tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCgg                 >  1:693838/1‑71 (MQ=255)                                                    tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCgg                 >  1:696829/1‑71 (MQ=255)                                                     ggATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAt                                       >  1:1299870/1‑48 (MQ=255)                                                       ataataGTCGGTTCTGTGCCGCTGTTGCTGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt              <  1:1677810/71‑1 (MQ=255)                                                       ataataGTCGGTTCTGTGCCGCTGTTGCGGatcatc                                                 <  1:2705328/36‑1 (MQ=255)                                                       ataataGTCGGTTCTGTGCCGCTGTTGCGGatcat                                                  >  1:1218715/1‑35 (MQ=255)                                                       ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt              <  1:1060286/71‑1 (MQ=255)                                                       ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt              <  1:817088/71‑1 (MQ=255)                                                       ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt              <  1:1867950/71‑1 (MQ=255)                                                       ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGt               >  1:3857410/1‑70 (MQ=255)                                                        taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                >  1:2391373/1‑52 (MQ=255)                                                        taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                >  1:3084155/1‑52 (MQ=255)                                                         aataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt            >  1:1556364/1‑71 (MQ=255)                                                          attgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt            <  1:2587251/67‑1 (MQ=255)                                                          ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg           <  1:676431/71‑1 (MQ=255)                                                          ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg           <  1:57239/71‑1 (MQ=255)                                                          ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg           <  1:3801900/71‑1 (MQ=255)                                                           ttgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGa          <  1:672465/69‑1 (MQ=255)                                                           taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                     >  1:1489284/1‑44 (MQ=255)                                                           taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                     >  1:3483684/1‑44 (MQ=255)                                                           taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                        >  1:2458016/1‑57 (MQ=255)                                                           taGTCGGTTCTGTGCAGCTGTTGCGGATCATCTCCAGCGTATTAcg                                   >  1:1114239/1‑46 (MQ=255)                                                            tgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgtg           >  1:3721245/2‑69 (MQ=255)                                                             gTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgtg           >  1:217582/1‑68 (MQ=255)                                                              tCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATa        <  1:1639213/70‑1 (MQ=255)                                                               cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa       >  1:1503227/1‑70 (MQ=255)                                                                ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa      <  1:3376541/70‑1 (MQ=255)                                                                ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAAt     >  1:166687/1‑71 (MQ=255)                                                                ggTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGa                                >  1:1047466/1‑44 (MQ=255)                                                                ggTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGAt     >  1:3355037/1‑71 (MQ=255)                                                                    cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt  >  1:322149/1‑70 (MQ=255)                                                                     tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCggg                <  1:1398041/55‑1 (MQ=255)                                                                     |                                                                     TCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGT  >  minE/2595147‑2595283

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