Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,178 C→G 14.5% V140V (GTG→GTC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1780CG14.5% 127.3 / 13.9 55V140V (GTG→GTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (20/27);  new base G (4/4);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

ACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAA  >  minE/2596110‑2596247
                                                                    |                                                                     
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)
         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: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:3698868/44‑1 (MQ=255)
                                                   aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                             <  1:848298/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: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)
                                                                    |                                                                     
ACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAA  >  minE/2596110‑2596247

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: