Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,788 A→G 6.7% N270N (AAT→AAC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7880AG 6.7% 211.4 / 6.2 74N270N (AAT→AACarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base G (2/3);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

AAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACC  >  minE/2595718‑2595857
                                                                      |                                                                     
aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                                                       >  1:1050198/1‑71 (MQ=255)
  accaccAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAtt                                                                     <  1:3981654/71‑1 (MQ=255)
   ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                                    <  1:3118844/71‑1 (MQ=255)
    caccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                                   >  1:1221211/1‑71 (MQ=255)
       caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                                >  1:3218807/1‑71 (MQ=255)
       caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                                <  1:873797/71‑1 (MQ=255)
        aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                                                       >  1:2292813/1‑63 (MQ=255)
         gATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                              >  1:479754/1‑71 (MQ=255)
           tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                                   >  1:1446015/1‑64 (MQ=255)
           tACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAAt                                                             <  1:3560186/70‑1 (MQ=255)
           tACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATg                                                            <  1:1272004/71‑1 (MQ=255)
            aCATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                           >  1:1143159/1‑71 (MQ=255)
             cATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                                <  1:1546797/65‑1 (MQ=255)
                 ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                      >  1:3715574/1‑71 (MQ=255)
                 ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                      >  1:1150191/1‑71 (MQ=255)
                 ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                      <  1:3724475/71‑1 (MQ=255)
                 ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgc                                                       <  1:1239249/70‑1 (MQ=255)
                  cgagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgc                                                       >  1:1707907/1‑69 (MQ=255)
                   gagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGtt                                                    >  1:2354371/1‑71 (MQ=255)
                    agcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcg                                                        <  1:542549/66‑1 (MQ=255)
                     gcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAtt                                                                     >  1:4146121/1‑52 (MQ=255)
                      tgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGttg                                                                    <  1:1255131/51‑2 (MQ=255)
                       gagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGc                                                >  1:2720152/1‑71 (MQ=255)
                             aTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                               <  1:3595693/66‑1 (MQ=255)
                               gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac                                             >  1:3410138/1‑66 (MQ=255)
                               gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                         >  1:2972399/1‑70 (MQ=255)
                                aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                       >  1:3746971/1‑71 (MQ=255)
                                 cAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                       >  1:3127947/1‑70 (MQ=255)
                                   aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                                   <  1:2015137/40‑1 (MQ=255)
                                   aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                    <  1:649158/71‑1 (MQ=255)
                                    tCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                    >  1:2263615/1‑70 (MQ=255)
                                     cTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                                    >  1:424177/1‑37 (MQ=255)
                                      tGCCAGGGGGCACCGCTCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                                 >  1:1346485/1‑71 (MQ=255)
                                      tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCAt                                                             >  1:2794700/1‑43 (MQ=255)
                                      tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                    >  1:339359/1‑68 (MQ=255)
                                      tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                                 >  1:1679962/1‑71 (MQ=255)
                                      tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                                 >  1:3684755/1‑71 (MQ=255)
                                       gCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcga                                              >  1:63436/1‑57 (MQ=255)
                                        ccAGTGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                               <  1:2647464/55‑1 (MQ=255)
                                        ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                               <  1:349721/55‑1 (MQ=255)
                                        ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                               <  1:1266369/55‑1 (MQ=255)
                                        ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                   <  1:1589338/67‑1 (MQ=255)
                                        ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAataata                                  <  1:4052833/68‑1 (MQ=255)
                                        ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCa                               <  1:959771/71‑1 (MQ=255)
                                         cAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAg                              <  1:936051/71‑1 (MQ=255)
                                           gggggCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAgc                             <  1:70213/70‑1 (MQ=255)
                                               gCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGcc                        <  1:3954906/71‑1 (MQ=255)
                                               gCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGcc                        <  1:472007/71‑1 (MQ=255)
                                                  ccGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGcca                                              >  1:16448/1‑44 (MQ=255)
                                                  ccGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaata                                     <  1:937985/55‑1 (MQ=255)
                                                  ccGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaata                                     <  1:902947/55‑1 (MQ=255)
                                                   cgcgCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcga                                              >  1:1090519/1‑45 (MQ=255)
                                                   cgcgCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                   >  1:3025188/1‑56 (MQ=255)
                                                   cgcgCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                   >  1:2578243/1‑56 (MQ=255)
                                                    gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGc                                                <  1:2715114/42‑1 (MQ=255)
                                                    gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                         >  1:2985067/1‑49 (MQ=255)
                                                        aGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGc               <  1:2579417/71‑1 (MQ=255)
                                                          gACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCaa             <  1:787428/71‑1 (MQ=255)
                                                            cTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATg           >  1:2107911/1‑71 (MQ=255)
                                                             tttACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                   >  1:2183796/1‑46 (MQ=255)
                                                             tttACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCaa             <  1:2987129/68‑1 (MQ=255)
                                                              ttACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGgc         >  1:1996529/1‑71 (MQ=255)
                                                               tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                       >  1:3192892/1‑40 (MQ=255)
                                                               tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                       >  1:666368/1‑40 (MQ=255)
                                                               tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                       >  1:3010115/1‑40 (MQ=255)
                                                               tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGgcg        >  1:3088745/1‑71 (MQ=255)
                                                               tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGgcg        >  1:1553068/1‑71 (MQ=255)
                                                                 cccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCt      <  1:2591/71‑1 (MQ=255)
                                                                 cccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCt      <  1:2719171/71‑1 (MQ=255)
                                                                     tATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                   <  1:1974810/38‑1 (MQ=255)
                                                                     tATTAATCGCATGACCGCGTTTAGCGACGACAaat                                      <  1:1171590/35‑1 (MQ=255)
                                                                     tATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAgcgc                           <  1:3293939/46‑1 (MQ=255)
                                                                     tATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGc               <  1:3817591/58‑1 (MQ=255)
                                                                     tATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGc               <  1:1692272/58‑1 (MQ=255)
                                                                     tATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGc               <  1:1201229/58‑1 (MQ=255)
                                                                     tATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCAcc  >  1:379657/1‑71 (MQ=255)
                                                                      aTTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                                 >  1:4092351/1‑39 (MQ=255)
                                                                      aTTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                                 >  1:72124/1‑39 (MQ=255)
                                                                      |                                                                     
AAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACC  >  minE/2595718‑2595857

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: