Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,484,897 (C)6→7 coding (524/3594 nt) evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,484,8911.C95.6% 152.1 / ‑0.6 45coding (518/3594 nt)evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (21/22);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAG  >  NC_000913/2484754‑2485023
                                                                                                                                           |                                                                                                                                     
ttCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                     >  2:281465/1‑139 (MQ=255)
        ctctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCCCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑CCCCCCAGa                                                                                                                             <  2:790139/139‑1 (MQ=255)
                aatGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTa                                                                                                                      <  1:849952/139‑1 (MQ=255)
                    aCATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtc                                                                                                                  >  1:564376/1‑139 (MQ=255)
                     cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtca                                                                                                                 >  2:361674/1‑139 (MQ=255)
                     cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtca                                                                                                                 >  1:1252456/1‑139 (MQ=255)
                       ttGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                           <  2:872395/127‑1 (MQ=255)
                       ttGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                           >  1:872395/1‑127 (MQ=255)
                        tgctgcAACCAAAACATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtcatca                                                                                                              >  2:904285/1‑139 (MQ=255)
                           tgcAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTGCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCatcaatc                                                                                                           <  2:1153975/139‑1 (MQ=255)
                               accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                     <  2:736196/108‑1 (MQ=255)
                               accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                     >  1:736196/1‑108 (MQ=255)
                                    accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                                  <  1:361674/139‑1 (MQ=255)
                                     ccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaaaa                                                                                                 <  1:904285/139‑1 (MQ=255)
                                                  ccTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                     <  1:100454/89‑1 (MQ=255)
                                                  ccTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑c                                                                                                                                     >  2:100454/1‑89 (MQ=255)
                                                     ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                           <  1:1269777/97‑1 (MQ=255)
                                                     ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACg                                                                                                                           >  2:1269777/1‑97 (MQ=255)
                                                     ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTaa                                                                                                                     <  1:435030/103‑1 (MQ=255)
                                                     ttCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTaa                                                                                                                     >  2:435030/1‑103 (MQ=255)
                                                      tCCGGCGCCGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                        >  1:566357/1‑84 (MQ=255)
                                                      tCCGGCGCCGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                        <  2:566357/84‑1 (MQ=255)
                                                         ggCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                             <  1:281465/139‑1 (MQ=255)
                                                          gcgcTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                             >  1:1332381/1‑139 (MQ=255)
                                                          gcgcTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACa                                                                            >  2:680556/1‑139 (MQ=255)
                                                                       cccTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCa                                                                                                    <  1:534971/102‑1 (MQ=255)
                                                                       cccTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCa                                                                                                    >  2:534971/1‑102 (MQ=255)
                                                                       accTTCACGGCTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTACGGTAGCAAAATAA‑‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAgg                                                               <  2:557574/138‑1 (MQ=255)
                                                                          ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTa                                                                                                                      <  1:166573/81‑1 (MQ=255)
                                                                          ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTa                                                                                                                      >  2:166573/1‑81 (MQ=255)
                                                                            cACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTa                                                          >  2:1286887/1‑139 (MQ=255)
                                                                               gACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                             <  2:573558/117‑1 (MQ=255)
                                                                               gACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                             >  1:573558/1‑117 (MQ=255)
                                                                               gACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                       <  1:1286887/139‑1 (MQ=255)
                                                                                aCTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                       <  2:1332381/139‑1 (MQ=255)
                                                                                      tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                               >  1:239602/1‑92 (MQ=255)
                                                                                      tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                               <  2:239602/92‑1 (MQ=255)
                                                                                         gACCGCTTACCTCACCAAAACCAATAAATATTGCTCGGGTAGCAAATTA‑‑cccccc                                                                                                                                <  2:1141665/55‑1 (MQ=255)
                                                                                         gACCGCTTACCTCACCAAAACCAATAAATATTGCTCGGGTAGCAAATTA‑‑cccccc                                                                                                                                >  1:1141665/1‑55 (MQ=255)
                                                                                           ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAg                                                                                                >  1:1345298/1‑86 (MQ=255)
                                                                                           ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAg                                                                                                <  2:1345298/86‑1 (MQ=255)
                                                                                                 aCCTCACCAAAACCAGTAGATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGcac                                     <  1:905534/139‑1 (MQ=255)
                                                                                                      aCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                   <  2:968562/104‑1 (MQ=255)
                                                                                                      aCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                   >  1:968562/1‑104 (MQ=255)
                                                                                                           aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                               <  2:508419/135‑1 (MQ=255)
                                                                                                           aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAt                               >  1:508419/1‑135 (MQ=255)
                                                                                                             ccAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACttt                         >  1:1342390/1‑139 (MQ=255)
                                                                                                                      atTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGgtagt                <  2:1252456/139‑1 (MQ=255)
                                                                                                                             gggTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCCTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAcatca         >  1:1128778/1‑139 (MQ=255)
                                                                                                                               gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt                                                                                      <  2:1287607/60‑1 (MQ=255)
                                                                                                                               gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt                                                                                      >  1:1287607/1‑60 (MQ=255)
                                                                                                                                  gCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTgg                    >  2:292588/1‑123 (MQ=255)
                                                                                                                                  gCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTgg                    <  1:292588/123‑1 (MQ=255)
                                                                                                                                     aaaTTACCCCCCCAGACTAGGTAATTCATCAATCATTTCCAAAATCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACcag  <  2:989153/139‑1 (MQ=255)
                                                                                                                                     aaaTTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACcag  <  2:1342390/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                     
TTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTA‑‑CCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAG  >  NC_000913/2484754‑2485023

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: