Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,184 C→G 10.3% A138A (GCG→GCC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1840CG10.3% 159.9 / 7.9 58A138A (GCG→GCCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (24/28);  new base G (4/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAAGCG  >  minE/2596114‑2596254
                                                                      |                                                                      
tatTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAAc                                                                        >  1:2097649/1‑71 (MQ=255)
   tCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                                                                      >  1:3485191/1‑70 (MQ=255)
     aCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCa                                                                    <  1:819963/70‑1 (MQ=255)
      ccAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgt                                                                  <  1:307056/71‑2 (MQ=255)
       cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                                                                 <  1:3783559/71‑1 (MQ=255)
        aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgt                                                                >  1:2497266/1‑71 (MQ=255)
         ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc                                                               >  1:3315800/1‑70 (MQ=255)
             ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAg                                                              >  1:1119177/1‑68 (MQ=255)
             ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                                                            <  1:2402450/70‑1 (MQ=255)
              gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGc                                                           <  1:2267886/70‑1 (MQ=255)
               agaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCt                                                         <  1:1856390/71‑1 (MQ=255)
               aaaCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACATAGGCCAGCGTCGTGccc                                                         >  1:548452/3‑70 (MQ=255)
                gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGt                                                                <  1:3057455/63‑1 (MQ=255)
                 aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                                       <  1:2091714/71‑4 (MQ=255)
                   aaTAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                                                          <  1:2481318/66‑1 (MQ=255)
                    aTAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                                    <  1:868664/71‑1 (MQ=255)
                        gCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAc                                                <  1:3213076/71‑1 (MQ=255)
                           gCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                                    >  1:1278145/1‑64 (MQ=255)
                           gCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAccc                                              >  1:2921311/1‑70 (MQ=255)
                            cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAcc                                               >  1:1064160/1‑68 (MQ=255)
                            cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                                            <  1:857920/71‑1 (MQ=255)
                            cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                                            <  1:988578/71‑1 (MQ=255)
                             aGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                                            >  1:2518166/1‑70 (MQ=255)
                              ggCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                                            >  1:2130549/1‑69 (MQ=255)
                                cTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc                                        <  1:3818098/71‑1 (MQ=255)
                                    cGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                                     >  1:2087774/1‑70 (MQ=255)
                                           ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCAt                               <  1:2851223/69‑1 (MQ=255)
                                            gCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc                            >  1:1089642/1‑71 (MQ=255)
                                            gCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc                            <  1:44490/71‑1 (MQ=255)
                                             cAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCg                           <  1:891881/71‑1 (MQ=255)
                                              aaTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGa                          >  1:1419448/1‑71 (MQ=255)
                                              aaTGAATCAGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGa                          >  1:1457080/1‑71 (MQ=255)
                                               aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                                    <  1:848298/44‑1 (MQ=255)
                                               aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                                    <  1:3698868/44‑1 (MQ=255)
                                               aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt                         >  1:3428531/1‑71 (MQ=255)
                                               aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt                         >  1:94696/1‑71 (MQ=255)
                                               aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt                         >  1:138630/1‑71 (MQ=255)
                                                 gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                                     >  1:3607990/1‑57 (MQ=255)
                                                  aaTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGAc                      <  1:257306/71‑1 (MQ=255)
                                                   aTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc                            >  1:2076347/1‑64 (MQ=255)
                                                      cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc                  <  1:1870091/71‑1 (MQ=255)
                                                          cgccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                                            <  1:2485358/41‑1 (MQ=255)
                                                             cATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt                         <  1:2019084/57‑1 (MQ=255)
                                                              aTCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc                                        >  1:860097/1‑41 (MQ=255)
                                                              aTCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaa          <  1:1472319/71‑1 (MQ=255)
                                                              aTCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCa           >  1:1787107/1‑70 (MQ=255)
                                                              aTCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCa           >  1:3835687/1‑70 (MQ=255)
                                                               tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                                            <  1:1069713/36‑1 (MQ=255)
                                                               tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGa                       >  1:115348/1‑57 (MQ=255)
                                                               tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaa         <  1:2241733/71‑1 (MQ=255)
                                                                 aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcagc                                    <  1:1064877/42‑1 (MQ=255)
                                                                 aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAAt       <  1:861728/71‑1 (MQ=255)
                                                                 aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAAt       <  1:1712926/71‑1 (MQ=255)
                                                                  cGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAAt       >  1:300887/1‑70 (MQ=255)
                                                                  cGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATa      <  1:2034473/71‑1 (MQ=255)
                                                                   gAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAAt       >  1:2590151/1‑69 (MQ=255)
                                                                   gAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAAt       >  1:2445719/1‑69 (MQ=255)
                                                                     aCGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAAgc   >  1:3256998/1‑71 (MQ=255)
                                                                     aCGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAAg    >  1:1548270/1‑70 (MQ=255)
                                                                      cGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAAgcg  <  1:869576/71‑1 (MQ=255)
                                                                      |                                                                      
TATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAAGCG  >  minE/2596114‑2596254

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: