Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,232 T→G 12.5% I22L (ATC→CTC) ‡ arsC ← arsenate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,2320TG12.5% 192.7 / 25.7 88I22L (ATC→CTC) ‡arsCarsenate reductase
Reads supporting (aligned to +/- strand):  ref base T (44/33);  new base G (5/6);  total (49/39)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTA  >  minE/2595163‑2595302
                                                                     |                                                                      
ccAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGat                                                                        >  1:1443439/1‑70 (MQ=255)
    ttCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtct                                                                   >  1:2661998/1‑71 (MQ=255)
          gcgcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                                             >  1:1073230/1‑71 (MQ=255)
                tGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAc                                                       >  1:3429468/1‑71 (MQ=255)
                  gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCa                                                                >  1:1874774/1‑60 (MQ=255)
                  gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                                     <  1:3148513/71‑1 (MQ=255)
                  gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                                     <  1:2682631/71‑1 (MQ=255)
                  gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGATGTTGCGGATCATCTCCAGCGTATTAcgc                                                     <  1:572767/71‑1 (MQ=255)
                   cGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgcg                                                    <  1:3269313/71‑1 (MQ=255)
                      aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                        >  1:833701/1‑64 (MQ=255)
                      aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                                 <  1:879926/71‑1 (MQ=255)
                      aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                                 >  1:3149622/1‑71 (MQ=255)
                      aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                                 >  1:3305580/1‑71 (MQ=255)
                        tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                               >  1:2816590/1‑71 (MQ=255)
                        tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                               >  1:3519560/1‑71 (MQ=255)
                        tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                               >  1:2440572/1‑71 (MQ=255)
                         ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                                              >  1:1903062/1‑71 (MQ=255)
                         ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                                              >  1:1793698/1‑71 (MQ=255)
                          tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                                              <  1:210215/67‑1 (MQ=255)
                          tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                                              >  1:2230615/4‑70 (MQ=255)
                          tCCAGATAATGGATAATAGTCGGTTCTGTGCCGTTGTTGCGGATCAtctc                                                                  <  1:124172/50‑1 (MQ=255)
                          tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                                              <  1:1160763/70‑1 (MQ=255)
                          tCCAGATAATGGATAATAGTCGGTTCTGTGACGCTGTTGCGGATCAtctc                                                                  <  1:2477670/50‑1 (MQ=255)
                           ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtctc                                                                  >  1:3559069/1‑49 (MQ=255)
                           ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                                               >  1:2373334/1‑52 (MQ=255)
                           ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                             <  1:1578844/70‑1 (MQ=255)
                           ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                             <  1:3564941/70‑1 (MQ=255)
                            gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGccac                                           <  1:1365175/70‑3 (MQ=255)
                            cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                           <  1:3867493/71‑1 (MQ=255)
                            cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                           <  1:545457/71‑1 (MQ=255)
                             aGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCcaca                                          <  1:3064865/71‑1 (MQ=255)
                             aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                        >  1:3419162/1‑57 (MQ=255)
                             aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                        >  1:3889924/1‑57 (MQ=255)
                             aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                            >  1:325666/1‑69 (MQ=255)
                             aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                            >  1:2546880/1‑69 (MQ=255)
                               aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAgg                                        >  1:31737/1‑71 (MQ=255)
                                tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                       >  1:2706900/1‑71 (MQ=255)
                                tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                       <  1:2711013/71‑1 (MQ=255)
                                 aaTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCa                                          >  1:2572014/1‑67 (MQ=255)
                                  aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGGGACGTGCCGCAGGcc                                      <  1:1000847/70‑1 (MQ=255)
                                  aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                            <  1:3754353/64‑1 (MQ=255)
                                   tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAt                                                          <  1:3019104/49‑1 (MQ=255)
                                   tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCgg                                    >  1:693838/1‑71 (MQ=255)
                                   tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCgg                                    >  1:696829/1‑71 (MQ=255)
                                    ggATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAt                                                          >  1:1299870/1‑48 (MQ=255)
                                      ataataGTCGGTTCTGTGCCGCTGTTGCTGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                 <  1:1677810/71‑1 (MQ=255)
                                      ataataGTCGGTTCTGTGCCGCTGTTGCGGatcatc                                                                    <  1:2705328/36‑1 (MQ=255)
                                      ataataGTCGGTTCTGTGCCGCTGTTGCGGatcat                                                                     >  1:1218715/1‑35 (MQ=255)
                                      ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                 <  1:817088/71‑1 (MQ=255)
                                      ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                 <  1:1060286/71‑1 (MQ=255)
                                      ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                 <  1:1867950/71‑1 (MQ=255)
                                      ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGt                                  >  1:3857410/1‑70 (MQ=255)
                                       taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                                   >  1:2391373/1‑52 (MQ=255)
                                       taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                                   >  1:3084155/1‑52 (MQ=255)
                                        aataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                               >  1:1556364/1‑71 (MQ=255)
                                         attgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                               <  1:2587251/67‑1 (MQ=255)
                                         ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                              <  1:3801900/71‑1 (MQ=255)
                                         ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                              <  1:676431/71‑1 (MQ=255)
                                         ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                              <  1:57239/71‑1 (MQ=255)
                                          ttgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGa                             <  1:672465/69‑1 (MQ=255)
                                          taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                        >  1:3483684/1‑44 (MQ=255)
                                          taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                        >  1:1489284/1‑44 (MQ=255)
                                          taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                           >  1:2458016/1‑57 (MQ=255)
                                          taGTCGGTTCTGTGCAGCTGTTGCGGATCATCTCCAGCGTATTAcg                                                      >  1:1114239/1‑46 (MQ=255)
                                           tgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgtg                              >  1:3721245/2‑69 (MQ=255)
                                            gTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgtg                              >  1:217582/1‑68 (MQ=255)
                                             tCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATa                           <  1:1639213/70‑1 (MQ=255)
                                              cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                          >  1:1503227/1‑70 (MQ=255)
                                               ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa                         <  1:3376541/70‑1 (MQ=255)
                                               ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAAt                        >  1:166687/1‑71 (MQ=255)
                                               ggTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGa                                                   >  1:1047466/1‑44 (MQ=255)
                                               ggTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGAt                        >  1:3355037/1‑71 (MQ=255)
                                                   cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                     >  1:322149/1‑70 (MQ=255)
                                                    tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCggg                                   <  1:1398041/55‑1 (MQ=255)
                                                      tgCCGCTGTTGCGTATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                              <  1:1139418/58‑1 (MQ=255)
                                                      tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                 >  1:2599270/1‑71 (MQ=255)
                                                        ccGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                   >  1:3530304/1‑67 (MQ=255)
                                                            tGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt            >  1:363305/1‑70 (MQ=255)
                                                              ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                       <  1:29745/57‑1 (MQ=255)
                                                              ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTT‑CTCat         <  1:2393005/70‑1 (MQ=255)
                                                               tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt            >  1:1855616/1‑67 (MQ=255)
                                                               tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt            >  1:2881427/1‑67 (MQ=255)
                                                               tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata        >  1:497094/1‑71 (MQ=255)
                                                               tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata        >  1:2731512/1‑71 (MQ=255)
                                                               tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata        >  1:2579112/1‑71 (MQ=255)
                                                               tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata        >  1:1013104/1‑71 (MQ=255)
                                                                 cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                          <  1:1213714/51‑1 (MQ=255)
                                                                  ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                          <  1:3132729/50‑1 (MQ=255)
                                                                  ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                   <  1:2820987/57‑1 (MQ=255)
                                                                  ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCat         <  1:3618691/67‑1 (MQ=255)
                                                                   gATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg    <  1:2906178/71‑1 (MQ=255)
                                                                    atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg    >  1:2867014/1‑70 (MQ=255)
                                                                    atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg    >  1:1012718/1‑70 (MQ=255)
                                                                    atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt   <  1:3667590/71‑1 (MQ=255)
                                                                    atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt   <  1:1984321/71‑1 (MQ=255)
                                                                     tcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa  >  1:1931336/1‑71 (MQ=255)
                                                                     |                                                                      
CCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTA  >  minE/2595163‑2595302

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: ‑