Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,018,109 A→T 10.9% Q167L (CAG→CTG)  PP_3545 → 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_0029474,018,1090AT10.9% 120.7 / 3.1 43Q167L (CAG→CTG) PP_3545two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (19/19);  new base T (2/3);  total (21/22)
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

GACGACGCACCACGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAGCTGCTTGAA  >  NC_002947/4017961‑4018257
                                                                                                                                                    |                                                                                                                                                    
gacgacGCACCACGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAg                                                                                                                                                     >  2:787286/1‑150 (MQ=255)
  cgacgCACCACGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGgc                                                                                                                                                   <  1:126731/150‑1 (MQ=255)
             gCCTCTACGAAATCCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGTTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATgg                                                                                                                                        <  1:119375/150‑1 (MQ=255)
                        aTTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTTCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGc                                                                                                                             <  1:672507/150‑1 (MQ=255)
                                  accCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCagaaga                                                                                                                                           >  1:219539/1‑126 (MQ=255)
                                  accCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCagaaga                                                                                                                                           <  2:219539/126‑1 (MQ=255)
                                    ccGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGt                                                                                                                         >  1:238165/1‑142 (MQ=255)
                                    ccGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGt                                                                                                                         <  2:238165/142‑1 (MQ=255)
                                    ccGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACcggc                                                                                                                 >  1:722830/1‑150 (MQ=255)
                                      gTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACcggcgg                                                                                                               <  1:333982/150‑1 (MQ=255)
                                        ggCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGa                                                                                                             <  1:110934/150‑1 (MQ=255)
                                         gCTCGATGCTCAGCAGCGAATTCAGGGCGGCTCCCTGTTCGCCTATGACATCACGCAGCGCGTGACTGAACAGGAACGCCTGCAGCAAGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGATCGGCCAGTGTACCGTCgggct                                                                                                            >  1:390875/1‑148 (MQ=255)
                                                     gcagcGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGAttt                                                                                                <  2:175181/150‑1 (MQ=255)
                                                            aTTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCaacaac                                                                                         <  1:706162/150‑1 (MQ=255)
                                                                 gggCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACctgct                                                                                    >  2:499472/1‑150 (MQ=255)
                                                                      gcTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTAAGCTGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCa                                                                               >  1:684821/1‑150 (MQ=255)
                                                                        tACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCTGGCGCag                                                                                                                                               >  1:606356/1‑84 (MQ=255)
                                                                        tACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCTGGCGCag                                                                                                                                               <  2:606356/84‑1 (MQ=255)
                                                                          ccTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCAt                                                                           >  2:734694/1‑150 (MQ=255)
                                                                              ttCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCtt                                                                       >  2:757287/1‑150 (MQ=255)
                                                                                 gCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGgc                                                                    >  1:561154/1‑150 (MQ=255)
                                                                                     aTGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGgcgcgc                                                                <  2:804953/150‑1 (MQ=255)
                                                                                      tGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCt                                                               <  2:722830/150‑1 (MQ=255)
                                                                                      tGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGAACCTCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCt                                                               <  1:488768/150‑1 (MQ=255)
                                                                                               gCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGc                                                      >  1:644003/1‑150 (MQ=255)
                                                                                                 aGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCg                                                    <  1:499472/150‑1 (MQ=255)
                                                                                                          ccGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCa                                                                               <  1:403443/114‑1 (MQ=255)
                                                                                                          ccGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCa                                                                               >  2:403443/1‑114 (MQ=255)
                                                                                                             aaCAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGGTGACGCATGATTTCAACAACCTGCTGGCCAGCATCACTGGCGCGATGGACCAGGCCGGGCAACGTCTgg                                        >  1:613016/1‑150 (MQ=255)
                                                                                                             aaCAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTgg                                        >  1:102298/1‑150 (MQ=255)
                                                                                                                aGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCg                                                    >  1:210094/1‑135 (MQ=255)
                                                                                                                aGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCg                                                    <  2:210094/135‑1 (MQ=255)
                                                                                                                     cGCCTGCAGCAGGCCGAAGAAGCACAGCAGGTAGCTCAGGAGATGGAGGCGGTCGGCCAGTTGACAGCCTCAATGGGCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCaa                                <  1:787286/150‑1 (MQ=255)
                                                                                                                      gccgggagcagGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCATTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGACCAGCATCTTTGGCGCGCTGGAACTGGCCGGGCAAAGTCTGGTGTAGCAAc                               <  2:121022/144‑1 (MQ=255)
                                                                                                                             gcagGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGc                        >  2:323647/1‑150 (MQ=255)
                                                                                                                             gcagGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGc                        >  2:651456/1‑150 (MQ=255)
                                                                                                                                         tcttcTGCGCCTGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAAc                                              <  2:412380/112‑1 (MQ=255)
                                                                                                                                         tcttcTGCGCCTGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAAc                                              >  1:412380/5‑116 (MQ=255)
                                                                                                                                            aCTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAgct         >  2:722553/1‑150 (MQ=255)
                                                                                                                                            aCTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAgct         <  1:323647/150‑1 (MQ=255)
                                                                                                                                              tGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAgctgc       <  2:102298/150‑1 (MQ=255)
                                                                                                                                                cgcCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAGCTGCtt     >  1:438877/1‑150 (MQ=255)
                                                                                                                                                   cAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAgctgttgtaa  <  2:613016/150‑7 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACGACGCACCACGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAGCTGCTTGAA  >  NC_002947/4017961‑4018257

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