Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,163 G→A 18.5% F145F (TTC→TTT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1630GA18.5% 117.4 / 19.4 54F145F (TTC→TTTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base A (5/5);  total (26/28)
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.39e-01

AAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAT  >  minE/2596095‑2596231
                                                                    |                                                                    
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)
                 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)
                                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:988578/71‑1 (MQ=255)
                                               cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                     <  1:857920/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:44490/71‑1 (MQ=255)
                                                               gCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc     >  1:1089642/1‑71 (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:138630/1‑71 (MQ=255)
                                                                  aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt  >  1:3428531/1‑71 (MQ=255)
                                                                  aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt  >  1:94696/1‑71 (MQ=255)
                                                                    gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag              >  1:3607990/1‑57 (MQ=255)
                                                                    |                                                                    
AAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAT  >  minE/2596095‑2596231

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: