Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,244,269 G→T 16.0% Q154K (CAG→AAG)  PP_1976 ← HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,244,2690GT16.0% 77.2 / 2.1 25Q154K (CAG→AAG) PP_1976HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base G (13/8);  new base T (2/2);  total (15/10)
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.41e-01

CGGCGACATGGCTGTCCGCCGACAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAGGGCGGCACGGGCGCGTTCGACCGCCAGCTCAAAGTCGGCGGCTTCCAG  >  NC_002947/2244134‑2244392
                                                                                                                                       |                                                                                                                           
cGGCGACATGGCTGTCCGCCGACAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTg                                                                                                                             >  1:245603/1‑136 (MQ=255)
   cGACATGGCTGTCCGCCGACAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGcc                                                                                                                          <  1:91953/136‑1 (MQ=255)
     aCATGGCTGTCCGCCGACAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCgg                                                                                                                        <  2:407471/136‑1 (MQ=255)
      cATGGCTGTCCGCCGACAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTTGCCGgc                                                                                                                       >  2:53829/1‑136 (MQ=255)
      cATGGCTGTCCGCCGACAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTTGCCGgc                                                                                                                       >  2:1953/1‑136 (MQ=255)
             gTCCGCCGACAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATc                                                                                                                <  1:265314/136‑1 (MQ=255)
                      cAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGcc                                                                                                                          >  2:196673/1‑117 (MQ=255)
                      cAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGcc                                                                                                                          <  1:196673/117‑1 (MQ=255)
                           tAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCttgttg                                                                                                  >  1:249357/1‑136 (MQ=255)
                               gACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTTGCCGGCAGCAATCAGGCTGCCTTGTTGGGTg                                                                                              <  1:1953/136‑1 (MQ=255)
                               gACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTTGCCGGCAGCAATCAGGCTGCCTTGTTGGGTg                                                                                              <  1:53829/136‑1 (MQ=255)
                               gACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTg                                                                                              <  2:36727/136‑1 (MQ=255)
                                    gCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAggcg                                                                                         <  1:302774/136‑1 (MQ=255)
                                     ctctTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAggcgg                                                                                        >  1:18434/1‑136 (MQ=255)
                                      tctTCAGAGCCAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGAGAACGGTCGCACGTGGCCTGGGTCGCCGCCAACTGCGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAggcggc                                                                                       >  2:249848/1‑136 (MQ=255)
                                        ttCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCtt                                                                                     >  1:522355/1‑136 (MQ=255)
                                              gaCAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCt                                                                               <  2:421927/136‑1 (MQ=255)
                                                      ccGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTa                                                                       <  2:249357/136‑1 (MQ=255)
                                                                           cACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAg                                                  >  2:192770/1‑136 (MQ=255)
                                                                           cACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAg                                                  >  2:460287/1‑136 (MQ=255)
                                                                                      tCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAGGGCGGCACggg                                       >  2:392968/1‑136 (MQ=255)
                                                                                                        aaCGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGGTTTCGGTTTGTGCGTACTCGGATTCGCGGGTGCCCAGGGCGGAACGGGCGCGTTCGACCGCCAGca                     >  2:416759/1‑135 (MQ=255)
                                                                                                                  gggTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAGGGCGGCACGGGCGCGTTCGACCGCCAGCTCAAAGTcggc           >  1:440775/1‑136 (MQ=255)
                                                                                                                  gggTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAGGGCGGCACGGGCGCGTTCGACCGCCAGCTCAAAGTcggc           >  1:242512/1‑136 (MQ=255)
                                                                                                                  gggTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAGGGCGGCACGGGCGCGTTCGACCGCCAGCTCAAAGTcggc           >  1:264831/1‑136 (MQ=255)
                                                                                                                           ccAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAGGGCGGCACGGGCGCGTTCGACCGCCAGCTCAAAGTCGGCGGCTTCCAg  >  1:165528/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                           
CGGCGACATGGCTGTCCGCCGACAGGGTAGTGACCCGCTCTTCAGAGACAAAACCCGGCTTGCGCAGCTTCTCGGCACGGCTCAGGTCCAGCCGTGAACGGTCGAACGTGGCCTGGGTCGCCGCCAACTGGGCCTGGCCGGCAGCAATCAGGCTGCCTTGTTGGGTGAGGCGGCTTTGCGCTTGTGCGTACTCGGCTTCGCGGGTGGCCAGGGCGGCACGGGCGCGTTCGACCGCCAGCTCAAAGTCGGCGGCTTCCAG  >  NC_002947/2244134‑2244392

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: