Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 319,217 A→T 14.8% I176N (ATC→AAC)  PP_0264 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947319,2170AT14.8% 63.5 / 5.6 27I176N (ATC→AAC) PP_0264sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base T (2/2);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

AGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGAAGTTTCGGCCGACGAAGCTGTCGTGC  >  NC_002947/319070‑319366
                                                                                                                                                   |                                                                                                                                                     
aGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATc                                                                                                                                                     <  1:347811/150‑1 (MQ=255)
        ggTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGTTCCGCTcgc                                                                                                                                              >  1:370907/1‑149 (MQ=255)
        ggTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGTTCCGCTcgc                                                                                                                                              >  1:371071/1‑149 (MQ=255)
        ggTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGTTCCGCTcgc                                                                                                                                              <  2:371071/149‑1 (MQ=255)
        ggTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGTTCCGCTcgc                                                                                                                                              <  2:370907/149‑1 (MQ=255)
         gTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCgg                                                                                                                                            <  1:22312/150‑1 (MQ=255)
                        cAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCcgcg                                                                                                                             >  2:346006/1‑150 (MQ=255)
                                    gAAGTCAGGATGACCACGCGATTGAGATCGGTCAGGTGCACGTTCTGCGGCGTACGGCTTCACAGGGTGTCGCTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACAATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAg                                                                                                                 <  2:352904/150‑1 (MQ=255)
                                                aCCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGcc                                                                                                     <  2:398106/150‑1 (MQ=255)
                                                                                   gggCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTg                                                                  >  1:284643/1‑150 (MQ=255)
                                                                                    ggCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTgg                                                                 >  2:440550/1‑150 (MQ=255)
                                                                                               ccACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGc                                                      <  1:346006/150‑1 (MQ=255)
                                                                                                cacaGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCt                                                     >  2:428675/1‑150 (MQ=255)
                                                                                                        ttCGGTGTGGTCGAGATCGACTGTGACCGCCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGCGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGaa                                             <  2:207652/150‑1 (MQ=255)
                                                                                                             tgtgGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTg                                               <  1:488638/143‑1 (MQ=255)
                                                                                                             tgtgGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTg                                               >  2:488638/1‑143 (MQ=255)
                                                                                                                 gTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTcgc                                                                                                                                              >  1:51840/1‑44 (MQ=255)
                                                                                                                 gTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTcgc                                                                                                                                              <  2:51840/44‑1 (MQ=255)
                                                                                                                          gACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGaa                           <  2:180721/150‑1 (MQ=255)
                                                                                                                          gACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGaa                           <  2:284643/150‑1 (MQ=255)
                                                                                                                          gACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGaa                           <  1:441494/150‑1 (MQ=255)
                                                                                                                          gACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGaa                           >  2:444436/1‑150 (MQ=255)
                                                                                                                                          aCGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGAAGTTTCGGCCGACGAAg           >  1:11894/1‑150 (MQ=255)
                                                                                                                                          aCGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGAAGTTTCGGCCGACGAAg           >  2:281085/1‑150 (MQ=255)
                                                                                                                                            gCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGtt                                                <  1:55188/111‑1 (MQ=255)
                                                                                                                                            gCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGtt                                                >  2:55188/1‑111 (MQ=255)
                                                                                                                                            gCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGAAGTTTCGGCCGACGAAGCt         >  1:292938/1‑150 (MQ=255)
                                                                                                                                               gacgaTCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGCGGGAGCGGCCAAGGTGTCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGAAGTTTCGGCCGACGAAGCTGTc      <  1:428675/150‑1 (MQ=255)
                                                                                                                                                  aaTCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGAAGTTTCGGCCGACGAAGCTGTCGTg   >  2:337172/2‑150 (MQ=255)
                                                                                                                                                   aTCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGAAGTTTCGGCCGACGAAGCTGTCGTGc  <  1:444436/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCACGCCATTGTGATCGGTCAGCAACAGCTGCTCGGGCGTACGGCCCCACAGGGTTTCGGTGTGGTCGAGATCGACTTTGACCACCAGCACGCCGACGATCCGCTCGCGGTCGCGCACGGCCGCGGCAAAGAAGTAGCCGCGCTTGGCCGAGGTGGTGCCCTGGCCGAAGAAGCGGCCAAGGTGGCCGGCCATGGCTTCGTTGAAGTAAGGGCGGAAGGCGAAGTTTCGGCCGACGAAGCTGTCGTGC  >  NC_002947/319070‑319366

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