Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,227	G→A	12.5%	R23R (CGC→CGT) ‡	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,227	0	G	A	12.5%	237.5 / 19.8	88	R23R (CGC→CGT) ‡	arsC	arsenate reductase
Reads supporting (aligned to +/- strand):  ref base G (48/29);  new base A (5/6);  total (53/35)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

GACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATA  >  minE/2595161‑2595296                                                                   |                                                                      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)             gcgcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                                       >  1:1073230/1‑71 (MQ=255)                   tGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAc                                                 >  1:3429468/1‑71 (MQ=255)                     gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCa                                                          >  1:1874774/1‑60 (MQ=255)                     gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                               <  1:2682631/71‑1 (MQ=255)                     gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                               <  1:3148513/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:3305580/1‑71 (MQ=255)                         aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                           >  1:3149622/1‑71 (MQ=255)                         aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                           <  1:879926/71‑1 (MQ=255)                           tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                         >  1:3519560/1‑71 (MQ=255)                           tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                         >  1:2816590/1‑71 (MQ=255)                           tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                         >  1:2440572/1‑71 (MQ=255)                            ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                                        >  1:1793698/1‑71 (MQ=255)                            ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                                        >  1:1903062/1‑71 (MQ=255)                             tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                                        >  1:2230615/4‑70 (MQ=255)                             tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                                        <  1:210215/67‑1 (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:1578844/70‑1 (MQ=255)                              ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                       <  1:3564941/70‑1 (MQ=255)                               gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGccac                                     <  1:1365175/70‑3 (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:3889924/1‑57 (MQ=255)                                aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                  >  1:3419162/1‑57 (MQ=255)                                aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                      >  1:325666/1‑69 (MQ=255)                                aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                      >  1:2546880/1‑69 (MQ=255)                                  aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAgg                                  >  1:31737/1‑71 (MQ=255)                                   tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                 <  1:2711013/71‑1 (MQ=255)                                   tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                 >  1:2706900/1‑71 (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:696829/1‑71 (MQ=255)                                      tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCgg                              >  1:693838/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:3084155/1‑52 (MQ=255)                                          taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                             >  1:2391373/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:3483684/1‑44 (MQ=255)                                             taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                  >  1:1489284/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)                                                         tgCCGCTGTTGCGTATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                        <  1:1139418/58‑1 (MQ=255)                                                         tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg           >  1:2599270/1‑71 (MQ=255)                                                           ccGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa             >  1:3530304/1‑67 (MQ=255)                                                               tGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt      >  1:363305/1‑70 (MQ=255)                                                                 ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                 <  1:29745/57‑1 (MQ=255)                                                                 ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTT‑CTCat   <  1:2393005/70‑1 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt      >  1:2881427/1‑67 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt      >  1:1855616/1‑67 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata  >  1:2731512/1‑71 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata  >  1:2579112/1‑71 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata  >  1:497094/1‑71 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata  >  1:1013104/1‑71 (MQ=255)                                                                   |                                                                      GACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATA  >  minE/2595161‑2595296

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