Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,916,386 G→A 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,3860GA87.6% 244.7 / 22.8 89E892K (GAA→AAA) barAhybrid sensory histidine kinase, in two‑component regulatory system with UvrY
Reads supporting (aligned to +/- strand):  ref base G (11/0);  new base A (43/35);  total (54/35)
Fisher's exact test for biased strand distribution p-value = 2.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGA  >  W3110S.gb/2916319‑2916432
                                                                   |                                              
taGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             <  1:2652200/71‑1 (MQ=255)
taGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             <  1:3796176/71‑1 (MQ=255)
taGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             <  1:2391056/71‑1 (MQ=255)
taGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             <  1:90143/71‑1 (MQ=255)
  gCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             >  1:2764182/1‑69 (MQ=255)
  gCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             >  1:2916467/1‑69 (MQ=255)
  gCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             >  1:3294502/1‑69 (MQ=255)
  gCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             >  1:4097073/1‑69 (MQ=255)
  gCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             >  1:4079039/1‑69 (MQ=255)
  gCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             >  1:2484302/1‑69 (MQ=255)
  gCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAaaaaag                                             >  1:3155807/1‑69 (MQ=255)
          ccACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTgg                                    >  1:2761361/1‑70 (MQ=255)
            aCGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAAc                                 <  1:4017359/71‑1 (MQ=255)
            aCGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAAc                                 <  1:2015982/71‑1 (MQ=255)
                       aTCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGAt                                        >  1:1813333/1‑53 (MQ=255)
                       aTCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGctgctg                      >  1:2877148/1‑71 (MQ=255)
                       aTCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGctgctg                      >  1:1672270/1‑71 (MQ=255)
                       aTCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGctgctg                      >  1:1618317/1‑71 (MQ=255)
                       aTCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGctgctg                      >  1:1603758/1‑71 (MQ=255)
                       aTCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGctgctg                      >  1:2619945/1‑71 (MQ=255)
                       aTCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGctgctg                      >  1:1247546/1‑71 (MQ=255)
                       aTCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGctgctg                      >  1:4080987/1‑71 (MQ=255)
                         ctctGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTgg                     >  1:1887188/1‑70 (MQ=255)
                         ctctGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTgg                     >  1:4100914/1‑70 (MQ=255)
                         ctctGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTgg                     >  1:3240644/1‑70 (MQ=255)
                         ctctGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTgg                     >  1:3407814/1‑70 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTgg                     <  1:738112/67‑1 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTgg                     <  1:37255/67‑1 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTgg                     <  1:796777/67‑1 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGa                    <  1:3925294/68‑1 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGa                    <  1:489059/68‑1 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGa                    <  1:3275801/68‑1 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGa                    <  1:1911387/68‑1 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGa                    <  1:2539256/68‑1 (MQ=255)
                            tGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACt                 >  1:1259557/1‑71 (MQ=255)
                              ccAACTGATCGAACAAGAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACt                 <  1:93907/69‑1 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACt                 <  1:3324990/69‑1 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACt                 <  1:3329689/69‑1 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACt                 <  1:316910/69‑1 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACt                 <  1:2676879/69‑1 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACt                 <  1:2738976/69‑1 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:1191124/1‑71 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:1216202/1‑71 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:401498/1‑71 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:3539403/1‑71 (MQ=255)
                              ccAACTGATCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACAGGAGCTGCTGGAACTGt               >  1:3978429/1‑71 (MQ=255)
                                     aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGAc          >  1:3240650/1‑69 (MQ=255)
                                     aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGAc          >  1:3273665/1‑69 (MQ=255)
                                     aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga        <  1:3745524/71‑1 (MQ=255)
                                     aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga        <  1:4107091/71‑1 (MQ=255)
                                     aTCGAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga        <  1:1715915/71‑1 (MQ=255)
                                        gAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga        <  1:2700372/68‑1 (MQ=255)
                                        gAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga        <  1:1388787/68‑1 (MQ=255)
                                        gAACAACAGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga        <  1:3306962/68‑1 (MQ=255)
                                         aacaacAGCTACGTAGCTGTACTAAAAAAGAAGATTTgg                                    <  1:3886567/39‑1 (MQ=38)
                                         aacaacAGCTACGTAGCGGTACTAAAAAAGAAGATTTgg                                    <  1:3108568/39‑1 (MQ=255)
                                         aacaacAGCTACGTAGCGGTACTAAAAAAGAAGATTTgg                                    <  1:2158939/39‑1 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGTAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:839903/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAgc                           >  1:3100323/1‑47 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:4028522/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:4008324/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:3517300/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:660924/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:2736126/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:3055471/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:21136/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:3110537/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgaga      >  1:3195184/1‑68 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACgag       >  1:1745443/1‑67 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:1108519/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:278334/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:2812107/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:2833865/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:593776/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:800699/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:3678738/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:3612950/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:1969210/1‑70 (MQ=255)
                                          acaacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    >  1:1834413/1‑70 (MQ=255)
                                           caacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGa  <  1:767130/71‑1 (MQ=255)
                                           caacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGa  <  1:2471094/71‑1 (MQ=255)
                                           caacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGa  <  1:784575/71‑1 (MQ=255)
                                           caacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGa  <  1:3583327/71‑1 (MQ=255)
                                           caacaGCTACGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGa  <  1:3482085/71‑1 (MQ=255)
                                                   aCGTAGCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACga        <  1:336422/57‑1 (MQ=255)
                                                        gCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    <  1:3443497/56‑1 (MQ=255)
                                                        gCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    <  1:460768/56‑1 (MQ=255)
                                                        gCGGTACTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATg    <  1:846937/56‑1 (MQ=255)
                                                              cTAAAAAAGAAGATTTGGAACCGGAGCTGCTGGAACt                 <  1:3363915/37‑1 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:961140/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:485577/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:480753/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:3606911/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:3501110/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:3407996/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:261091/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:2350297/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:1884860/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:1252903/1‑37 (MQ=255)
                                                                aaaGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGt               >  1:1140289/1‑37 (MQ=255)
                                                                   |                                              
TAGCGGTGTGCCACGTATGAAGAATCTCTGCCAACTGATCGAACAACAGCTACGTAGCGGTACTAAAGAAGAAGATTTGGAACCGGAGCTGCTGGAACTGTTGGACGAGATGGA  >  W3110S.gb/2916319‑2916432

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: