Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,852 C→T 22.0% intergenic (+19/+27) 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
*NC_0009133,171,8520CT22.0% 100.5 / 25.9 50intergenic (+19/+27)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (17/22);  new base T (5/6);  total (22/28)
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.35e-01

ACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAA  >  NC_000913/3171704‑3172000
                                                                                                                                                    |                                                                                                                                                    
aCCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGc                                                                                                                                                      <  3:226519/149‑1 (MQ=255)
 cccGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGcc                                                                                                                                                     <  4:70351/149‑1 (MQ=255)
 cccGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGcc                                                                                                                                                     <  3:212677/149‑1 (MQ=255)
        aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                              >  1:514738/1‑149 (MQ=255)
        aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                              >  3:676321/1‑149 (MQ=255)
        aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTAGTAACTTCTCTGTCTATGAGAGTCGttttaa                                                                                                                                              >  1:162108/1‑145 (MQ=255)
         aaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                             <  3:377416/149‑1 (MQ=255)
         aaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                             <  4:663630/149‑1 (MQ=255)
          aCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACg                                                                                                                                            >  3:589229/1‑149 (MQ=255)
          aCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACg                                                                                                                                            >  2:636550/1‑149 (MQ=255)
               tACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTAGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGctct                                                                                                                                       <  2:162108/149‑1 (MQ=255)
                 ccGGCAGCGGGCTTGGGGTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                                                                     <  4:479803/149‑1 (MQ=255)
                 ccGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCa                                                                                                                                     <  2:358718/149‑1 (MQ=255)
                  cGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                                                                    <  4:589229/149‑1 (MQ=255)
                   ggCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                                                                   <  1:187406/149‑1 (MQ=255)
                   ggCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                                                                   <  1:636550/149‑1 (MQ=255)
                   ggCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                                                                   <  3:108182/149‑1 (MQ=255)
                      aGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGATTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                                                                                >  4:48236/1‑149 (MQ=255)
                      aGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                                                                                <  1:183922/149‑1 (MQ=255)
                         gggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                              <  1:60935/132‑1 (MQ=255)
                         gggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                              >  2:60935/1‑132 (MQ=255)
                          ggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                            <  4:676321/149‑1 (MQ=255)
                           gCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                                                           >  2:41394/1‑149 (MQ=255)
                                  cTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                                                                    >  1:476623/1‑149 (MQ=255)
                                            tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                          >  4:486550/1‑149 (MQ=255)
                                               aGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                                                                       >  3:160605/1‑149 (MQ=255)
                                                  gAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                        <  4:612823/113‑1 (MQ=255)
                                                  gAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                        >  3:612823/1‑113 (MQ=255)
                                                   aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                                                                   <  2:476623/149‑1 (MQ=255)
                                                          aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                            <  4:636385/149‑1 (MQ=255)
                                                          aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                            >  3:636385/1‑149 (MQ=255)
                                                              tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                        <  4:160605/149‑1 (MQ=255)
                                                                 atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                                                                     >  1:616951/1‑149 (MQ=255)
                                                                            gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                                                          <  3:48236/149‑1 (MQ=255)
                                                                                          tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                            <  2:616951/149‑1 (MQ=255)
                                                                                          tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                            <  3:486550/149‑1 (MQ=255)
                                                                                               aaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                       <  3:587163/149‑1 (MQ=255)
                                                                                                      tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                                                                          >  4:444385/1‑59 (MQ=255)
                                                                                                      tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                                                                          <  3:444385/59‑1 (MQ=255)
                                                                                                                    aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                                                          <  2:605652/125‑1 (MQ=255)
                                                                                                                    aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTCTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                                                          >  1:605652/1‑125 (MQ=255)
                                                                                                                           cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTcc                           <  2:514738/149‑1 (MQ=255)
                                                                                                                            tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTccc                          >  2:673367/1‑149 (MQ=255)
                                                                                                                             ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCg                         <  1:41394/149‑1 (MQ=255)
                                                                                                                                          tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACAtata                                               >  4:137593/1‑114 (MQ=255)
                                                                                                                                          tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACAtata                                               <  3:137593/114‑1 (MQ=255)
                                                                                                                                             aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt            <  2:217285/146‑1 (MQ=255)
                                                                                                                                             aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt            >  1:217285/1‑146 (MQ=255)
                                                                                                                                              tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta        >  4:3824/1‑149 (MQ=255)
                                                                                                                                                agagTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAAc                                 <  1:428686/122‑1 (MQ=255)
                                                                                                                                                agagTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAAc                                 >  2:428686/1‑122 (MQ=255)
                                                                                                                                                 gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa       >  1:231469/1‑147 (MQ=255)
                                                                                                                                                 gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa       <  2:231469/147‑1 (MQ=255)
                                                                                                                                                    ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataa  >  2:194346/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAA  >  NC_000913/3171704‑3172000

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑