Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 526,368 C→T 18.3% I10I (ATC→ATT ybhL → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE526,3680CT18.3% 229.3 / 32.1 141I10I (ATC→ATTybhLpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (56/59);  new base T (2/24);  total (58/83)
Fisher's exact test for biased strand distribution p-value = 1.01e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 145 total aligned reads displayed.

TATCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTGGATG  >  minE/526301‑526422
                                                                   |                                                      
tATCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCg                                                       <  1:891817/69‑1 (MQ=255)
 aTCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCTCACGTTCTGATTCAATCGt                                                      <  1:1008918/69‑1 (MQ=255)
 aTCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGt                                                      <  1:1587201/69‑1 (MQ=255)
 aTCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGt                                                      <  1:1939621/69‑1 (MQ=255)
 aTCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGt                                                      <  1:377210/69‑1 (MQ=255)
 aTCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGt                                                      <  1:1012583/69‑1 (MQ=255)
  tCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGt                                                      >  1:1782234/1‑68 (MQ=255)
    cTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACa                                                   >  1:3212815/1‑69 (MQ=255)
    cTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACa                                                   >  1:1002545/1‑69 (MQ=255)
    cTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACa                                                   >  1:1894578/1‑69 (MQ=255)
       aaGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAAc                                                 >  1:659976/1‑68 (MQ=255)
            atTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCgg                                            >  1:406092/1‑68 (MQ=255)
             ttCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggc                                          >  1:606969/1‑69 (MQ=255)
             ttCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggc                                          >  1:3037776/1‑69 (MQ=255)
             ttCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggc                                          >  1:2232551/1‑69 (MQ=255)
             ttCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggc                                          >  1:2218868/1‑69 (MQ=255)
             ttCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggc                                          >  1:2754206/1‑69 (MQ=255)
             ttCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggc                                          >  1:1317530/1‑69 (MQ=255)
             ttCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggc                                          >  1:1361609/1‑69 (MQ=255)
               ccTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggctg                                        >  1:2965391/1‑69 (MQ=255)
               ccTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggctg                                        >  1:2133421/1‑69 (MQ=255)
               ccTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggctg                                        >  1:31191/1‑69 (MQ=255)
               ccTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGggctg                                        >  1:2113875/1‑69 (MQ=255)
                    tatCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCtt                                    >  1:2474156/1‑68 (MQ=255)
                    tatCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCtt                                    >  1:1288089/1‑68 (MQ=255)
                        ttCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCaaa                               >  1:2298287/1‑69 (MQ=255)
                        ttCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCaaa                               >  1:2800902/1‑69 (MQ=255)
                        ttCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCaaa                               >  1:662938/1‑69 (MQ=255)
                        ttCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCaaa                               >  1:2882197/1‑69 (MQ=255)
                        ttCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCaa                                <  1:918791/68‑1 (MQ=255)
                        ttCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCaa                                <  1:2194692/68‑1 (MQ=255)
                        ttCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCaa                                <  1:36275/68‑1 (MQ=255)
                        ttCAGGAGAAAGTCATGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCaa                                <  1:2851144/68‑1 (MQ=255)
                        ttCAGGAGAAAGTCATGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCaa                                <  1:267476/68‑1 (MQ=255)
                        ttCAGGAGAAAGTCATGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCaa                                <  1:3043035/68‑1 (MQ=255)
                        ttCAGGAGAAAGTCATGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCaa                                <  1:2501637/68‑1 (MQ=255)
                        ttCAGGAGAAAGTCATGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCaa                                <  1:1168550/68‑1 (MQ=255)
                        ttCAGGAGAAAGTCATGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCaa                                <  1:1341706/68‑1 (MQ=255)
                         tCAGGAGAAAGTCATGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCaa                                <  1:3005804/67‑1 (MQ=255)
                         tCAGGAGAAAGTCATGGACAGATTCCCACGTTCAGATTCACATGTACAACCCCGGGCTGGCTTTCaa                                <  1:2788130/67‑1 (MQ=255)
                            ggAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCaa                                <  1:620441/64‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:1859307/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:3054621/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:76480/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:976940/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2870895/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:3155818/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:3179561/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2739517/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:617232/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:3259301/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:33457/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:1260855/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:988782/68‑1 (MQ=255)
                             gagaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2332476/68‑1 (MQ=255)
                              agaTCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2061198/67‑1 (MQ=255)
                                aaagTCATGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCaa                                <  1:112854/57‑1 (MQ=255)
                                  cGTCATGGACAGATTCGCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           >  1:2604962/1‑63 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTAc                                                    >  1:2548494/1‑38 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTAc                                                    >  1:2469142/1‑38 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAATCTta                           >  1:549683/1‑63 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:669159/63‑1 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           >  1:3193694/1‑63 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2922812/63‑1 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2947520/63‑1 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:142404/63‑1 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           >  1:1932659/1‑63 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           >  1:1319142/1‑63 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2680626/63‑1 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           >  1:261494/1‑63 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:1350681/63‑1 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2975541/63‑1 (MQ=255)
                                  cGTCATGGACAGATTCCCACGTTCTGATTAATCGTAc                                                    >  1:3275551/1‑37 (MQ=38)
                                      aTGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCaa                 <  1:517456/69‑1 (MQ=255)
                                      aTGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCaa                 <  1:811006/69‑1 (MQ=255)
                                      aTGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCaa                 <  1:2013531/69‑1 (MQ=255)
                                      aTGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCaa                 <  1:584273/69‑1 (MQ=255)
                                      aTGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCaa                 <  1:567071/69‑1 (MQ=255)
                                      aTGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCaa                 <  1:336155/69‑1 (MQ=255)
                                       ttgACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCaa                 <  1:252803/66‑1 (MQ=255)
                                       tGGACAGATTCCCACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCaa                 <  1:2280362/68‑1 (MQ=255)
                                          aCAGATTCCGACGTTCTGATTCAATTGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCGCAAGTCt             >  1:1122886/1‑69 (MQ=255)
                                           cAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTta                           <  1:2578411/54‑1 (MQ=255)
                                            aGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTAt           >  1:458216/1‑69 (MQ=255)
                                             gATTCCCACGTTCTGATTCAATCTTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATg          >  1:818813/1‑69 (MQ=255)
                                             gATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATg          >  1:388711/1‑69 (MQ=255)
                                             gATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATg          >  1:150780/1‑69 (MQ=255)
                                             gATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATg          >  1:1584778/1‑69 (MQ=255)
                                             gATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATg          >  1:2100281/1‑69 (MQ=255)
                                             gATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATg          >  1:2089501/1‑69 (MQ=255)
                                             gATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATg          >  1:2812402/1‑69 (MQ=255)
                                             gATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATg          >  1:1927601/1‑69 (MQ=255)
                                                 cccACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTg      <  1:1686671/69‑1 (MQ=255)
                                                 cccACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTg      <  1:776621/69‑1 (MQ=255)
                                                 cccACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTg      <  1:3164949/69‑1 (MQ=255)
                                                 cccACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTg      <  1:2525913/69‑1 (MQ=255)
                                                 cccACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTg      <  1:1190679/69‑1 (MQ=255)
                                                                   cGTACAGCCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTGGa    >  1:437556/1‑53 (MQ=255)
                                                                   cGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTGGATg  >  1:2821467/1‑55 (MQ=255)
                                                                   cGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTGGATg  >  1:2333295/1‑55 (MQ=255)
                                                                   |                                                      
TATCCTTAAGATATTCCTTATATCTTCAGGAGATCGTCATGGACAGATTCCCACGTTCTGATTCAATCGTACAACCCCGGGCTGGCTTGCAAACTTATATGGCTCAAGTCTATGGCTGGATG  >  minE/526301‑526422

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: