Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,916,386 G→A 100% E892K (GAA→AAA)  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
*W3110S.gb2,916,3860GA77.8% 49.4 / 12.0 27E892K (GAA→AAA) barAhybrid sensory histidine kinase, in two‑component regulatory system with UvrY
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (6/15);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

CGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGC  >  W3110S.gb/2916332‑2916450
                                                      |                                                                
cGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGATTTGGAAcc                                                  >  1:1126026/1‑71 (MQ=255)
  tATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCgg                                                <  1:442949/71‑1 (MQ=255)
                      tGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGctgctg                                        >  1:194725/1‑59 (MQ=255)
                      tGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGAc                            >  1:111141/1‑71 (MQ=255)
                      tGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGAc                            >  1:1674602/1‑71 (MQ=255)
                        aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGAc                            >  1:1058682/1‑69 (MQ=255)
                        aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGAc                            >  1:852775/1‑69 (MQ=255)
                        aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGAc                            >  1:1296899/1‑69 (MQ=255)
                        aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga                          <  1:820219/71‑1 (MQ=255)
                        aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga                          <  1:1955416/71‑1 (MQ=255)
                        aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCAGGAACTGTTGGACga                          <  1:480987/71‑1 (MQ=255)
                         tCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga                          <  1:133056/70‑1 (MQ=255)
                         tCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga                          <  1:1101086/70‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:1399222/71‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:903439/71‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:595040/71‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:1849722/71‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:396484/71‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:230063/71‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:1836496/71‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:157657/71‑1 (MQ=255)
                               aacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:1544058/71‑1 (MQ=255)
                                   gCTACGTAGCGGTACTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGAt                   <  1:1514287/67‑1 (MQ=255)
                                               tACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga                          <  1:895244/48‑1 (MQ=255)
                                                 cTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGc  >  1:1053581/1‑70 (MQ=255)
                                                 cTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGc  >  1:67827/1‑70 (MQ=255)
                                                 cTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGc  >  1:932247/1‑70 (MQ=255)
                                                      |                                                                
CGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGC  >  W3110S.gb/2916332‑2916450

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