Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,170,020 G→T 7.3% E65* (GAA→TAA)  qseB → quorum sensing DNA‑binding response regulator in two‑component regulatory system with QseC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,170,0200GT 7.3% 179.1 / 5.3 55E65* (GAA→TAA) qseBquorum sensing DNA‑binding response regulator in two‑component regulatory system with QseC
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base T (2/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

ATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCTTTTGCGTTGATAGAAGTCGCCGCCAGGCTGGAAGC  >  NC_000913/3169873‑3170168
                                                                                                                                                   |                                                                                                                                                    
aTCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGa                                                                                                                                                     >  1:291579/1‑149 (MQ=255)
 tCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGaa                                                                                                                                                    >  3:101286/1‑149 (MQ=255)
  cAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAAt                                                                                                                                                   <  2:306926/149‑1 (MQ=255)
   aaaaCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATg                                                                                                                                                  >  4:82132/1‑149 (MQ=255)
    aaaCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATgg                                                                                                                                                 <  4:101286/149‑1 (MQ=255)
    aaaCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATgg                                                                                                                                                 >  2:187570/1‑149 (MQ=255)
     aaCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGc                                                                                                                                                <  4:241048/149‑1 (MQ=255)
     aaCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGc                                                                                                                                                <  1:187570/149‑1 (MQ=255)
     aaCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGc                                                                                                                                                <  2:208538/149‑1 (MQ=255)
         ggCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCgaga                                                                                                                                            >  2:80890/1‑149 (MQ=255)
            cTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCACGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGaaaa                                                                                                                                         <  4:276961/149‑1 (MQ=255)
             ttAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGaaaaa                                                                                                                                        <  2:468970/149‑1 (MQ=255)
                 tAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCTGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTc                                                                                                                                    >  4:53177/1‑149 (MQ=255)
                     aTGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCTAATGGCGAGAAAAAGGTCAGCg                                                                                                                                >  3:318085/1‑149 (MQ=255)
                     aTGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCTAATGGCGAGAAAAAGGTCAGCg                                                                                                                                >  1:376000/1‑149 (MQ=255)
                     aTGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCg                                                                                                                                >  3:46828/1‑149 (MQ=255)
                                cGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTAc                                                                                                                     <  3:53177/149‑1 (MQ=255)
                                       tGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCt                                                                                                              <  1:80890/149‑1 (MQ=255)
                                          tttACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGAc                                                                                                           >  4:390759/1‑149 (MQ=255)
                                                          aGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTg                                                                                           >  2:290516/1‑149 (MQ=255)
                                                           gggAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTgg                                                                                          >  1:450887/1‑149 (MQ=255)
                                                           gggAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTgg                                                                                          >  4:123259/1‑149 (MQ=255)
                                                           gggAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTgg                                                                                          >  2:174487/1‑149 (MQ=255)
                                                           gggAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTgg                                                                                          >  2:283767/1‑149 (MQ=255)
                                                                 agagGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCTAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAAc                                                                                    <  2:376000/149‑1 (MQ=255)
                                                                 agagGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCTAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGACGCGCTGGCGGAAc                                                                                    <  4:318085/149‑1 (MQ=255)
                                                                     gcgcTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACgtgt                                                                                >  1:327984/1‑149 (MQ=255)
                                                                          ttATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAg                                                                           <  2:291579/149‑1 (MQ=255)
                                                                           tataGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAgg                                                                          <  1:283767/149‑1 (MQ=255)
                                                                            ataGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAggg                                                                         >  4:212178/1‑149 (MQ=255)
                                                                                  cACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGaa                                                                            <  3:403159/140‑1 (MQ=255)
                                                                                  cACCTTATGATGCGGTGATACTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGAACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGaa                                                                            >  4:403159/1‑140 (MQ=255)
                                                                                   aCCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTc                                                                  >  4:197080/1‑149 (MQ=255)
                                                                                      ttATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTgg                                                               >  3:209768/1‑149 (MQ=255)
                                                                                      ttATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTgg                                                               <  4:209768/149‑1 (MQ=255)
                                                                                      ttATGATGCGGCGATCCTGGATTTAACCTTGCCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAATAGGTGAGCGTGAGCCGGTACTGATTCTGACCGTGCGCGATGCGATGGCGGAACGTGTATAAGGGCTGCGTCTgg                                                               <  4:295011/149‑1 (MQ=255)
                                                                                                gTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGa                                                     <  2:450887/149‑1 (MQ=255)
                                                                                                  gATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGAtt                                                   <  1:290516/149‑1 (MQ=255)
                                                                                                   aTCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTa                                                  <  1:174487/149‑1 (MQ=255)
                                                                                                     ccTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATc                                                <  3:82132/149‑1 (MQ=255)
                                                                                                      cTGGATTTAACCTTACCAGGCAGGGATGGTCGCGATATTTTGCAAGAATGGCGAGGAAAAGGTCAGCGTGAGCCGGAACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCt                                               >  2:96468/1‑149 (MQ=255)
                                                                                                       tGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCtg                                              >  4:202958/1‑149 (MQ=255)
                                                                                                        ggATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCtgt                                             <  3:390759/149‑1 (MQ=255)
                                                                                                                 cTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTg                                                                <  1:142630/121‑1 (MQ=255)
                                                                                                                 cTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTg                                                                >  2:142630/1‑121 (MQ=255)
                                                                                                                 cTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCtt                                    >  3:126046/1‑149 (MQ=255)
                                                                                                                   tACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCtttt                                  <  2:327984/149‑1 (MQ=255)
                                                                                                                   tACCAGGCAAGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCtttt                                  >  3:375548/1‑149 (MQ=255)
                                                                                                                       aGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCt                                               <  2:452481/132‑1 (MQ=255)
                                                                                                                       aGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCt                                               >  1:452481/1‑132 (MQ=255)
                                                                                                                                      cgATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCTTTTGCGTTGATAGAAGTcgccg               >  3:465848/1‑149 (MQ=255)
                                                                                                                                       gATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCTTTTGCGTTGATAGAAGTcgccgc              >  2:448674/1‑149 (MQ=255)
                                                                                                                                             ttGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCt                                                          >  3:208323/1‑99 (MQ=255)
                                                                                                                                             ttGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCt                                                          <  4:208323/99‑1 (MQ=255)
                                                                                                                                               gcgcGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCTTTTGCGTTGATAGAAGTCGCCGCCAGGCTgg      <  4:46828/149‑1 (MQ=255)
                                                                                                                                                 gcgAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCTTTTGCGTTGATAGAAGTCGCCGCCAGGCTGGaa    >  3:360945/1‑149 (MQ=255)
                                                                                                                                                   gAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCTTTTGCGTTGATAGAAGTCGCCGCCAGGCTGGAAGc  >  4:126381/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATTTTGCGCGAATGGCGAGAAAAAGGTCAGCGTGAGCCGGTACTGATCCTGACCGCGCGCGATGCGCTGGCGGAACGTGTAGAAGGGCTGCGTCTGGGAGCTGACGATTATCTGTGTAAACCTTTTGCGTTGATAGAAGTCGCCGCCAGGCTGGAAGC  >  NC_000913/3169873‑3170168

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