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.4 / ‑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

GGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGT  >  NC_000913/418818‑419097
                                                                                                                                                |                                                                                                                                         
ggTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTGCGAACCCGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACg                                                                                                                                           <  1:221425/143‑1 (MQ=255)
    tCTGGCAACGAAGATCGGCTACGCAGTGCGATTTCGCACCTGGTCTATACCGCCGTGACTGCAACGCTGGAAGGCACGCATATCACCGTACCCTGGCAGGGAGTGCTGATCGGTGCCGAATTTAGCGTTGAAGATAACgaccg                                                                                                                                       <  2:617199/143‑5 (MQ=255)
      tGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGCAGGCACGCATATCACCGTACGCTGGCAGCGAGCGCCGCACGGTGCCGAATTTAGCGGTGAAGATGACgtccggg                                                                                                                                     >  2:212094/1‑137 (MQ=255)
      tGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCAAGCCCAATTAGAGATTGTTGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCggg                                                                                                                                    >  1:927819/1‑143 (MQ=255)
       ggCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGGGACGAAGTTAGCGTTGAAGATAACGACCggga                                                                                                                                   >  2:305381/1‑142 (MQ=255)
             gAAGATCAGCTCCGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATGTAGCGTTGAAGATAACGACCGGGCACTGCa                                                                                                                             >  2:250711/1‑143 (MQ=255)
             gAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCGGCACGGTGCCGAATTTAGCGTTGAAGATGACGACCGGGCATTGCa                                                                                                                             >  1:80808/1‑143 (MQ=255)
                  tCATCTACTGAGAGCGAGTACGAACCTGGTCTATACCGCCGTGAATCATACGCCGGAAGGCACGCAGATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGa                                                                                                                        <  2:669226/143‑1 (MQ=255)
                          gCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGAGCCGAACGGTGCCGAACTTAGCGTGGAAGCAAACGACCGGGCAATGCACCGGAGCAGATTc                                                                                                                >  2:878584/1‑143 (MQ=255)
                                   tttCGAAGCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGa                                                                                                       <  1:200531/143‑1 (MQ=255)
                                   tttCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATCTCACCGTACGCTGGCAGCGCGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCCTCGCACCCGCGCCTCTCCCGCGCCtca                                                                                                       >  2:919675/1‑141 (MQ=255)
                                   tttCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGGGCATATTCCGCGCCTGa                                                                                                       >  1:391813/1‑143 (MQ=255)
                                    ttCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGGGCCGAATTTAGCGTTGACGATAACGACCGGGCGTCGCACCGGGGCATATTCCgcgcatcag                                                                                                      >  1:913262/1‑138 (MQ=255)
                                       gAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGa                                                                                                   <  2:919437/143‑1 (MQ=255)
                                       gAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGa                                                                                                   >  1:243838/1‑143 (MQ=255)
                                            tGGTCTATAACGCAGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGtt                                                                                              <  2:773151/143‑1 (MQ=255)
                                                cTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTa                                                                                          <  1:618197/143‑1 (MQ=255)
                                                 tataACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGCCGACGAGTGCCACACGGCGCAGTTGTTAGCGTTGAAGATCACGCCCGGGCATTGCACCGGCGCATCTTCCGCGCCTTGCCGAGCGtttctct                                                                                         >  2:329196/1‑139 (MQ=255)
                                                  ataACGCCGTGAATCATACGCCGTATTGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATGTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAACATATTCCGCGCCTGACCGAGCGCTTTTATc                                                                                        >  2:116240/1‑143 (MQ=255)
                                                   taACGCCGTGGATCATAGGGCGGAAGGAACGCATACCACCGTACGGTGGCAGCGAGTGCAGCTCGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGACCATATTCCGCGCCTGACCGAGCGTTTTTATcg                                                                                       <  2:505448/143‑1 (MQ=255)
                                                   gaaCGCCGTGAGTCATACGCCGGAAGGCACGCAGATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATcg                                                                                       <  2:585731/142‑1 (MQ=255)
                                                    aaCGCCGTGAATCATACGCCGGAAGGCACGAGCACAGGCGGACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATcgc                                                                                      <  1:723494/143‑1 (MQ=255)
                                                       gCCTTGAATCATACGCCGCAATGCACGCATATCACCGCACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGtt                                                                                   <  1:25390/143‑1 (MQ=255)
                                                       gCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGtt                                                                                   <  1:860344/143‑1 (MQ=255)
                                                       gCCGTGAATCAGACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGtt                                                                                   <  1:247136/143‑1 (MQ=255)
                                                        ccGTGAATCCTACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGGTAACGACCGGGCCTTGCACCGGAGCATATTCCGCGCCTGACCGAGCGACTTTATCGCGTTg                                                                                  >  2:216695/1‑143 (MQ=255)
                                                            gAACCATACGCCGGAAGGCACGCAGATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGGAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATaa                                                                              <  1:343126/143‑1 (MQ=255)
                                                               tCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGGATTGCAGCGGAGCGTATTCCGCGCCTGGCCGAGCGTGTTTATCGCGTTGATAAAgc                                                                           >  1:5601/1‑143 (MQ=255)
                                                                    cGCCGGAAGGCACGCATATCACCGTACGCTGGCAGAGAGTGTGGGGTGGTGGGGGATCAAGAGATCACACTTAAC‑‑AACTATCTTTGCACCGGAGCATGTTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGtt                                                                      >  1:931682/1‑143 (MQ=255)
                                                                                   cATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGTGTTGATAAAGCGCGTTCCCGGCAAACcggc                                                        >  1:687923/1‑143 (MQ=255)
                                                                                        cACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAg                                                   >  1:101709/1‑143 (MQ=255)
                                                                                        cACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTCACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCTCGACACACCGGCGGTAg                                                   >  2:405170/1‑143 (MQ=255)
                                                                                                 cgggcggcCAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAgg                                          <  2:80808/141‑1 (MQ=255)
                                                                                                           gTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATc                                <  1:908369/143‑1 (MQ=255)
                                                                                                              ccGCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTGGGGTTAGCGATCGTg                             >  1:196063/1‑143 (MQ=255)
                                                                                                               ggCACGGTGCCGAATTCAGCGTTTAAGTTAACGGCCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTGTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGa                            <  2:85275/142‑1 (MQ=255)
                                                                                                                gCACGGTGCCGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGaa                           >  1:760329/1‑143 (MQ=255)
                                                                                                                        ccGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCtg                   >  1:661283/1‑143 (MQ=255)
                                                                                                                        ccGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGAGTAGAGTTAGCGAGCGTGAAACATGCtg                   >  2:44654/1‑143 (MQ=255)
                                                                                                                        ccGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGGTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGTCGGTAGCGGATGAGGGTTAGCGATCGTGAAACATGCtg                   >  2:220771/1‑143 (MQ=255)
                                                                                                                         cGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCtgt                  <  1:790997/143‑1 (MQ=255)
                                                                                                                         cGAATTTAGCGTTGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCtgt                  >  1:633740/1‑143 (MQ=255)
                                                                                                                                  gggtGAAGATAACGACCGGGCATTGCACCGGAGGATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAAtcatca         <  1:632810/140‑1 (MQ=255)
                                                                                                                                  cggtGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAAtcatca         <  1:725241/140‑1 (MQ=255)
                                                                                                                                     tGAAGATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGa      >  2:502802/1‑143 (MQ=255)
                                                                                                                                          gATAACGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGt  >  1:529962/1‑143 (MQ=255)
                                                                                                                                                |                                                                                                                                         
GGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGT  >  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: