Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,215 T→G 10.2% T27T (ACA→ACC) ‡ arsC ← arsenate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,2150TG10.2% 232.6 / 21.3 98T27T (ACA→ACC) ‡arsCarsenate reductase
Reads supporting (aligned to +/- strand):  ref base T (41/47);  new base G (3/7);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 111 total aligned reads displayed.

TATCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGT  >  minE/2595145‑2595283
                                                                      |                                                                    
tatCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCt                                                                      <  1:1297752/71‑1 (MQ=255)
  tCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCtgt                                                                    <  1:2903383/71‑1 (MQ=255)
  tCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCtgt                                                                    <  1:4015992/71‑1 (MQ=255)
        aTGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCt                                                              <  1:1599215/71‑1 (MQ=255)
        aTGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTGCTGTGCCGCt                                                              <  1:5434576/71‑1 (MQ=255)
          gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTg                                                             >  1:3059578/1‑70 (MQ=255)
          gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTg                                                             >  1:3610266/1‑70 (MQ=255)
          gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTg                                                             >  1:4004157/1‑70 (MQ=255)
          gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGt                                                            >  1:2506547/1‑71 (MQ=255)
          gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGt                                                            >  1:1113835/1‑71 (MQ=255)
           aGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGtt                                                           >  1:3375841/1‑71 (MQ=255)
               tGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCg                                                        >  1:420386/1‑70 (MQ=255)
               tGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCg                                                        >  1:3540696/1‑70 (MQ=255)
                     gTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGatcat                                                  >  1:5034946/1‑70 (MQ=255)
                      ttCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGc                                                         <  1:3169537/62‑1 (MQ=255)
                      ttCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtct                                                <  1:4188503/71‑1 (MQ=255)
                      ttCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtct                                                <  1:1235344/71‑1 (MQ=255)
                        cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCt                                                                      >  1:1001450/1‑47 (MQ=255)
                        cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCgg                                                       >  1:5314449/1‑62 (MQ=255)
                        cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtctc                                               >  1:3720909/1‑70 (MQ=255)
                        cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtctc                                               >  1:156286/1‑70 (MQ=255)
                         aTCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCgg                                                       >  1:398630/1‑61 (MQ=255)
                         aTCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTcc                                              <  1:677122/70‑1 (MQ=255)
                           cgcgcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGc                                                         <  1:3772907/57‑1 (MQ=255)
                            gcgcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                          <  1:1662411/71‑1 (MQ=255)
                              gcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                        <  1:4507579/71‑1 (MQ=255)
                               cgTTGGCGGAGTTTCCAGATAATGGATAATAGTCTGTTCTGTGCCGCTGTTGCGGatc                                                    <  1:5663473/58‑1 (MQ=255)
                               cgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGatc                                                    <  1:1422093/58‑1 (MQ=255)
                                 ttGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCtgt                                                                    <  1:727457/40‑1 (MQ=255)
                                 ttGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCgg                                                       <  1:3189027/53‑1 (MQ=255)
                                  tGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCtg                                                                     >  1:1582131/1‑38 (MQ=255)
                                   ggCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAc                                    <  1:1590286/70‑1 (MQ=255)
                                    gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCg                                                        <  1:1617702/49‑1 (MQ=255)
                                    gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                  <  1:2444658/71‑1 (MQ=255)
                                    gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcg                                   <  1:772138/70‑1 (MQ=255)
                                      ggAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                <  1:3161115/71‑1 (MQ=255)
                                      ggAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                <  1:26413/71‑1 (MQ=255)
                                        aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                              >  1:214205/1‑71 (MQ=255)
                                         gTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGt                             <  1:3291164/71‑1 (MQ=255)
                                          tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTg                                                             <  1:5458828/38‑1 (MQ=255)
                                          tttCCAGATAATAGATAATGGTCGGTTCTGTACCGCTGTTGCGGATCATCTCCAGCGTATTACGTGACGt                             <  1:1088262/70‑1 (MQ=255)
                                           ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                          <  1:4389756/56‑1 (MQ=255)
                                           ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGAc                               <  1:32204/67‑1 (MQ=255)
                                            tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                           >  1:3430818/4‑70 (MQ=255)
                                            tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGtt                                                           >  1:3198932/1‑38 (MQ=255)
                                            tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGat                                                     >  1:3112716/1‑44 (MQ=255)
                                            tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGc                                           >  1:3515282/1‑54 (MQ=255)
                                            tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGc                                           >  1:5279192/1‑54 (MQ=255)
                                            tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                          >  1:532932/1‑71 (MQ=255)
                                             ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                         >  1:2515271/1‑71 (MQ=255)
                                              gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGccac                        <  1:5354453/70‑3 (MQ=255)
                                              gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGcca                         <  1:4622529/69‑2 (MQ=255)
                                              gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGcca                         <  1:3918693/69‑2 (MQ=255)
                                              gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGcca                         <  1:1903951/69‑2 (MQ=255)
                                              cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                  <  1:2921762/61‑1 (MQ=255)
                                              cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                  <  1:2770303/61‑1 (MQ=255)
                                                 aTAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCg                                          <  1:4125509/50‑1 (MQ=255)
                                                 aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                          <  1:2008825/50‑1 (MQ=255)
                                                 aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                          <  1:1369154/50‑1 (MQ=255)
                                                    aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCa                       >  1:4934723/1‑66 (MQ=255)
                                                    aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCa                       >  1:3825895/1‑66 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                            <  1:3649309/42‑1 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                            <  1:2043411/42‑1 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                            <  1:1751469/42‑1 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                            <  1:589468/42‑1 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                            <  1:559705/42‑1 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                            <  1:3691019/42‑1 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                        >  1:631560/1‑46 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                        >  1:1638362/1‑46 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGt                             >  1:446952/1‑57 (MQ=255)
                                                       gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGt               >  1:2619882/1‑71 (MQ=255)
                                                        ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt              >  1:2615485/1‑71 (MQ=255)
                                                         taataGTCGGTTCTGTGCCCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt            <  1:552968/71‑1 (MQ=255)
                                                           attgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt            >  1:765989/4‑70 (MQ=255)
                                                           ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt            <  1:1234973/70‑1 (MQ=255)
                                                            taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                     >  1:2146850/1‑44 (MQ=255)
                                                            taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                     >  1:2270696/1‑44 (MQ=255)
                                                            taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                          <  1:5451187/55‑1 (MQ=255)
                                                            taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                          <  1:3382434/55‑1 (MQ=255)
                                                              gTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAtt                                      >  1:5068052/1‑41 (MQ=255)
                                                              gTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCCTAtt                                      >  1:4003727/1‑41 (MQ=255)
                                                              gTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGa          >  1:1025140/1‑69 (MQ=255)
                                                               tCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa       <  1:930599/71‑1 (MQ=255)
                                                               tCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa       <  1:1132343/71‑1 (MQ=255)
                                                               tCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATa        <  1:3370188/70‑1 (MQ=255)
                                                                cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                        >  1:4416215/1‑37 (MQ=255)
                                                                cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa      >  1:5471374/1‑71 (MQ=255)
                                                                 ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAAt     >  1:5296498/1‑71 (MQ=255)
                                                                 ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAAt     <  1:1609557/71‑1 (MQ=255)
                                                                  gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                         <  1:5081422/50‑1 (MQ=255)
                                                                  gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                         <  1:3248690/50‑1 (MQ=255)
                                                                  gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                         <  1:410810/50‑1 (MQ=255)
                                                                   ttCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg   >  1:5562081/1‑71 (MQ=255)
                                                                   ttCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg   <  1:5046617/71‑1 (MQ=255)
                                                                     cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt  >  1:2132299/1‑70 (MQ=255)
                                                                     cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt  <  1:1929783/70‑1 (MQ=255)
                                                                     cTGTGCCGCTGTTGCAGATCATCTCCAGCGTATTACGCGACGTGCCGCa                       >  1:3609444/1‑49 (MQ=255)
                                                                      tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                <  1:4425046/39‑1 (MQ=255)
                                                                      tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCg                  <  1:5660910/53‑1 (MQ=255)
                                                                      tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCg                  <  1:2059901/53‑1 (MQ=255)
                                                                      |                                                                    
TATCGGCAATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGT  >  minE/2595145‑2595283

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