Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,548 G→C 27.7% G350G (GGC→GGG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5480GC27.7% 243.7 / 112.8 159G350G (GGC→GGGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (62/53);  new base C (28/16);  total (90/69)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 164 total aligned reads displayed.

TCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCG  >  minE/2595482‑2595614
                                                                  |                                                                  
 cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                                                >  1:763433/1‑70 (MQ=255)
 cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                                                >  1:3396011/1‑70 (MQ=255)
  ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                               <  1:1722350/70‑1 (MQ=255)
  ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                               <  1:5299032/70‑1 (MQ=255)
  ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                               <  1:3949922/70‑1 (MQ=255)
  ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                               <  1:3367170/70‑1 (MQ=255)
   taaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCaa                                                             >  1:3565030/2‑71 (MQ=255)
   cAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                             <  1:3422265/71‑1 (MQ=255)
    aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                            <  1:833011/71‑1 (MQ=255)
    aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                            <  1:4455755/71‑1 (MQ=255)
     aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                           <  1:334155/71‑1 (MQ=255)
     aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                                                           <  1:5437133/71‑2 (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:4587681/70‑1 (MQ=255)
         acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                         >  1:2101087/2‑69 (MQ=255)
          cAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACag                                                      <  1:1179477/71‑2 (MQ=255)
          cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                            <  1:84031/65‑1 (MQ=255)
            gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                                    <  1:2046393/71‑1 (MQ=255)
            gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                                    >  1:4631505/1‑71 (MQ=255)
            gCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcg                                                    <  1:2923227/71‑1 (MQ=255)
             ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgc                                                   <  1:3262228/71‑1 (MQ=255)
                aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                         >  1:3782538/1‑62 (MQ=255)
                aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                         >  1:2181547/1‑62 (MQ=255)
                aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                         >  1:2176540/1‑62 (MQ=255)
                aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                <  1:2079645/71‑1 (MQ=255)
                aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                <  1:5071896/71‑1 (MQ=255)
                 tCACATTGGCATAAATCATTGCCTCTTTGACGCCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCAcc                                                >  1:2467144/1‑70 (MQ=255)
                 tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                         >  1:2221376/1‑61 (MQ=255)
                 tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                         >  1:2594187/1‑61 (MQ=255)
                 tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                               <  1:2781269/71‑1 (MQ=255)
                  cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctgc                                                                   >  1:1351659/1‑50 (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:466127/70‑1 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                             <  1:3254161/71‑1 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                             >  1:3459158/1‑71 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                             <  1:5472608/71‑1 (MQ=255)
                   acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                             >  1:157761/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)
                        ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                             >  1:3142306/1‑50 (MQ=255)
                        ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                                         <  1:1472548/70‑1 (MQ=255)
                        ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                                         <  1:4262823/70‑1 (MQ=255)
                             aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc                                       >  1:1249439/1‑67 (MQ=255)
                             aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc                                       >  1:1088319/1‑67 (MQ=255)
                             aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                   >  1:2560031/1‑71 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                                      >  1:2762815/1‑51 (MQ=255)
                              aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                  >  1:5035245/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:1911825/1‑71 (MQ=255)
                              aaCCATAGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                  >  1:3215243/1‑55 (MQ=255)
                                 cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                               <  1:4131968/39‑1 (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)
                                     gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                           >  1:285528/1‑71 (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)
                                           ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                 >  1:1388473/1‑43 (MQ=255)
                                            tGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCat                     <  1:278853/70‑1 (MQ=255)
                                             gATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                   >  1:900912/1‑71 (MQ=255)
                                             gATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                   >  1:467426/1‑71 (MQ=255)
                                              aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                         >  1:5610185/1‑64 (MQ=255)
                                                 aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                     <  1:5301797/65‑1 (MQ=255)
                                                 aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATaga               <  1:1441744/71‑1 (MQ=255)
                                                 aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                <  1:4980414/70‑1 (MQ=255)
                                                 aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                <  1:1061604/70‑1 (MQ=255)
                                                  cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa              <  1:3890566/71‑1 (MQ=255)
                                                  cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa              <  1:464555/71‑1 (MQ=255)
                                                   gCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCatca                   >  1:1893540/1‑63 (MQ=255)
                                                   gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                           >  1:2674055/1‑41 (MQ=255)
                                                   gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                           >  1:2004573/1‑41 (MQ=255)
                                                   gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                           >  1:1560367/1‑41 (MQ=255)
                                                   gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaag             <  1:1003764/71‑1 (MQ=255)
                                                    ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                      <  1:2483012/45‑1 (MQ=255)
                                                    ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga            <  1:1625058/71‑1 (MQ=255)
                                                          gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                   >  1:3029323/1‑42 (MQ=255)
                                                          gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                   >  1:5695721/1‑42 (MQ=255)
                                                          gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                   >  1:1026562/1‑42 (MQ=255)
                                                           ccGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGtt                        >  1:4158798/1‑52 (MQ=255)
                                                           ccGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                     <  1:1134433/55‑1 (MQ=255)
                                                              tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCgg                                 <  1:3779870/40‑1 (MQ=255)
                                                                ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCat                              >  1:569487/1‑41 (MQ=255)
                                                                  |                                                                  
TCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCG  >  minE/2595482‑2595614

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: