Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,557 G→C 9.9% S347S (TCC→TCG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5570GC 9.9% 284.3 / 16.8 100S347S (TCC→TCGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base C (8/2);  total (54/46)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 108 total aligned reads displayed.

AATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCG  >  minE/2595487‑2595626
                                                                      |                                                                     
aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                       <  1:2535050/71‑1 (MQ=255)
aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                       >  1:4115897/1‑71 (MQ=255)
 aTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                       >  1:360019/1‑70 (MQ=255)
  tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                       <  1:766684/69‑1 (MQ=255)
  tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                       <  1:2573050/69‑1 (MQ=255)
  tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                       <  1:3067505/69‑1 (MQ=255)
    gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                                    >  1:1839690/1‑70 (MQ=255)
    gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                                    <  1:643250/70‑1 (MQ=255)
       gCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcg                                                                <  1:3030406/71‑1 (MQ=255)
        ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                                                >  1:3492090/1‑70 (MQ=255)
          aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                             >  1:250142/1‑71 (MQ=255)
          aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                             >  1:489177/1‑71 (MQ=255)
          aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                             >  1:1982458/1‑71 (MQ=255)
           aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                            <  1:1189527/71‑1 (MQ=255)
            tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                                           <  1:3039604/71‑1 (MQ=255)
            tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                                           <  1:3010525/71‑1 (MQ=255)
              acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                         >  1:3792116/1‑71 (MQ=255)
              acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                         >  1:1472393/1‑71 (MQ=255)
               caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                         <  1:2885897/70‑1 (MQ=255)
                aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                            >  1:3042308/1‑66 (MQ=255)
                aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                                       <  1:2423512/71‑1 (MQ=255)
                  tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg                                                      <  1:1339796/70‑1 (MQ=255)
                   ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                                    >  1:843990/1‑71 (MQ=255)
                   ggCATAAACCATCGCTTCTTTGATAACGCCAGAAGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                                    >  1:2676502/1‑71 (MQ=255)
                     cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                         <  1:920206/64‑1 (MQ=255)
                     cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                  <  1:2098857/71‑1 (MQ=255)
                     cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                  <  1:2578237/71‑1 (MQ=255)
                     cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                  <  1:668832/71‑1 (MQ=255)
                       tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                       <  1:495801/48‑1 (MQ=255)
                       tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                                       <  1:103480/48‑1 (MQ=255)
                       tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                         <  1:190327/62‑1 (MQ=255)
                       tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                  <  1:3149685/69‑1 (MQ=255)
                       tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                                >  1:3616433/1‑71 (MQ=255)
                       tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                                >  1:2218321/1‑71 (MQ=255)
                        aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                               >  1:12138/1‑71 (MQ=255)
                         aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                                       >  1:357494/1‑62 (MQ=255)
                         aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                              >  1:3678812/1‑71 (MQ=255)
                         aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                              >  1:3014557/1‑71 (MQ=255)
                         aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                              <  1:310282/71‑1 (MQ=255)
                         aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                              >  1:3349748/1‑71 (MQ=255)
                         aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                              >  1:3892744/1‑71 (MQ=255)
                         aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                              >  1:81717/1‑71 (MQ=255)
                         aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGAACAACGCGCCAACCAGTACCGTCg                                              <  1:2289317/71‑1 (MQ=255)
                          aCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                                    <  1:3931714/48‑1 (MQ=255)
                          aCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                                    <  1:837545/48‑1 (MQ=255)
                            cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCa                                           <  1:2187272/71‑1 (MQ=255)
                             aTCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                                    <  1:2311672/45‑1 (MQ=255)
                              tCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                                       >  1:2786869/1‑57 (MQ=255)
                              tCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCata                                         <  1:3855249/71‑1 (MQ=255)
                               cGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCatat                                        >  1:1477518/1‑71 (MQ=255)
                                gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                                                  >  1:2183550/1‑44 (MQ=255)
                                gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                                       >  1:3478741/1‑71 (MQ=255)
                                gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                                       >  1:2203847/1‑71 (MQ=255)
                                 cttcttTGATAACGCCAGATGCCGTGCTTCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                                      <  1:1520909/71‑1 (MQ=255)
                                 cttcttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                                      <  1:1238127/71‑1 (MQ=255)
                                      ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCa                                  >  1:1295209/1‑70 (MQ=255)
                                       tGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCat                                 <  1:1604755/70‑1 (MQ=255)
                                         aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                <  1:2527304/69‑1 (MQ=255)
                                           aaCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                                     >  1:2260972/1‑62 (MQ=255)
                                            aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCatcac                              <  1:2500424/68‑4 (MQ=255)
                                             cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa                          <  1:477469/71‑1 (MQ=255)
                                              gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCat                                          <  1:424490/54‑1 (MQ=255)
                                               ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                        <  1:1733267/71‑1 (MQ=255)
                                                  gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                     <  1:273273/71‑1 (MQ=255)
                                                  gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                     <  1:897302/71‑1 (MQ=255)
                                                  gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                     <  1:902455/71‑1 (MQ=255)
                                                   tcgCCGTACTCCCGTCAATCGACAGCGCACCAATCATCACCGTCGGCATATTGTTAATAACTGAtgaga                      >  1:1648825/3‑69 (MQ=255)
                                                   aTGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                        <  1:2344369/67‑1 (MQ=255)
                                                    tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagaga                      >  1:265722/1‑68 (MQ=255)
                                                    tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa                   <  1:3914127/71‑1 (MQ=255)
                                                    tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGCCATATTGTTCATAATAGAAGAGAgg                    >  1:3063473/1‑70 (MQ=255)
                                                    cgCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGAtgagaga                    >  1:301409/2‑69 (MQ=255)
                                                     gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAgg                    >  1:2012762/1‑69 (MQ=255)
                                                     gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa                   <  1:2279910/70‑1 (MQ=255)
                                                        gTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCata                                         >  1:1538798/1‑45 (MQ=255)
                                                         tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataata                             >  1:1029990/1‑56 (MQ=255)
                                                           ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                >  1:3794525/1‑51 (MQ=255)
                                                           ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt            >  1:1287844/1‑71 (MQ=255)
                                                           ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt            <  1:2412370/71‑1 (MQ=255)
                                                                aTCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt            <  1:3900776/66‑1 (MQ=255)
                                                                 tCAATGGACAACGCGCCAACCAGTACCGTCGGGATATTGTTCata                                >  1:2555150/1‑45 (MQ=255)
                                                                 tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                >  1:393022/1‑45 (MQ=255)
                                                                 tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                     <  1:2739145/56‑1 (MQ=255)
                                                                 tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGa          >  1:3784368/1‑67 (MQ=255)
                                                                 tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAgg        <  1:2350959/69‑1 (MQ=255)
                                                                 tCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTc                                   >  1:3653462/1‑42 (MQ=19)
                                                                  cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGaa                  >  1:4045008/1‑58 (MQ=255)
                                                                  cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGaa                  >  1:1597676/1‑58 (MQ=255)
                                                                  cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa                   >  1:2060411/1‑57 (MQ=255)
                                                                  cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa                   >  1:2348929/1‑57 (MQ=255)
                                                                   aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                               >  1:3300282/1‑44 (MQ=255)
                                                                   aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCa                                  >  1:1483563/1‑41 (MQ=255)
                                                                   aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGa          >  1:740986/1‑65 (MQ=255)
                                                                   aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATc    >  1:3072922/1‑71 (MQ=255)
                                                                   aaTCGACAGCGCACCAATCAGCACCGTCGGCATAtt                                       >  1:1908294/1‑36 (MQ=16)
                                                                   aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGAtgagagaa                   >  1:698983/1‑54 (MQ=255)
                                                                   aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGAtgagagaa                   >  1:2820743/1‑54 (MQ=255)
                                                                     tGGATAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAgg                <  1:581612/57‑1 (MQ=255)
                                                                     tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCg  >  1:1849194/1‑71 (MQ=255)
                                                                     tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAAAAGAAGAGAGGAAGGCGGTGAGGAATCCg  >  1:4035344/1‑71 (MQ=255)
                                                                      |                                                                     
AATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCG  >  minE/2595487‑2595626

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: