Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,683,180 C→T 10.8% Q426Q (CAG→CAA glrK ← sensor protein kinase regulating glmY sRNA in two‑component system with response regulator GlrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,683,1800CT10.8% 92.2 / 2.5 37Q426Q (CAG→CAAglrKsensor protein kinase regulating glmY sRNA in two‑component system with response regulator GlrR
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (2/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGC  >  CP009273/2683042‑2683327
                                                                                                                                          |                                                                                                                                                   
tttCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATCCGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGACCTTGACCGCCCCTTTTCGCtggt                                                                                                                                                 <  1:465790/143‑1 (MQ=255)
          ggTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTg                                                                                                                                       >  1:61024/1‑143 (MQ=255)
           gTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCCCCGCCCCTTTTCGCTGGTGGCTTCCCTGa                                                                                                                                      >  1:135381/1‑143 (MQ=255)
             aaTTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGaaa                                                                                                                                    <  1:501692/143‑1 (MQ=255)
               ttCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGtt                                                                                                                            >  2:508925/1‑149 (MQ=255)
                   aTGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGaa                                                                                                                        >  2:540186/1‑149 (MQ=255)
                   aTGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATAAAGTTCCCCTTGCATAAGGCGAATACAATCCCTGGCAATGCTTAATCCAAAACCGCTGCCCTTCACCGCCCCTTTTCGCTGCTGGCTTCCCTGAACAACggg                                                                                                                              >  1:32994/1‑143 (MQ=255)
                        gAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCa                                                                                                                   <  2:98973/149‑1 (MQ=255)
                                  tCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGgc                                                                                                               >  1:489347/1‑143 (MQ=255)
                                  tCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTc                                                                                                         <  2:68595/149‑1 (MQ=255)
                                             cgcgtcgACTAGTTACATTTCCCCTTGCATACGGCGAATACCAGCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCCCCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGaa                                                                                              <  2:272603/147‑1 (MQ=255)
                                                        gATACCTTTCCCCTTGCCTACGGCGAATCCAACCCCTGGCCCTGCTTCCGGCCATACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCt                                                                                   <  2:303044/149‑1 (MQ=255)
                                                        gATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCt                                                                                   <  2:365877/149‑1 (MQ=255)
                                                              gtcccccGTGCCTCCTGCGAATGCAAGACCTGGCAATGCTGAACCCCAGTCCGCGGCCCTTCACCGCCCCGTTTCGCTGGTGGCTGGCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCt                                                                                   <  1:371700/140‑1 (MQ=255)
                                                                 ccccTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTa                                                                                >  1:572560/1‑143 (MQ=255)
                                                                     ttGCATACGGCGAATACAATCCCTGGCAATGCTTAGTCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATtga                                                                            >  1:353686/1‑143 (MQ=255)
                                                                     ttGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGTTGGTGGCTTCCATGACAAAAAGGCCAGCAGAGCAGGGCGCGTTCCTCTTGCGGAATGGGCGGGCCTGTAttta                                                                            >  1:459950/1‑141 (MQ=255)
                                                                     ttGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGTGGGTGGCTTCCCTGCCAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATtga                                                                            >  1:459785/1‑143 (MQ=255)
                                                                        cATACGGCGAATACAAGCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCaa                                                                   <  2:255579/149‑1 (MQ=255)
                                                                             ggCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATaa                                                              >  2:267625/1‑149 (MQ=255)
                                                                              gCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATaaa                                                             <  2:543436/149‑1 (MQ=255)
                                                                               cggaTACCATCCCTCGCACTGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGTTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAc                                                            <  2:511740/146‑1 (MQ=255)
                                                                               cGAATCCCATCGCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGTTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAc                                                            <  2:511551/149‑1 (MQ=255)
                                                                                   tACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGTCATCAatacag                                                              <  1:263296/143‑4 (MQ=255)
                                                                                      aaTCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCa                                                     <  2:289540/149‑1 (MQ=255)
                                                                                        tCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGGGCCGGTATTCCTGACATCAATATAAACCcg                                                         >  1:29603/1‑143 (MQ=255)
                                                                                           cTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCgcgc                                                      >  1:110352/1‑143 (MQ=255)
                                                                                           cTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAACCgcgc                                                      >  1:143349/1‑143 (MQ=255)
                                                                                                 aTGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCCATATAAACCCGCGCACCATGTAAAcg                                          >  2:291397/1‑148 (MQ=255)
                                                                                                     ttAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGCCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTaaa                                            <  1:95318/143‑1 (MQ=255)
                                                                                                      tAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATCTAAACCCGCGCACCATGTACACtgctgc                                     >  2:536039/1‑149 (MQ=255)
                                                                                                               aCCGCTGCCCTTCCCCGCCCCTTTCCGCCGGTGGCTTCTCCTAAAAAAGGGTTCGACGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGCATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGaa                                  <  1:400725/143‑1 (MQ=255)
                                                                                                                  gCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGt                         >  2:340044/1‑149 (MQ=255)
                                                                                                                    cgCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGCCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCaa                             <  1:32883/142‑1 (MQ=255)
                                                                                                                           cACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGAtt                <  2:534923/149‑1 (MQ=255)
                                                                                                                                    gtttCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATACACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTCCCGGATTCAg             >  1:7518/2‑143 (MQ=255)
                                                                                                                                         gCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGc  >  2:271499/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGC  >  CP009273/2683042‑2683327

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