Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 484,413 C→T 12.4% R73R (CGC→CGT)  tolB → periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE484,4130CT12.4% 263.1 / 22.7 97R73R (CGC→CGT) tolBperiplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (51/34);  new base T (9/3);  total (60/37)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 101 total aligned reads displayed.

CACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAT  >  minE/484343‑484481
                                                                      |                                                                    
cACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATcgc                                                                      <  1:3783020/71‑1 (MQ=255)
     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCt                                                                   >  1:2156131/1‑69 (MQ=255)
     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:973481/1‑71 (MQ=255)
     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:2346452/1‑71 (MQ=255)
     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:3990787/1‑71 (MQ=255)
     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:1441452/1‑71 (MQ=255)
     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:1293722/1‑71 (MQ=255)
      aaGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTc                                                                  <  1:2251711/69‑1 (MQ=255)
      aaGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTc                                                                  <  1:2526762/69‑1 (MQ=255)
      aaGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTc                                                                  <  1:183243/69‑1 (MQ=255)
          tatTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGc                                                            >  1:5032444/1‑71 (MQ=255)
            ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:5024136/71‑1 (MQ=255)
            ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:3934714/71‑1 (MQ=255)
            ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:1911562/71‑1 (MQ=255)
            ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:3968939/71‑1 (MQ=255)
             tGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:415750/70‑1 (MQ=255)
                  gcATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCAcagcag                                                    <  1:2230236/71‑1 (MQ=255)
                   cATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACagcagc                                                   <  1:3470395/71‑1 (MQ=255)
                    aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACagcagc                                                   <  1:4497911/70‑1 (MQ=255)
                    aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  <  1:3116447/71‑1 (MQ=255)
                    aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  <  1:304255/71‑1 (MQ=255)
                        ttGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGt                                              <  1:134081/71‑1 (MQ=255)
                              ctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTgc                                         >  1:3200793/1‑70 (MQ=255)
                              ctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTgc                                         <  1:1512489/70‑1 (MQ=255)
                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:1405030/1‑69 (MQ=255)
                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:2568564/1‑69 (MQ=255)
                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:1248139/1‑69 (MQ=255)
                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:417388/1‑69 (MQ=255)
                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:3322310/1‑69 (MQ=255)
                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:5124175/1‑69 (MQ=255)
                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:3794876/69‑1 (MQ=255)
                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:3697408/69‑1 (MQ=255)
                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:4171059/69‑1 (MQ=255)
                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 >  1:1386852/1‑71 (MQ=255)
                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 >  1:3111809/1‑71 (MQ=255)
                                      cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4077991/57‑1 (MQ=255)
                                      cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4765968/57‑1 (MQ=255)
                                      cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4174140/57‑1 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCAcagca                                                     >  1:2256892/1‑48 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:3124706/1‑57 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:3605645/1‑57 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:1793924/1‑57 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:2316276/1‑57 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:4326772/1‑57 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:188709/1‑57 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:2189797/1‑57 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              <  1:3891415/71‑1 (MQ=255)
                                        tAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              <  1:3475548/71‑1 (MQ=255)
                                         aaCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTg                                                             >  1:2156493/1‑39 (MQ=255)
                                         aaCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              >  1:3781454/1‑70 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:474664/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:5114070/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:767597/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3999399/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:1080455/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:1439391/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:2354194/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:248188/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3026030/1‑43 (MQ=255)
                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3312105/1‑43 (MQ=255)
                                            aGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:4912169/71‑1 (MQ=255)
                                                gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:3290992/67‑1 (MQ=255)
                                                gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:4653797/67‑1 (MQ=255)
                                                gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:5161206/67‑1 (MQ=255)
                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:531551/54‑1 (MQ=255)
                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:4182551/54‑1 (MQ=255)
                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:5189153/54‑1 (MQ=255)
                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACaa                            <  1:4624919/62‑1 (MQ=255)
                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACaa                            <  1:172070/62‑1 (MQ=255)
                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCa                   >  1:217982/1‑71 (MQ=255)
                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCa                   >  1:572602/1‑71 (MQ=255)
                                                     tttAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAc                           >  1:3086278/1‑61 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  >  1:3363271/1‑35 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:744349/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:627231/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:706034/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1408754/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:76171/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:537950/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1571011/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1582753/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:105250/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1794182/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1921216/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:2102796/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:398656/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:3733995/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:378018/1‑71 (MQ=255)
                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:3799750/1‑71 (MQ=255)
                                                          tCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCg            >  1:536670/1‑71 (MQ=255)
                                                            cGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGc           <  1:1703316/70‑1 (MQ=255)
                                                              ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 <  1:1191845/46‑1 (MQ=255)
                                                              ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:3202075/71‑1 (MQ=255)
                                                              ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:3689644/71‑1 (MQ=255)
                                                              ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:58573/71‑1 (MQ=255)
                                                                aGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTgg      >  1:3463048/1‑71 (MQ=255)
                                                                aGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTgg      >  1:3944476/1‑71 (MQ=255)
                                                                    cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAt  >  1:2375499/1‑71 (MQ=255)
                                                                    cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAt  >  1:802516/1‑71 (MQ=255)
                                                                    cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGAATGGTCCGCACTGGGCAt  >  1:1718268/1‑71 (MQ=255)
                                                                      |                                                                    
CACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAT  >  minE/484343‑484481

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