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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1810GA10.7% 154.5 / 8.0 56F139F (TTC→TTTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (21/29);  new base A (4/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATA  >  minE/2596114‑2596250
                                                                   |                                                                     
tatTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAAc                                                                    >  1:2097649/1‑71 (MQ=255)
   tCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGaa                                                                     <  1:3465117/67‑1 (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:94696/1‑71 (MQ=255)
                                               aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt                     >  1:138630/1‑71 (MQ=255)
                                               aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt                     >  1:3428531/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:3835687/1‑70 (MQ=255)
                                                              aTCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCa       >  1:1787107/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:1712926/71‑1 (MQ=255)
                                                                 aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAAt   <  1:861728/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)
                                                                   |                                                                     
TATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATA  >  minE/2596114‑2596250

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: