Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,445,467 G→A 52.6% T638T (ACC→ACT rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,445,4670GA52.6% 1.1 / 66.5 76T638T (ACC→ACTrcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (8/28);  new base A (20/20);  total (28/48)
Fisher's exact test for biased strand distribution p-value = 1.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TCTGAATCTACCGAGATATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGC  >  minE/1445399‑1445527
                                                                    |                                                            
tCTGAATCTACCGAGATATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCg                                                              >  1:1912097/1‑69 (MQ=255)
tCTGAATCTACCGAGATATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCg                                                              >  1:3101644/1‑69 (MQ=255)
tCTGAATCTACCGAGATATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCg                                                              >  1:3323047/1‑69 (MQ=255)
             agaTGTCGCCGTCCATCATTCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:117279/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:1981432/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:2099521/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:1770541/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:2798629/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:3239913/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:1493801/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:3370533/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:129129/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:3429381/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:1232976/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:122930/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:433293/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:1167173/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:1206609/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:1214378/60‑1 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:947967/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:3396730/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:2132491/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:2184313/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:2510476/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:2572795/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:689303/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:2988392/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:3031890/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:447936/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:3258170/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:1658471/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:1647124/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:1940991/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:155698/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:176970/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaagt                                                 >  1:1256618/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaag                                                  >  1:12253/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaag                                                  >  1:674549/1‑67 (MQ=255)
             agaTGTCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGaa                                                   >  1:652851/1‑67 (MQ=255)
             agaTATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAAt                                                 <  1:1229136/69‑1 (MQ=255)
             agaTATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAAt                                                 <  1:2578417/69‑1 (MQ=255)
             agaTATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAAt                                                 <  1:2486571/69‑1 (MQ=255)
             agaTATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAAt                                                 <  1:2063475/69‑1 (MQ=255)
             agaTATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAAt                                                 <  1:1998596/69‑1 (MQ=255)
                  tCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGc                                            >  1:2224455/1‑69 (MQ=255)
                    gCCGTCCATCATGCTTATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:1759294/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:1029195/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:487795/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:749939/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:140988/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:756594/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:3223888/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:3209021/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:3010945/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:2154455/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:2983609/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:2337538/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:2351374/69‑1 (MQ=255)
                    gCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:1820746/69‑1 (MQ=255)
                    gCCGTACATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTg                                          <  1:440221/69‑1 (MQ=255)
                         ccATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAAt                                                 <  1:1866771/57‑1 (MQ=255)
                          catcatGCTTATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAAt                                                 <  1:2901253/56‑1 (MQ=255)
                          catcatGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGcc                                    >  1:2677972/1‑69 (MQ=255)
                          catcatGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGcc                                    >  1:1668550/1‑69 (MQ=255)
                          catcatGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCAGTcc                                                          <  1:1522312/47‑1 (MQ=255)
                                      cAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGt                                         >  1:3371530/1‑52 (MQ=255)
                                      cAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGAAGTTACGCTGTACGcc                                    <  1:1505801/57‑1 (MQ=39)
                                      cAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGAAGTTACGCTGTACGcc                                    <  1:1228291/57‑1 (MQ=39)
                                      cAGTTTTTCACAAATCGCCAGACCCAGACCAGTCCCCTGGAAGTTACGCTGTACGcc                                    <  1:1829841/57‑1 (MQ=39)
                                               acaAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAgg                <  1:3180406/68‑1 (MQ=255)
                                               acaAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAgg                <  1:807758/68‑1 (MQ=255)
                                               acaAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAgg                <  1:1863530/68‑1 (MQ=255)
                                               acaAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAgg                <  1:2150730/68‑1 (MQ=255)
                                               acaAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAgg                <  1:2343380/68‑1 (MQ=255)
                                               acaAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAgg                <  1:333780/68‑1 (MQ=255)
                                                            cccagaccGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAgcgc  >  1:656123/1‑69 (MQ=255)
                                                            cccagaccGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAgcgc  >  1:593817/1‑69 (MQ=255)
                                                            cccagaccGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAgcgc  >  1:3427397/1‑69 (MQ=255)
                                                            cccagaccGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAgcgc  >  1:1284972/1‑69 (MQ=255)
                                                                    |                                                            
TCTGAATCTACCGAGATATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGC  >  minE/1445399‑1445527

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: