Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,291,910 G→A 13.8% H273H (CAC→CAT)  PP_3758 ← two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,291,9100GA13.8% 75.5 / 3.8 29H273H (CAC→CAT) PP_3758two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (14/11);  new base A (2/2);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGAAACATCGAGCATG  >  NC_002947/4291775‑4292036
                                                                                                                                       |                                                                                                                              
cTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGAtg                                                                                                                                <  1:150528/136‑1 (MQ=255)
  gCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGAtgtg                                                                                                                              >  2:125178/1‑136 (MQ=255)
      gAACACCCGCACCTACACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGtt                                                                                                                          >  1:392740/1‑136 (MQ=255)
        acacCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGAtgt                                                                                                                               >  1:421103/1‑129 (MQ=255)
        acacCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGAtgt                                                                                                                               <  2:421103/129‑1 (MQ=255)
                     cACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGcc                                                                                                           >  2:501159/1‑136 (MQ=255)
                     cACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGcc                                                                                                           >  2:181407/1‑136 (MQ=255)
                           aTTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGcc                                                                                                     <  2:392740/136‑1 (MQ=255)
                           aTTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGcc                                                                                                     <  2:289588/136‑1 (MQ=255)
                                gcTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGaa                                                                                                >  2:310558/1‑136 (MQ=255)
                                       gTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTcc                                                                                         <  2:44981/136‑1 (MQ=255)
                                       gTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTcc                                                                                         >  2:496994/1‑136 (MQ=255)
                                                 gcattgGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCAc                                                                               >  2:420954/1‑136 (MQ=255)
                                                   attgGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGc                                                                             >  1:465765/1‑136 (MQ=255)
                                                          aacaGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCt                                                                      <  1:308159/136‑1 (MQ=255)
                                                               gTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCAATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGccc                                                                 >  1:406113/1‑136 (MQ=255)
                                                               gTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATATGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCAc                                                                               <  2:428271/122‑1 (MQ=255)
                                                               gTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATATGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCAc                                                                               >  1:428271/1‑122 (MQ=255)
                                                                                     cGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTg                                                  >  2:108009/1‑129 (MQ=255)
                                                                                     cGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTg                                                  <  1:108009/129‑1 (MQ=255)
                                                                                                ccACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTa                                >  2:105908/1‑136 (MQ=255)
                                                                                                        ccTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCg                        <  2:465765/136‑1 (MQ=255)
                                                                                                              gTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATgcg                  >  1:200159/1‑136 (MQ=255)
                                                                                                              gTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATgcg                  >  2:76745/1‑136 (MQ=255)
                                                                                                                 cTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCCCCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGa               >  1:206873/1‑136 (MQ=255)
                                                                                                                  ttCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCAATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGaa              <  2:406113/136‑1 (MQ=255)
                                                                                                                   gcGCACCGCTTCACTTCGATGTGAGTTTCCAGGGCCATGGCCTGCGCCGAGAAGTTATCCACCAGCCCACGCACCACCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGaaa             <  2:200159/135‑1 (MQ=255)
                                                                                                                     gCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGAAACa           >  1:99836/1‑136 (MQ=255)
                                                                                                                       aaCGCTGCAATTCGATATGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGAAACATc         >  1:256158/1‑136 (MQ=255)
                                                                                                                        aCGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGAAACATCg        <  1:501159/136‑1 (MQ=255)
                                                                                                                        aCGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGAAACATCg        <  1:181407/136‑1 (MQ=255)
                                                                                                                        aCGCTGCAATTCGATATGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGAAACATCg        <  2:256158/136‑1 (MQ=255)
                                                                                                                              cAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGAAACATCGAGCATg  <  1:440043/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                              
CTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATTGCATTGGACAACAGGTTGGCCACTACTTGTTCAATCCGGAATTCGTCCCACTCCCCCTTCAGTGCTTCGCAACGCTGCAATTCGATGTGTGTTTCCAGGGCCATGGCCTGCGCCGTGAAGTTTTCCACCAGCCCACGCACCAGCTGGCCCAGGTCAAACAGTTTGGGCCGCAGCGAAAGTTTACCGGTCCGGATGCGCGAAACATCGAGCATG  >  NC_002947/4291775‑4292036

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: ‑