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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,0940AG10.1% 236.8 / 13.5 89F168F (TTT→TTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (32/48);  new base G (6/3);  total (38/51)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGGCAACAATCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGA  >  minE/2596027‑2596164
                                                                    |                                                                      
gCGTGGCAACAATCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCaa                                                                       <  1:1658216/69‑1 (MQ=255)
     gCAACAATCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATc                                                                <  1:4271830/71‑1 (MQ=255)
        acaaTCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGc                                                             >  1:4656980/1‑71 (MQ=255)
           aTCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGa                                                          <  1:350982/71‑1 (MQ=255)
              gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                        <  1:4620175/70‑1 (MQ=255)
              gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                        <  1:4626850/70‑1 (MQ=255)
                ggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGa                                                     <  1:1529366/71‑1 (MQ=255)
                  cgATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                   >  1:5208722/1‑71 (MQ=255)
                     tatCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGa                                                          <  1:2569275/61‑1 (MQ=255)
                      atCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCa                                               >  1:3750295/1‑71 (MQ=255)
                       tCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                   >  1:3468187/1‑66 (MQ=255)
                       tCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                   >  1:4212338/1‑66 (MQ=255)
                       tCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                   >  1:4246305/1‑66 (MQ=255)
                       tCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTc                                                >  1:2542745/1‑69 (MQ=255)
                         cACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGaa                                                                   <  1:1592809/48‑1 (MQ=255)
                         cACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAAt                                                                 >  1:3275166/1‑50 (MQ=255)
                         cACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAcc                                             <  1:4238865/70‑1 (MQ=255)
                         cACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCa                                            <  1:2734764/71‑1 (MQ=255)
                          aCCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGaaa                                                                  >  1:35870/1‑48 (MQ=255)
                           ccGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGaa                                                                   <  1:5404506/46‑1 (MQ=255)
                           ccGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAAt                                                                 >  1:1408805/1‑48 (MQ=255)
                           ccGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGa                                                          <  1:245251/55‑1 (MQ=255)
                            cGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGt                                         <  1:3055166/71‑1 (MQ=255)
                             ggCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                        >  1:878055/1‑55 (MQ=255)
                             ggCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGtt                                        <  1:2563728/71‑1 (MQ=255)
                              gCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                       <  1:3655638/71‑1 (MQ=255)
                               gatcATAACGGAGGCGTACTGCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGt                                         >  1:1748718/4‑68 (MQ=255)
                                 ccATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATc                                                                <  1:1598296/43‑1 (MQ=255)
                                 ccATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATc                                                                <  1:3048243/43‑1 (MQ=255)
                                 ccATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATc                                                                <  1:642444/43‑1 (MQ=255)
                                 ccATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATc                                                                <  1:2693380/43‑1 (MQ=255)
                                 ccATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                       <  1:3571826/52‑1 (MQ=255)
                                  cATAACGGAGGCGTACTGCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAgaaa                                   >  1:2358490/1‑69 (MQ=255)
                                  cATAACGGAGGCGTACTGCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAgaaa                                   >  1:3207204/1‑69 (MQ=255)
                                  cATAACGGAGGCGTACTGCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAgaa                                    >  1:4646796/1‑69 (MQ=255)
                                    tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                       <  1:2312640/49‑1 (MQ=255)
                                      aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                       >  1:467677/1‑63 (MQ=255)
                                      aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                >  1:2641538/1‑70 (MQ=255)
                                      aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAAt                               <  1:249723/71‑1 (MQ=255)
                                      aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAAt                               <  1:1330378/71‑1 (MQ=255)
                                        cGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAg                                                              <  1:2248939/38‑1 (MQ=255)
                                        cGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGa                                                     <  1:1179556/47‑1 (MQ=255)
                                        cGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAAt                               >  1:3443319/1‑69 (MQ=255)
                                            gCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCgg                         >  1:1004157/1‑71 (MQ=255)
                                             cGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGc                        <  1:4355585/71‑1 (MQ=255)
                                               tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                       <  1:2006210/38‑1 (MQ=255)
                                               tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACg                                                      <  1:3665311/39‑1 (MQ=255)
                                               tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACg                                                      <  1:5605021/39‑1 (MQ=255)
                                               tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACg                                                      <  1:3508099/39‑1 (MQ=255)
                                               tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAg                                           <  1:3006230/50‑1 (MQ=255)
                                               tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAg                                           <  1:1461946/50‑1 (MQ=255)
                                               tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGACATTCACCAg                                           <  1:1082108/50‑1 (MQ=255)
                                                aTTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAgg                                          >  1:349147/1‑50 (MQ=255)
                                                aTTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                              <  1:1196447/62‑1 (MQ=255)
                                                 ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                       <  1:2061554/36‑1 (MQ=255)
                                                 ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                              <  1:223291/45‑1 (MQ=255)
                                                 ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                       >  1:4434754/1‑52 (MQ=255)
                                                 ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                       >  1:2173524/1‑52 (MQ=255)
                                                 ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGAAATTCACAAGGTTGGAGACAATAAGCGGCAgg                     >  1:2493236/1‑70 (MQ=255)
                                                  tCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                       <  1:3576795/35‑1 (MQ=255)
                                                  tCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                   >  1:2091012/1‑39 (MQ=255)
                                                  tCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAg                            >  1:5291843/1‑62 (MQ=255)
                                                   cgcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                              >  1:5292188/1‑43 (MQ=255)
                                                    gcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                   <  1:1122080/37‑1 (MQ=255)
                                                     cgcaaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTgg                 >  1:3387305/4‑71 (MQ=255)
                                                     cgcaaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTgg                 >  1:3393996/4‑71 (MQ=255)
                                                     cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                <  1:1877612/71‑1 (MQ=255)
                                                     cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                >  1:1196194/1‑71 (MQ=255)
                                                     cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                >  1:351590/1‑71 (MQ=255)
                                                     cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                <  1:3513256/71‑1 (MQ=255)
                                                      gAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGtt                                        <  1:4460923/46‑1 (MQ=255)
                                                      gAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGtt                                        <  1:519162/46‑1 (MQ=255)
                                                      gAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg               <  1:3270281/71‑1 (MQ=255)
                                                       caaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACa                                 <  1:5531407/52‑1 (MQ=255)
                                                        aaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCgg                         <  1:3444343/60‑1 (MQ=255)
                                                        aaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCgg                         <  1:11626/60‑1 (MQ=255)
                                                        aaaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTTGAGACAATAAGCGGCAGGCTGGCgg              >  1:2759706/1‑71 (MQ=255)
                                                          aGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCa                                            <  1:5092252/39‑1 (MQ=255)
                                                          aGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTa            >  1:4325507/1‑71 (MQ=255)
                                                           gCCGAGGCCAAATAAATCAGCGGAAACGATATTCACCAGGtt                                        <  1:4819783/42‑1 (MQ=255)
                                                                ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                     >  1:1260618/1‑40 (MQ=255)
                                                                ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                     >  1:3288686/1‑40 (MQ=255)
                                                                ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                     >  1:3708164/1‑40 (MQ=255)
                                                                ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                     >  1:4005825/1‑40 (MQ=255)
                                                                ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                     >  1:5066925/1‑40 (MQ=255)
                                                                ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                     >  1:4091978/1‑40 (MQ=255)
                                                                 gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                       <  1:15690/37‑1 (MQ=255)
                                                                 gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                              <  1:303257/46‑1 (MQ=255)
                                                                  ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                              <  1:5329182/45‑1 (MQ=255)
                                                                  ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt    <  1:61551/71‑1 (MQ=255)
                                                                    aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCa      <  1:254872/67‑1 (MQ=255)
                                                                    aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa  <  1:2147117/71‑1 (MQ=255)
                                                                    |                                                                      
GCGTGGCAACAATCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGA  >  minE/2596027‑2596164

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: