Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,051,593 G→A 9.5% R98H (CGC→CAC)  PP_2664 → 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_0029473,051,5930GA 9.5% 115.1 / 2.9 42R98H (CGC→CAC) PP_2664two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base A (2/2);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GCGTTCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGATGCTGGGCCAGAGCCTGCGCGAGGCCCACAG  >  NC_002947/3051448‑3051734
                                                                                                                                                 |                                                                                                                                             
gcgTTCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCt                                                                                                                                            <  5:2766/149‑1 (MQ=255)
   ttCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGAt                                                                                                                                         <  2:11619/149‑1 (MQ=255)
      cAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTAc                                                                                                                                      >  1:126895/1‑149 (MQ=255)
      cAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTAc                                                                                                                                      >  1:199398/1‑149 (MQ=255)
      cAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTAc                                                                                                                                      <  7:270830/149‑1 (MQ=255)
              cgtggtgCTGGCCGCGCAGGTACGCGAACGCACCGAGGCTTTGAACCCGGCCATGGCCGAGCTCAAGGCCTGCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACg                                                                                                                              <  2:126895/148‑1 (MQ=255)
              aGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAAc                                                                                                                               <  1:197776/148‑1 (MQ=255)
                 ggtgCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTcc                                                                                                                           <  5:238393/149‑1 (MQ=255)
                  gtgCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTcc                                                                                                                           <  2:37823/148‑1 (MQ=255)
                                ggTACGCGAACGCACCGATGCGTTGAACCAGGCCATGACCGAGCTCACGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGACGGCCCCCCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGCTCGCt                                                                                                             >  6:38106/1‑148 (MQ=255)
                                     ttcaaCGCATCGGTGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGTCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCt                                                                                                        >  3:160271/4‑148 (MQ=255)
                                        aaCGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAAc                                                                                                    <  8:176676/149‑1 (MQ=255)
                                               cGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATc                                                                                             <  7:263202/149‑1 (MQ=255)
                                                    cGTTGAACCATGCCATGGCCGAGCTCAAGGCCGTCACCGTCATGCTCAGCGTGTCACGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCGGGCACTGATCGCTTAGCTGAACGCGGATCTGGTc                                                                                        <  4:99591/149‑1 (MQ=255)
                                                    cGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTc                                                                                        <  2:199398/149‑1 (MQ=255)
                                                     gTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTc                                                                                        >  1:237540/1‑148 (MQ=255)
                                                           ccAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTg                                                                                           >  2:108984/1‑139 (MQ=255)
                                                           ccAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTg                                                                                           <  1:108984/139‑1 (MQ=255)
                                                            cAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAAtt                                                                                >  7:134/1‑149 (MQ=255)
                                                                ccATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCAc                                                                             >  7:223599/1‑148 (MQ=255)
                                                                 cATGGCCGAGCTCCAGGCCGTCAACCCCCTGGTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATc                                                                                                                <  4:151050/112‑1 (MQ=255)
                                                                   tGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACcaac                                                                         <  3:41858/149‑1 (MQ=255)
                                                                     gCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACcaac                                                                         >  7:178269/1‑147 (MQ=255)
                                                                      ccGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAg                                                                      >  4:65280/1‑149 (MQ=255)
                                                                         aGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGtt                                                                   >  6:174428/1‑149 (MQ=255)
                                                                          gCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGttt                                                                  <  8:178269/149‑1 (MQ=255)
                                                                             cAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCACCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTAc                                                                >  3:211484/1‑148 (MQ=255)
                                                                             cAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCACCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTAc                                                                >  3:149953/1‑148 (MQ=255)
                                                                                       aaCCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGt                                                      <  2:237540/148‑1 (MQ=255)
                                                                                              tGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTgg                                              >  8:272430/1‑149 (MQ=255)
                                                                                                cTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGcc                                            >  1:20454/1‑149 (MQ=255)
                                                                                                   aGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGt                                                                                         >  6:173030/1‑101 (MQ=255)
                                                                                                   aGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGt                                                                                         <  5:173030/101‑1 (MQ=255)
                                                                                                          cccGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCACCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTg                                  <  4:211484/149‑1 (MQ=255)
                                                                                                          cccGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCACCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTg                                  <  4:149953/149‑1 (MQ=255)
                                                                                                            cGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGAt                                >  7:84017/1‑149 (MQ=255)
                                                                                                                 gCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGATGCTgg                           <  8:223599/149‑1 (MQ=255)
                                                                                                                       cgcgaCGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGGCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGATGCTGGGCCaga                     <  5:38106/146‑1 (MQ=255)
                                                                                                                           aCGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGATGCTGGGCCAGAGCCt                 <  4:160271/149‑1 (MQ=255)
                                                                                                                                ccaccaGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGATGCTGGGCCAGAGCCTGcgcg            >  3:25596/1‑149 (MQ=255)
                                                                                                                                   ccagcaccagATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGATGCTGGGCCAGAGCCTGCGCGAgg         <  2:20454/149‑1 (MQ=255)
                                                                                                                                          cagATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGATGCTGGGCCAGAGCCTGCGCGAGGCCCACAg  >  7:256928/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GCGTTCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTGATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATGGCGTGATGCTGGGCCAGAGCCTGCGCGAGGCCCACAG  >  NC_002947/3051448‑3051734

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