Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,230	G→A	12.3%	I22I (ATC→ATT) ‡	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,230	0	G	A	12.3%	237.5 / 19.0	89	I22I (ATC→ATT) ‡	arsC	arsenate reductase
Reads supporting (aligned to +/- strand):  ref base G (48/29);  new base A (5/6);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAG  >  minE/2595161‑2595300                                                                      |                                                                       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: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:879926/71‑1 (MQ=255)                         aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                               >  1:3305580/1‑71 (MQ=255)                         aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                               >  1:3149622/1‑71 (MQ=255)                           tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                             >  1:2440572/1‑71 (MQ=255)                           tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                             >  1:2816590/1‑71 (MQ=255)                           tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                             >  1:3519560/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:3419162/1‑57 (MQ=255)                                aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                      >  1:3889924/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: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: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:1867950/71‑1 (MQ=255)                                         ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                               <  1:817088/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:3801900/71‑1 (MQ=255)                                            ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                            <  1:57239/71‑1 (MQ=255)                                            ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                            <  1:676431/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)                                                         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:1855616/1‑67 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt          >  1:2881427/1‑67 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata      >  1:2731512/1‑71 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata      >  1:497094/1‑71 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata      >  1:2579112/1‑71 (MQ=255)                                                                  tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata      >  1:1013104/1‑71 (MQ=255)                                                                    cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                        <  1:1213714/51‑1 (MQ=255)                                                                     ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                        <  1:3132729/50‑1 (MQ=255)                                                                     ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                 <  1:2820987/57‑1 (MQ=255)                                                                     ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCat       <  1:3618691/67‑1 (MQ=255)                                                                      gATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg  <  1:2906178/71‑1 (MQ=255)                                                                      |                                                                       GACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAG  >  minE/2595161‑2595300

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