Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,551 A→G 26.7% D349D (GAT→GAC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5510AG26.7% 253.8 / 108.6 165D349D (GAT→GACarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (65/56);  new base G (28/16);  total (93/72)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 165 total aligned reads displayed.

TCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAG  >  minE/2595482‑2595619
                                                                     |                                                                    
 cctaaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTc                                                                    <  1:2033291/68‑1 (MQ=255)
 cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                                                     >  1:763433/1‑70 (MQ=255)
  ctaaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTc                                                                    >  1:2891717/3‑70 (MQ=255)
  ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                                    <  1:3367170/70‑1 (MQ=255)
   cAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                                  <  1:3422265/71‑1 (MQ=255)
    aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                                 <  1:833011/71‑1 (MQ=255)
    aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                                 <  1:4455755/71‑1 (MQ=255)
    aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                                 >  1:2927406/1‑71 (MQ=255)
     aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                                 >  1:2903195/1‑70 (MQ=255)
       tcacAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                              >  1:2392320/4‑71 (MQ=255)
       tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                                              >  1:4288879/1‑71 (MQ=255)
         gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                             <  1:3386210/70‑1 (MQ=255)
         gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                             <  1:3255432/70‑1 (MQ=255)
         acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGAc                                                             <  1:1526969/69‑1 (MQ=255)
         acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGAc                                                             <  1:3159032/69‑1 (MQ=255)
         acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACa                                                            <  1:214923/70‑1 (MQ=255)
          cAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACag                                                           <  1:1179477/71‑2 (MQ=255)
          cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                                 <  1:84031/65‑1 (MQ=255)
            gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                                         >  1:4631505/1‑71 (MQ=255)
                aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                              >  1:3782538/1‑62 (MQ=255)
                aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                              >  1:2181547/1‑62 (MQ=255)
                aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                              >  1:2777726/1‑62 (MQ=255)
                aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                              >  1:1608917/1‑62 (MQ=255)
                aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                     <  1:2079645/71‑1 (MQ=255)
                aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                     <  1:5071896/71‑1 (MQ=255)
                 tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                              >  1:4010419/1‑61 (MQ=255)
                 tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                              >  1:2594187/1‑61 (MQ=255)
                 tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                                    <  1:2781269/71‑1 (MQ=255)
                  cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATgg                                                               >  1:5564043/1‑59 (MQ=255)
                   acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                                                       >  1:223661/1‑65 (MQ=255)
                   acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                                                       >  1:3832121/1‑65 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                                   >  1:3466013/1‑70 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                  <  1:3254161/71‑1 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                  <  1:5472608/71‑1 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                  >  1:157761/1‑71 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                  >  1:3459158/1‑71 (MQ=255)
                    caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAcc                                                 <  1:1634087/71‑1 (MQ=255)
                    caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                  <  1:1966156/70‑1 (MQ=255)
                     aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                                              >  1:168534/1‑57 (MQ=255)
                     aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAc                                                          <  1:2983404/61‑1 (MQ=255)
                     aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                                <  1:4849085/71‑1 (MQ=255)
                      ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCaa                                                                  >  1:79368/1‑52 (MQ=255)
                      ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                                 >  1:3608310/1‑53 (MQ=255)
                      ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                                 >  1:5379357/1‑53 (MQ=255)
                      ttGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg                                               <  1:3909527/71‑1 (MQ=255)
                       tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                                  >  1:2678707/1‑51 (MQ=255)
                       tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                                              >  1:1689067/1‑55 (MQ=255)
                       tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                                   >  1:3055542/1‑66 (MQ=255)
                       tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                                   >  1:3920017/1‑66 (MQ=255)
                        ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                                  >  1:3308327/1‑50 (MQ=255)
                         gCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                             <  1:2396974/70‑1 (MQ=255)
                            tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                <  1:135212/48‑1 (MQ=255)
                            tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                                           >  1:5394825/1‑53 (MQ=255)
                             aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc                                            >  1:1249439/1‑67 (MQ=255)
                             aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc                                            >  1:1088319/1‑67 (MQ=255)
                             aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                        >  1:2560031/1‑71 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                       >  1:5035245/1‑55 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                       >  1:3736814/1‑55 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                       >  1:611440/1‑55 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                        <  1:4005571/70‑1 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                       >  1:4762948/1‑71 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                       >  1:4571849/1‑71 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                       >  1:190643/1‑71 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGAGCTGCCATCAATGGACAACGCGCCAAcc                                                 >  1:4243777/1‑61 (MQ=255)
                                 cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                                    <  1:769114/39‑1 (MQ=255)
                                  aTCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                                  >  1:4003632/1‑40 (MQ=255)
                                     gcctcttTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCatat                                 >  1:300896/4‑70 (MQ=255)
                                     gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                                 >  1:5639749/1‑38 (MQ=255)
                                      cctcttTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTg                               <  1:4234961/69‑1 (MQ=255)
                                         cttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                                                        <  1:5218190/43‑1 (MQ=255)
                                         cttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                              <  1:1743118/69‑1 (MQ=255)
                                          tttGATAACGCCAGATGCCGTGCTGCCATCAATGGACa                                                            >  1:2753541/1‑38 (MQ=255)
                                           ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                      >  1:367765/1‑43 (MQ=255)
                                            tGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                         >  1:389031/1‑55 (MQ=255)
                                             gATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                        >  1:900912/1‑71 (MQ=255)
                                             gATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                        >  1:467426/1‑71 (MQ=255)
                                              aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                              >  1:5610185/1‑64 (MQ=255)
                                               cgaCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTc                            <  1:2852098/63‑1 (MQ=255)
                                                 aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                          <  1:5301797/65‑1 (MQ=255)
                                                 aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATaga                    <  1:1441744/71‑1 (MQ=255)
                                                 aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATaga                    <  1:1448564/71‑1 (MQ=255)
                                                 aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                     <  1:1061604/70‑1 (MQ=255)
                                                 aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                     <  1:4980414/70‑1 (MQ=255)
                                                  cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa                   <  1:464555/71‑1 (MQ=255)
                                                   gCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCatca                        >  1:1893540/1‑63 (MQ=255)
                                                   gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                                >  1:2004573/1‑41 (MQ=255)
                                                   gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaag                  <  1:1003764/71‑1 (MQ=255)
                                                    ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                           <  1:2483012/45‑1 (MQ=255)
                                                       gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag              <  1:5283863/71‑1 (MQ=255)
                                                          gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                        >  1:5695721/1‑42 (MQ=255)
                                                          gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                        >  1:1026562/1‑42 (MQ=255)
                                                          gCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGt                              >  1:2675742/1‑52 (MQ=255)
                                                           ccGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGtt                             >  1:4158798/1‑52 (MQ=255)
                                                           ccGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                          <  1:1134433/55‑1 (MQ=255)
                                                           ccGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                          <  1:4648604/55‑1 (MQ=255)
                                                              tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCgg                                      <  1:2187126/40‑1 (MQ=255)
                                                              tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCg       <  1:2098620/71‑1 (MQ=255)
                                                                ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCat                                   >  1:569487/1‑41 (MQ=255)
                                                                   acATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAAAAGAAGAGAGGAAGGCGGTGAg  >  1:2566778/2‑71 (MQ=255)
                                                                    cATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                        >  1:3032727/1‑48 (MQ=255)
                                                                     |                                                                    
TCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAG  >  minE/2595482‑2595619

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: