Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,345,870 C→T 11.4% A225T (GCC→ACC)  PP_2942 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,345,8700CT11.4% 84.0 / 2.9 35A225T (GCC→ACC) PP_2942two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (15/16);  new base T (2/2);  total (17/18)
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.36e-01

TGGCGTGGGTAGTCGCGATAGCGGTCGCGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCTGTGTTTCAGCTCTTCGTAGGCCA  >  NC_002947/3345722‑3346007
                                                                                                                                                     |                                                                                                                                         
tGGCGTGGGTAGTCGCGATAGCGGTCGCGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCg                                                                                                                                          <  2:309396/150‑1 (MQ=255)
 ggCGTGGGTAGTTGCGATAGCGGTCGCGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGt                                                                                                                                         <  1:398954/150‑1 (MQ=255)
 ggCGTGGGTAGTCGCGATAGCGGTCGCGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGt                                                                                                                                         <  2:5168/150‑1 (MQ=255)
 ggCGTGGGTAGTCGCGATAGCGGTCGCGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGt                                                                                                                                         <  2:330000/150‑1 (MQ=255)
 ggCGTGGGTAGTCGCGATAGCGGTCGAGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGt                                                                                                                                         <  2:222856/150‑1 (MQ=255)
       ggTAGTCGCGATAGCGGTCGCGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAGCGATCATGCTATCCGCCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCAAGATCGCGGCTGATCGCCGCGTCGAAGGCGTGGAccc                                                                                                                                   >  2:43947/1‑150 (MQ=255)
           gTCGCGATAGCGGTCGCGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACa                                                                                                                               <  2:206472/150‑1 (MQ=255)
                             gtcgcTATGCATCCACTGGTCCGAAGGGGCAACGATCATGCTCTACGGCAACCTACCAAACATAGTACTGAACAGCGTCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTc                                                                                                             <  1:43947/150‑1 (MQ=255)
                                gcTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGCAGCAACGCACCAGTTCgat                                                                                                          <  2:513664/150‑1 (MQ=255)
                                   gtgCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGc                                                                                                                     >  1:94635/1‑136 (MQ=255)
                                   gtgCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGc                                                                                                                     <  2:94635/136‑1 (MQ=255)
                                                tCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGACTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATATGCTGATTGGTc                                                                                          >  2:284968/1‑150 (MQ=255)
                                                 ccGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTc                                                                                                  <  2:451114/141‑1 (MQ=255)
                                                 ccGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTc                                                                                                  >  1:451114/1‑141 (MQ=255)
                                                              gATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAgg                                                                            >  2:130583/1‑150 (MQ=255)
                                                                 cATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAg                                                                         <  2:229130/150‑1 (MQ=255)
                                                                   tGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCCACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGcc                                                                       >  1:17583/1‑150 (MQ=255)
                                                                       aTCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAg                                                                   >  2:472514/1‑150 (MQ=255)
                                                                               aaCTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGTGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCa                                                           >  1:386158/1‑150 (MQ=255)
                                                                               aaCTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCa                                                           >  2:164090/1‑150 (MQ=255)
                                                                                    aaTTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGGGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATc                                                      >  1:447796/1‑150 (MQ=255)
                                                                                         gCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGaa                                                 <  1:472514/150‑1 (MQ=255)
                                                                                           cTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGTGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAaca                                               <  2:386158/150‑1 (MQ=255)
                                                                                           cTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGa                                                                                                <  1:83009/101‑1 (MQ=255)
                                                                                           cTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGa                                                                                                >  2:83009/1‑101 (MQ=255)
                                                                                            tATTTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGTGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAacac                                              >  2:380174/1‑150 (MQ=255)
                                                                                                                 ccATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGc                          >  1:19535/1‑150 (MQ=255)
                                                                                                                         aGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCTGTGTTTc                  <  2:19535/150‑1 (MQ=255)
                                                                                                                            tCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCTGTGTTTCAGc               >  2:131004/1‑150 (MQ=255)
                                                                                                                                ggCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCt                         >  2:245886/1‑136 (MQ=255)
                                                                                                                                ggCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCt                         <  1:245886/136‑1 (MQ=255)
                                                                                                                                 gCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCTGTGTTTCAGCTCTTc          <  1:130583/150‑1 (MQ=255)
                                                                                                                                    gATCGCCTCGTCCAAGGTGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCTGTGTTTCAGCTCTTCGTa       <  1:380174/150‑1 (MQ=255)
                                                                                                                                      tCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGGTTGGACTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCTGTGTTTCAGCTCTTCGTAgg     >  1:168340/1‑150 (MQ=255)
                                                                                                                                         ccTCGTCCAAGGCGTGGACCCTACAGTAGAAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCTGTGTTTCAGCTCTTCGTAGGCCa  >  1:239822/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
TGGCGTGGGTAGTCGCGATAGCGGTCGCGGTCGCTGTGCATCAACTGGTCCGAAGGGGCAACGATCATGCTATCCGGCCAACTCAATTTGCCTATGTTGTACAGCGCCGCAGCCATGGCCAGATCGCGGCTGATCGCCTCGTCCAAGGCGTGGACCCTACAGTAGCAACGCACCAGTTCGATGATCTGTCGATTGGTCTGCTTGGCTGGAGGCAGCCGCAGATTAGCCAACATCGAGAACACCTCGGTCCCGGTCACATAGCTGTGTTTCAGCTCTTCGTAGGCCA  >  NC_002947/3345722‑3346007

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: