Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,886,458 A→G 33.3% T640T (ACA→ACG barA → hybrid sensory histidine kinase, in two‑component regulatory system with UvrY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,886,4580AG33.3% 73.1 / 52.3 96T640T (ACA→ACGbarAhybrid sensory histidine kinase, in two‑component regulatory system with UvrY
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base G (10/22);  total (40/56)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGAT  >  minE/1886393‑1886524
                                                                 |                                                                  
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:2867342/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:3164234/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:1950162/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:3214731/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:1635109/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:1385188/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:3317120/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:460596/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:2277450/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:2448624/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:1175396/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:752975/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:678896/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:2356675/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGc                                                                 >  1:723286/1‑69 (MQ=255)
gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGATGAAACCATTAACACCTACACGc                                                                 >  1:209191/1‑69 (MQ=255)
 tCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGcc                                                                <  1:2382651/69‑1 (MQ=255)
      gCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTtg                                                           <  1:2239520/69‑1 (MQ=255)
            aaaCTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTtgcctg                                                       <  1:3333412/67‑1 (MQ=255)
            aaaCTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTtgcctg                                                       <  1:2891397/67‑1 (MQ=255)
            aaaCTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTtgcctg                                                       <  1:420822/67‑1 (MQ=255)
            aaaCTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTtgcctg                                                       <  1:279035/67‑1 (MQ=255)
            aaaCTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTtgcctg                                                       <  1:1099806/67‑1 (MQ=255)
            aaaCTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTtgcctg                                                       <  1:1213841/67‑1 (MQ=255)
            aaaCTCAAGCAAGATGGTATCGGCGCGTGTCCGCTGAAACCATTAACACCTACACGCCTGTtgcctg                                                       <  1:2095530/67‑1 (MQ=255)
                           ggTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:2967943/69‑1 (MQ=255)
                           ggTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:1903604/69‑1 (MQ=255)
                           ggTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:3144744/69‑1 (MQ=255)
                           ggTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:1016376/69‑1 (MQ=255)
                           ggTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:1986022/69‑1 (MQ=255)
                           ggTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:3249181/69‑1 (MQ=255)
                            gTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTtgcctg                                                       <  1:1239609/51‑1 (MQ=255)
                            gTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:3312959/68‑1 (MQ=255)
                            gTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:1386143/68‑1 (MQ=255)
                              aTCGGCGCGTGTCTGCTGAAACCATTAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGtcat                                   <  1:3096685/69‑1 (MQ=255)
                              aTCGGCGCGTGTCTGCTGAAACCATTAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGtcat                                   <  1:3124107/69‑1 (MQ=255)
                              aTCGGCGCGTGTCTGCTGAAACCATTAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGtcat                                   <  1:1075248/69‑1 (MQ=255)
                              aTCGGCGCGTGTCTGCTGAAACCATTAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGtcat                                   <  1:775360/69‑1 (MQ=255)
                              aTCGGCGCGTGTCTGCTGAAACCATTAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGtcat                                   <  1:867339/69‑1 (MQ=255)
                              aTCGGCGCGTGTCTGCTGAAACCATTAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGtcat                                   <  1:239624/69‑1 (MQ=255)
                                    gcgTGTCTGCTGAAACCATTAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGtcat                                   <  1:1958330/63‑1 (MQ=255)
                                          ctgctgAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGt                                      <  1:37205/54‑1 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  <  1:2538935/69‑1 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  <  1:2947058/69‑1 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  <  1:317822/69‑1 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  <  1:3176316/69‑1 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  >  1:2416872/1‑69 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  >  1:1581125/1‑69 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  >  1:161824/1‑69 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  <  1:2008658/69‑1 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  <  1:1187058/69‑1 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  <  1:2212349/69‑1 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  >  1:2258933/1‑69 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  >  1:10171/1‑69 (MQ=255)
                                               gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCt                  <  1:1045466/69‑1 (MQ=255)
                                                   ccATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTg              <  1:2551079/69‑1 (MQ=255)
                                                   ccATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTg              <  1:1982318/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:1240590/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:337763/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:352694/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:3284754/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:1280492/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:2531791/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:1895612/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:10044/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:602692/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:922288/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:1065199/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:2863490/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:745618/69‑1 (MQ=255)
                                                      ttAACACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTATTGCCt           <  1:440434/69‑1 (MQ=255)
                                                      ttAAAACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTGTCATCACAAACAAAACACGCTTTTGCCt           <  1:689859/69‑1 (MQ=255)
                                                        aaCACGTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:2898071/1‑39 (MQ=37)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:512951/1‑39 (MQ=38)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:997750/1‑39 (MQ=38)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:357094/1‑39 (MQ=38)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:1360380/1‑39 (MQ=38)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:3280675/1‑39 (MQ=38)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:2561926/1‑39 (MQ=38)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:2016002/1‑39 (MQ=38)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:266751/1‑39 (MQ=38)
                                                        aaCACCTACGCGCCTGTTGCCTGCCCTGACAGAATTTTg                                       >  1:202002/1‑39 (MQ=38)
                                                         acacCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTg              <  1:1061388/63‑1 (MQ=255)
                                                            ccTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAAcc     <  1:1955496/69‑1 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:2181997/1‑50 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:884061/1‑50 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:826598/1‑50 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:2229162/1‑50 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:722600/1‑50 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:133643/1‑50 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:3140972/1‑50 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:1299557/1‑50 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTa        <  1:660037/64‑1 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTa        <  1:483617/64‑1 (MQ=255)
                                                              tACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTa        <  1:2552570/64‑1 (MQ=255)
                                                              tACACGACTGTTGCCTGCCCTGACGGAATTTTGTCATCACaaacaaaaca                      >  1:439474/1‑50 (MQ=255)
                                                                cacGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGAt  >  1:2209107/1‑68 (MQ=255)
                                                                cacGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGAt  >  1:2383965/1‑68 (MQ=255)
                                                                cacGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGAt  >  1:2488722/1‑68 (MQ=255)
                                                                 |                                                                  
GTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGAT  >  minE/1886393‑1886524

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: