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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1690GT19.2% 134.1 / 15.3 52A143A (GCC→GCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base T (5/5);  total (24/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.95e-01

CGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATC  >  minE/2596105‑2596238
                                                                |                                                                     
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)
               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: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:3698868/44‑1 (MQ=255)
                                                        aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                    <  1:848298/44‑1 (MQ=255)
                                                        aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         >  1:3428531/1‑71 (MQ=255)
                                                        aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         >  1:138630/1‑71 (MQ=255)
                                                        aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         >  1:94696/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)
                                                                |                                                                     
CGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATC  >  minE/2596105‑2596238

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: