Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 418,961 Δ1 bp coding (1073/1296 nt) phoR → sensory histidine kinase in two‑component regulatory system with PhoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913418,9600G.97.7% 179.8 / ‑3.2 43coding (1072/1296 nt)phoRsensory histidine kinase in two‑component regulatory system with PhoB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (25/17);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGT  >  NC_000913/418818‑419097
                                                                                                                                                |                                                                                                                                         
ggTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTGCGAACCCGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAACg                                                                                                                                           <  1:221425/143‑1 (MQ=255)
    tCTGGCAACGAAGATCGGCTACGCAGTGCGATTTCGCACCTGGTCTATACCGCCGTGACTGCAACGCTGGAAGGCACG‑CATATCACCGTACCCTGGCAGGGAGTGCTGATCGGTGCCGAATTTAGCGTTGA‑AGATAACgaccg                                                                                                                                       <  2:617199/143‑5 (MQ=255)
      tGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGCAGGCACG‑CATATCACCGTACGCTGGCAGCGAGCGCCGCACGGTGCCGAATTTAGCGGTGA‑AGATGACgtccggg                                                                                                                                     >  2:212094/1‑137 (MQ=255)
      tGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCAAGCCCAATTAGAG‑ATTGTTGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCggg                                                                                                                                    >  1:927819/1‑143 (MQ=255)
       ggCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGGGACGAAGTTAGCGTTGA‑AGATAAC‑GACCggga                                                                                                                                   >  2:305381/1‑142 (MQ=255)
             gAAGATCAGCTCCGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATGTAGCGTTGA‑AGATAAC‑GACCGGGCACTGCa                                                                                                                             >  2:250711/1‑143 (MQ=255)
             gAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCGGCACGGTGCCGAATTTAGCGTTGA‑AGATGAC‑GACCGGGCATTGCa                                                                                                                             >  1:80808/1‑143 (MQ=255)
                  tCATCTACTGAGAGCGAGTACGAACCTGGTCTATACCGCCGTGAATCATACGCCGGAAGGCACG‑CAGATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGa                                                                                                                        <  2:669226/143‑1 (MQ=255)
                          gCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGAGCCGAACGGTGCCGAACTTAGCGTGGA‑AGCAAAC‑GACCGGGCAATGCACCGGAGCAGATTc                                                                                                                >  2:878584/1‑143 (MQ=255)
                                   tttCGAAGCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGa                                                                                                       <  1:200531/143‑1 (MQ=255)
                                   tttCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATCTCACCGTACGCTGGCAGCGCGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCCTCGCACCCGCGCCTCTCCCGCGCCtca                                                                                                       >  2:919675/1‑141 (MQ=255)
                                   tttCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGGGCATATTCCGCGCCTGa                                                                                                       >  1:391813/1‑143 (MQ=255)
                                    ttCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGGGCCGAATTTAGCGTTGA‑CGATAAC‑GACCGGGCGTCGCACCGGGGCATATTCCgcgcatcag                                                                                                      >  1:913262/1‑138 (MQ=255)
                                       gAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGa                                                                                                   <  2:919437/143‑1 (MQ=255)
                                       gAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGa                                                                                                   >  1:243838/1‑143 (MQ=255)
                                            tGGTCTATAACGCAGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGtt                                                                                              <  2:773151/143‑1 (MQ=255)
                                                cTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTa                                                                                          <  1:618197/143‑1 (MQ=255)
                                                 tataACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGCCGACGAGTGCCACACGGCGCAGTTGTTAGCGTTGA‑AGATCAC‑GCCCGGGCATTGCACCGGCGCATCTTCCGCGCCTTGCCGAGCGtttctct                                                                                         >  2:329196/1‑139 (MQ=255)
                                                  ataACGCCGTGAATCATACGCCGTATTGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATGTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAACATATTCCGCGCCTGACCGAGCGCTTTTATc                                                                                        >  2:116240/1‑143 (MQ=255)
                                                   taACGCCGTGGATCATAGGGCGGAAGGAACG‑CATACCACCGTACGGTGGCAGCGAGTGCAGCTCGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGACCATATTCCGCGCCTGACCGAGCGTTTTTATcg                                                                                       <  2:505448/143‑1 (MQ=255)
                                                   gaaCGCCGTGAGTCATACGCCGGAAGGCACG‑CAGATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATcg                                                                                       <  2:585731/142‑1 (MQ=255)
                                                    aaCGCCGTGAATCATACGCCGGAAGGCACG‑AGCACAGGCGGACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATcgc                                                                                      <  1:723494/143‑1 (MQ=255)
                                                       gCCTTGAATCATACGCCGCAATGCACG‑CATATCACCGCACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGtt                                                                                   <  1:25390/143‑1 (MQ=255)
                                                       gCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGtt                                                                                   <  1:860344/143‑1 (MQ=255)
                                                       gCCGTGAATCAGACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGtt                                                                                   <  1:247136/143‑1 (MQ=255)
                                                        ccGTGAATCCTACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGGTAAC‑GACCGGGCCTTGCACCGGAGCATATTCCGCGCCTGACCGAGCGACTTTATCGCGTTg                                                                                  >  2:216695/1‑143 (MQ=255)
                                                            gAACCATACGCCGGAAGGCACG‑CAGATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGG‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATaa                                                                              <  1:343126/143‑1 (MQ=255)
                                                               tCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGGATTGCAGCGGAGCGTATTCCGCGCCTGGCCGAGCGTGTTTATCGCGTTGATAAAgc                                                                           >  1:5601/1‑143 (MQ=255)
                                                                    cGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGAGAGTGTGGGGTGGTGGGGGATCAAGAGATCACACTTAAC‑‑AACTATCTTTGCACCGGAGCATGTTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGtt                                                                      >  1:931682/1‑143 (MQ=255)
                                                                                   cATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGTGTTGATAAAGCGCGTTCCCGGCAAACcggc                                                        >  1:687923/1‑143 (MQ=255)
                                                                                        cACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAg                                                   >  1:101709/1‑143 (MQ=255)
                                                                                        cACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTCACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCTCGACACACCGGCGGTAg                                                   >  2:405170/1‑143 (MQ=255)
                                                                                                 cgggcggcCAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAgg                                          <  2:80808/141‑1 (MQ=255)
                                                                                                           gTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATc                                <  1:908369/143‑1 (MQ=255)
                                                                                                              ccGCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTGGGGTTAGCGATCGTg                             >  1:196063/1‑143 (MQ=255)
                                                                                                               ggCACGGTGCCGAATTCAGCGTTTA‑AGTTAACGG‑CCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTGTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGa                            <  2:85275/142‑1 (MQ=255)
                                                                                                                gCACGGTGCCGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGaa                           >  1:760329/1‑143 (MQ=255)
                                                                                                                        ccGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCtg                   >  1:661283/1‑143 (MQ=255)
                                                                                                                        ccGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGAGTAGAGTTAGCGAGCGTGAAACATGCtg                   >  2:44654/1‑143 (MQ=255)
                                                                                                                        ccGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGGTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGTCGGTAGCGGATGAGGGTTAGCGATCGTGAAACATGCtg                   >  2:220771/1‑143 (MQ=255)
                                                                                                                         cGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCtgt                  <  1:790997/143‑1 (MQ=255)
                                                                                                                         cGAATTTAGCGTTGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCtgt                  >  1:633740/1‑143 (MQ=255)
                                                                                                                                  gggtGA‑AGATAAC‑GACCGGGCATTGCACCGGAGGATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAAtcatca         <  1:632810/140‑1 (MQ=255)
                                                                                                                                  cggtGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAAtcatca         <  1:725241/140‑1 (MQ=255)
                                                                                                                                     tGA‑AGATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGa      >  2:502802/1‑143 (MQ=255)
                                                                                                                                          gATAAC‑GACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGt  >  1:529962/1‑143 (MQ=255)
                                                                                                                                                |                                                                                                                                         
GGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACG‑CATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGA‑AGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGT  >  NC_000913/418818‑419097

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