Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,545 G→A 8.5% S351S (AGC→AGT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5450GA 8.5% 293.1 / 12.8 94S351S (AGC→AGTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (40/46);  new base A (6/2);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 104 total aligned reads displayed.

TTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGA  >  minE/2595476‑2595609
                                                                     |                                                                
tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                                 >  1:2636073/1‑71 (MQ=255)
 ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                                >  1:2552317/1‑71 (MQ=255)
 ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                                >  1:1971017/1‑71 (MQ=255)
 ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                                >  1:1271173/1‑71 (MQ=255)
 ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGACAGATGCCGtgc                                                                 >  1:2625970/1‑70 (MQ=255)
    ggTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctgc                                                              <  1:1404519/70‑1 (MQ=255)
        ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                          <  1:4034034/70‑1 (MQ=255)
        ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                          >  1:4091261/1‑70 (MQ=255)
          aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                       >  1:3418224/1‑71 (MQ=255)
          aaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                       >  1:4046826/1‑71 (MQ=255)
           aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                      <  1:2535050/71‑1 (MQ=255)
           aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                      >  1:4115897/1‑71 (MQ=255)
            aTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                      >  1:360019/1‑70 (MQ=255)
             tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                      <  1:3067505/69‑1 (MQ=255)
             tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                      <  1:2573050/69‑1 (MQ=255)
             tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                      <  1:766684/69‑1 (MQ=255)
               gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                   <  1:643250/70‑1 (MQ=255)
               gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                   >  1:1839690/1‑70 (MQ=255)
                  gCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcg                                               <  1:3030406/71‑1 (MQ=255)
                   ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                               >  1:3492090/1‑70 (MQ=255)
                     aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                            >  1:489177/1‑71 (MQ=255)
                     aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                            >  1:250142/1‑71 (MQ=255)
                     aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                            >  1:1982458/1‑71 (MQ=255)
                      aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                           <  1:1189527/71‑1 (MQ=255)
                       tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                          <  1:3010525/71‑1 (MQ=255)
                       tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                          <  1:3039604/71‑1 (MQ=255)
                        cacaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                       >  1:578681/1‑57 (MQ=255)
                        cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                                 >  1:1219385/1‑47 (MQ=255)
                         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                        >  1:800071/1‑55 (MQ=255)
                         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                        >  1:402603/1‑55 (MQ=255)
                         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                        >  1:3450223/1‑55 (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)
                              ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                                 >  1:1094884/1‑41 (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:2447130/71‑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)
                                   aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                        >  1:3523999/1‑45 (MQ=255)
                                   aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                              >  1:12138/1‑71 (MQ=255)
                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                      >  1:357494/1‑62 (MQ=255)
                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:81717/1‑71 (MQ=255)
                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:3014557/1‑71 (MQ=255)
                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:3678812/1‑71 (MQ=255)
                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:3349748/1‑71 (MQ=255)
                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:3892744/1‑71 (MQ=255)
                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             <  1:310282/71‑1 (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)
                                        aTCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCCGCat                         <  1:1576686/71‑1 (MQ=255)
                                         tCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                      >  1:2786869/1‑57 (MQ=255)
                                         tCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCata                        <  1:3855249/71‑1 (MQ=255)
                                          cGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCatat                       >  1:1477518/1‑71 (MQ=255)
                                           gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCa                                                         >  1:2743605/1‑36 (MQ=255)
                                           gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                       >  1:3742028/1‑38 (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:3165685/68‑4 (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)
                                                           cAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga       <  1:57482/70‑1 (MQ=255)
                                                             gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag    <  1:902455/71‑1 (MQ=255)
                                                             gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag    <  1:897302/71‑1 (MQ=255)
                                                             gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag    <  1:273273/71‑1 (MQ=255)
                                                              tcgCCGTACTCCCGTCAATCGACAGCGCACCAATCATCACCGTCGGCATATTGTTAATAACTGAtgaga     >  1:1648825/3‑69 (MQ=255)
                                                              aTGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga       <  1:2344369/67‑1 (MQ=255)
                                                               tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:2814061/1‑44 (MQ=255)
                                                               tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagaga     >  1:265722/1‑68 (MQ=255)
                                                               tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa  <  1:3914127/71‑1 (MQ=255)
                                                               tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGCCATATTGTTCATAATAGAAGAGAgg   >  1:3063473/1‑70 (MQ=255)
                                                               cgCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGAtgagag  >  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)
                                                                     |                                                                
TTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGA  >  minE/2595476‑2595609

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: