Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,959,068 T→C 11.1% I120V (ATT→GTT)  PP_4363 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,959,0680TC11.1% 94.6 / 4.8 36I120V (ATT→GTT) PP_4363two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base C (2/2);  total (17/19)
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.41e-01

TTGTCGAACACCGCCTTGGCGACCCGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGACGTCAGGAAAATGATCGGCACCAGCGCTTCGCCAGC  >  NC_002947/4958947‑4959200
                                                                                                                         |                                                                                                                                    
ttGTCGAACACCGCCTTGGCGACCCGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGa                                                                                                                        <  2:119872/136‑1 (MQ=255)
   tCGAACACCGCCTTGGCGACCCGCTGCTCGTTAAGCAGGTGATTGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCa                                                                                                                     <  2:258234/136‑1 (MQ=255)
        cacCGCCTTGGCGACCCGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACcgcg                                                                                                                <  1:474957/136‑1 (MQ=255)
         acCGCCTTGGCGACCCGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCgcgc                                                                                                               <  1:290781/136‑1 (MQ=255)
                   cGACCCGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCtt                                                                                                     >  2:327341/1‑136 (MQ=255)
                        cGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAACCTTGGCAGCCAGGATCa                                                                                                                     <  2:237262/115‑1 (MQ=255)
                        cGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAACCTTGGCAGCCAGGATCa                                                                                                                     >  1:237262/1‑115 (MQ=255)
                        cGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAACCTTGGCAGCCAGGATCa                                                                                                                     >  1:240406/1‑115 (MQ=255)
                        cGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAACCTTGGCAGCCAGGATCa                                                                                                                     <  2:240406/115‑1 (MQ=255)
                                     gCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGcc                                                                                   <  2:616237/136‑1 (MQ=255)
                                     gCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGcc                                                                                   <  2:237522/136‑1 (MQ=255)
                                     gCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGcc                                                                                   <  1:273599/136‑1 (MQ=255)
                                             tggtgatgCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCt                                                                           <  1:532951/136‑1 (MQ=255)
                                               gtgatgCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTATGGCTTGGCCATGAAATCATCGCCCCCCGCCTcc                                                                         <  2:353783/136‑1 (MQ=255)
                                                  atgCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAgg                                                                      <  2:315444/136‑1 (MQ=255)
                                                      cGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGAt                                                                                                                       >  1:164695/1‑83 (MQ=255)
                                                      cGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGAt                                                                                                                       <  2:164695/83‑1 (MQ=255)
                                                      cGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCg                                                                              <  2:142746/124‑1 (MQ=255)
                                                      cGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCg                                                                              >  1:142746/1‑124 (MQ=255)
                                                      cGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAAc                                                                  >  1:464997/1‑136 (MQ=255)
                                                        ggTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGc                                                                <  1:349371/136‑1 (MQ=255)
                                                           gATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCAcc                                                             >  1:584171/1‑136 (MQ=255)
                                                           gATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCAcc                                                             >  2:175607/1‑136 (MQ=255)
                                                           gATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCAcc                                                             >  2:482073/1‑136 (MQ=255)
                                                               tGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGcc                                                         <  1:409739/136‑1 (MQ=255)
                                                                              aGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCACCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCttcattca                                          >  2:321793/1‑136 (MQ=255)
                                                                                   cGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTTGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGAc                                     >  1:59340/1‑136 (MQ=255)
                                                                                   cGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGa                                                                                                                        >  2:311243/1‑53 (MQ=255)
                                                                                   cGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGa                                                                                                                        <  1:311243/53‑1 (MQ=255)
                                                                                   cGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGAc                                     >  1:659471/1‑136 (MQ=255)
                                                                                       gCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGACGTCa                                 >  1:442293/1‑136 (MQ=255)
                                                                                       gCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGACGTCa                                 >  1:557515/1‑136 (MQ=255)
                                                                                       gCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGACGTCa                                 >  1:135294/1‑136 (MQ=255)
                                                                                         cTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGACGTCAgg                               <  2:557515/136‑1 (MQ=255)
                                                                                              aGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGACGTCAGGAAAAt                          <  1:313179/136‑1 (MQ=255)
                                                                                                                      gAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGCCAACGCACCACCCCCTCTTCTTCATTCAGCGACGTCAGGAAAAAGATCGGCACCAGCGCTACgccagc  >  2:492433/1‑136 (MQ=255)
                                                                                                                         |                                                                                                                                    
TTGTCGAACACCGCCTTGGCGACCCGCTGCTCGTTAAGCAGGTGATGGTGATGCCGGGTGATCTGGTCACGCTGTTCCAGTACCGTTGCCTGCAAGCGACGCAGGCGGTCCATGGCGCGAATCTTGGCAGCCAGGATCACCGCGCTGTAGGGCTTGGCCATGAAATCATCGCCCCCCGCCTCCAGGCAACGCACCAGCCCCTCTTCTTCATTCAGCGACGTCAGGAAAATGATCGGCACCAGCGCTTCGCCAGC  >  NC_002947/4958947‑4959200

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