Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,791 A→G 5.5% I269I (ATT→ATC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7910AG 5.5% 206.1 / 3.4 73I269I (ATT→ATCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base G (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAA  >  minE/2595721‑2595861
                                                                      |                                                                      
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)
      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:3724475/71‑1 (MQ=255)
              ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                          >  1:1150191/1‑71 (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)
                    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:3684755/1‑71 (MQ=255)
                                   tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                                     >  1:1679962/1‑71 (MQ=255)
                                    gCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcga                                                  >  1:63436/1‑57 (MQ=255)
                                     ccAGTGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                                   <  1:2647464/55‑1 (MQ=255)
                                     ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                                   <  1:1266369/55‑1 (MQ=255)
                                     ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                                   <  1:349721/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:472007/71‑1 (MQ=255)
                                            gCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGcc                            <  1:3954906/71‑1 (MQ=255)
                                               ccGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGcca                                                  >  1:16448/1‑44 (MQ=255)
                                               ccGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaata                                         <  1:902947/55‑1 (MQ=255)
                                               ccGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaata                                         <  1:937985/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:2719171/71‑1 (MQ=255)
                                                              cccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCt          <  1:2591/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:1692272/58‑1 (MQ=255)
                                                                  tATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGc                   <  1:3817591/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)
                                                                    ttAATCGCATGACCGCGTTTAGCGACGACAAataata                                      <  1:1686498/37‑1 (MQ=255)
                                                                    ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                                    >  1:68454/1‑39 (MQ=255)
                                                                    ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGcc                            <  1:2274176/47‑1 (MQ=255)
                                                                    ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGc                   <  1:2450024/56‑1 (MQ=255)
                                                                    ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCg     >  1:426050/1‑70 (MQ=255)
                                                                      aaTCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCt                     >  1:1506822/1‑52 (MQ=255)
                                                                      aaTCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGaa  >  1:709439/1‑71 (MQ=255)
                                                                      aaTCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGaa  >  1:1154954/1‑71 (MQ=255)
                                                                      aaTCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGaa  >  1:1099487/1‑71 (MQ=255)
                                                                      |                                                                      
CCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAA  >  minE/2595721‑2595861

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: