Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,260	G→A	5.4%	C12C (TGC→TGT)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,260	0	G	A	5.4%	310.6 / 5.2	93	C12C (TGC→TGT)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand):  ref base G (42/46);  new base A (3/2);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACG  >  minE/2595195‑2595328                                                                  |                                                                     tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                                                 <  1:2994010/71‑1 (MQ=255) tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                                                 <  1:48403/71‑1 (MQ=255)  aaTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGcc                                                                >  1:1606428/1‑71 (MQ=255)        taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                           <  1:1621966/70‑1 (MQ=255)        taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtg                                                          <  1:2720924/71‑1 (MQ=255)         aataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                                         >  1:3549476/1‑71 (MQ=255)          attgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                                         <  1:3246237/67‑1 (MQ=255)          ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGtgttgt                                                         <  1:387680/70‑1 (MQ=255)             gTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATa                                                     <  1:2374799/71‑1 (MQ=255)                ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                                         >  1:905731/1‑64 (MQ=255)                ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa                                                   >  1:252935/1‑70 (MQ=255)                 gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                 <  1:744594/71‑1 (MQ=255)                 gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                 <  1:1268239/71‑1 (MQ=255)                 gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                 <  1:990271/71‑1 (MQ=255)                   tCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                                                     >  1:2423743/1‑49 (MQ=255)                   tCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                                >  1:3776337/1‑70 (MQ=255)                   tCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                                >  1:115969/1‑70 (MQ=255)                    cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                               <  1:1623907/70‑1 (MQ=255)                    cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                               <  1:1919791/70‑1 (MQ=255)                       tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                           >  1:1754070/1‑71 (MQ=255)                       tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                           >  1:1913872/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)                             tGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATa                                                     <  1:859934/55‑1 (MQ=255)                              gTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATa                                                     <  1:4040832/54‑1 (MQ=255)                              gTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCa                                    >  1:3673576/1‑71 (MQ=255)                              gTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCa                                    >  1:2620261/1‑71 (MQ=255)                               ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                           <  1:1070953/47‑1 (MQ=255)                               ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                                         <  1:2420007/49‑1 (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)                                    gATCATCTCCAGCGTATTACGCGACGTGCCGCAGGcc                                                                <  1:3254149/37‑1 (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)                                       catcTCCAGCGTATTACGCGACGTGCCGCAGGCCggg                                                             >  1:1778797/1‑37 (MQ=255)                                       catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGa                                                       >  1:4017217/1‑43 (MQ=255)                                        atcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATa                                                     <  1:453207/44‑1 (MQ=255)                                        atcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAAt                                                  <  1:3366685/47‑1 (MQ=255)                                        atcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                               <  1:1638371/50‑1 (MQ=255)                                         tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                                    <  1:1790447/44‑1 (MQ=255)                                         tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATa                                                     <  1:147999/43‑1 (MQ=255)                                          ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                               >  1:3276253/1‑48 (MQ=255)                                          ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCa                               >  1:3938327/1‑64 (MQ=255)                                          ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                              >  1:1876939/1‑65 (MQ=255)                                          ctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                              >  1:1224926/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)                                              aGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGa                                                       >  1:1117714/1‑36 (MQ=255)                                              aGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATa                                                     <  1:1287238/38‑1 (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:546729/41‑1 (MQ=255)                                                   aTTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                             <  1:38501/41‑1 (MQ=255)                                                      aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa                 >  1:2202706/1‑66 (MQ=255)                                                      aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg             <  1:2030771/70‑1 (MQ=255)                                                      aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg             >  1:1639357/1‑70 (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:885400/71‑1 (MQ=255)                                                          gACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga        <  1:77634/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:1790879/1‑63 (MQ=255)                                                                 cGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag         >  1:2236511/1‑63 (MQ=255)                                                                  |                                                                     TAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACG  >  minE/2595195‑2595328

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