Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,165,217 A→T 15.3% intergenic (+25/+21) qseC → / ← ygiZ quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,165,2170AT15.3% 53.5 / 4.5 25intergenic (+25/+21)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (12/9);  new base T (2/2);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

CCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATA  >  CP009273/3165080‑3165358
                                                                                                                                         |                                                                                                                                             
ccGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAct                                                                                                                                        >  2:156164/1‑145 (MQ=255)
           cTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAttt                                                                                                                             <  2:303034/145‑1 (MQ=255)
               ggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagagaa                                                                                                                         <  1:283561/145‑8 (MQ=38)
                             cAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATCGATTTTACTCATAGCAACAtata                                                                                                           >  1:65248/1‑145 (MQ=255)
                             cAGCGAATCGCCAAATTGCATGGCATGAATGTTAAATTTGGGAATGCGGAACAAGGTGGAAATGCGGCGAAGGTCAGCTGCTTACATCTCTGTCTATGAGAGCCGTGGAAACGACTCTCATAGATTTTACTAATAGCAAAataga                                                                                                           >  2:310180/1‑143 (MQ=255)
                                     cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTACTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcc                                                                                                   <  1:353621/145‑1 (MQ=255)
                                      gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGGTTTGAGGCGAAGGTAAGCTGGTAACTTCGCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                                                                  <  1:100767/145‑1 (MQ=255)
                                               catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATACCAAAATATAAACCGTCCCCAAAAAAGc                                                                                         >  2:284446/1‑145 (MQ=255)
                                                      tgAATGTTGAATTTTGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                                                                  <  2:44376/145‑1 (MQ=255)
                                                                   tGGGAATGCGGAAAAAGGTCGATTTCACGGCAAGGTAAGCTGCTAACTTCTCTGTCTATGAGAGCCGTTATAACTACTATCATAGATTTTAATCCTAACTAAACATACACCGCCCATAAAAAAGCCACCAACCACAAACCAGACa                                                                     <  1:38648/145‑1 (MQ=255)
                                                                     ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATATATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                                                                   <  2:329244/145‑1 (MQ=255)
                                                                      gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAt                                                                  >  2:161153/1‑145 (MQ=255)
                                                                         tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTACTAGCAAAATATAAACAGACCCCAAATAATCCACCAACCTCAAAACAGACAAATAcc                                                               >  2:231854/1‑145 (MQ=255)
                                                                                           tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAc                                             <  1:310180/145‑1 (MQ=255)
                                                                                               gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACACACCAGACAAATACAGAGCGAGAATACATATACCCaa                                         >  2:356695/1‑145 (MQ=255)
                                                                                               gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaa                                         >  2:135751/1‑145 (MQ=255)
                                                                                                       aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAGGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                 <  2:149185/145‑1 (MQ=255)
                                                                                                                ttagtaaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCg                        >  2:422573/10‑145 (MQ=38)
                                                                                                                       tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACAAGACAAAAAACGAGAGAGAATACATATACGCTAAACCTGAACATTTCCCGCATtaat                 >  2:124986/1‑144 (MQ=255)
                                                                                                                           tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATATCCCGCATAAAGAGtt             >  2:411335/1‑145 (MQ=255)
                                                                                                                            aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGACAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGttt            >  2:56903/1‑145 (MQ=255)
                                                                                                                             tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt           <  1:135751/145‑1 (MQ=255)
                                                                                                                                 agCCGTTAAAACGACTCTCACAGATTCTACCAATAGCACAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta       <  1:43251/145‑1 (MQ=255)
                                                                                                                                    cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaaga    >  2:278669/1‑145 (MQ=255)
                                                                                                                                      ttAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagata  >  2:361620/1‑145 (MQ=255)
                                                                                                                                         |                                                                                                                                             
CCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATA  >  CP009273/3165080‑3165358

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