Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,157 G→T 28.8% A147A (GCC→GCA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1570GT28.8% 100.7 / 28.2 52A147A (GCC→GCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base T (6/9);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCAT  >  minE/2596087‑2596225
                                                                      |                                                                    
cGAGGCCAAAGAAATCAGCGTAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCg                                                                      <  1:3030136/71‑1 (MQ=255)
  aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                                    <  1:1731717/71‑1 (MQ=255)
    gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCaa                                                                  >  1:2084086/1‑71 (MQ=255)
     ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                                 <  1:2061612/71‑1 (MQ=255)
       gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAtaa                                                               <  1:1599426/70‑3 (MQ=255)
        aagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATaaa                                                              <  1:2122001/71‑1 (MQ=255)
        aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa                                                              >  1:333578/1‑71 (MQ=255)
        aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa                                                              >  1:3449102/1‑71 (MQ=255)
         aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt                                                             >  1:3510425/1‑71 (MQ=255)
                  cGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgc                                                     >  1:70191/1‑70 (MQ=255)
                   ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGTTATCGGCAATGAATCCGGCCGCCa                                                   >  1:3568100/1‑71 (MQ=255)
                      gaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAtga                                                <  1:2985791/70‑3 (MQ=255)
                       aCGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgcc                                                       >  1:56520/1‑63 (MQ=255)
                       aCGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCAc                                               <  1:1297938/71‑1 (MQ=255)
                       aCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAtaa                                                               >  1:2573727/1‑53 (MQ=255)
                         gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAt                                                                 <  1:2702354/51‑1 (MQ=255)
                         gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                             <  1:2419478/71‑3 (MQ=255)
                         gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAca                                              <  1:1542783/70‑2 (MQ=255)
                           tatTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAAc                                           >  1:2097649/1‑71 (MQ=255)
                             ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCAc                                               <  1:487693/65‑1 (MQ=255)
                              tCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGaa                                            <  1:3465117/67‑1 (MQ=255)
                              tCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                                         >  1:3485191/1‑70 (MQ=255)
                                aCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCg                                                          <  1:3889235/51‑1 (MQ=255)
                                aCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCa                                       <  1:819963/70‑1 (MQ=255)
                                 ccAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATcc                                                           >  1:991115/1‑49 (MQ=255)
                                 ccAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgt                                     <  1:307056/71‑2 (MQ=255)
                                  cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgcc                                                    <  1:416266/55‑1 (MQ=255)
                                  cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                                    <  1:3783559/71‑1 (MQ=255)
                                   aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTgca                                                       >  1:3436126/1‑50 (MQ=255)
                                   aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgt                                   >  1:2497266/1‑71 (MQ=255)
                                    ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc                                  >  1:3315800/1‑70 (MQ=255)
                                     gTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                                 <  1:1555969/39‑1 (MQ=255)
                                     gTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGc                                                                    <  1:955960/36‑1 (MQ=12)
                                        ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAg                                 >  1:1119177/1‑68 (MQ=255)
                                        ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                               <  1:2402450/70‑1 (MQ=255)
                                         gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGc                              <  1:2267886/70‑1 (MQ=255)
                                          agaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCgg                                                         <  1:58925/42‑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)
                                                                      |                                                                    
CGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCAT  >  minE/2596087‑2596225

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: