Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,190,826 A→G 12.5% E114E (GAA→GAG)  xcpX → type II secretion system protein K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,190,8260AG12.5% 75.6 / 4.7 32E114E (GAA→GAG) xcpXtype II secretion system protein K
Reads supporting (aligned to +/- strand):  ref base A (18/10);  new base G (2/2);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

GGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAGGGCGCCGAGAGCGCCCAG  >  NC_002947/1190708‑1190952
                                                                                                                      |                                                                                                                              
ggCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGtt                                                                                                               <  1:226988/136‑1 (MQ=255)
  cAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGc                                                                                                             >  1:235094/1‑136 (MQ=255)
       cTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGttac                                                                                                        <  1:123652/136‑1 (MQ=255)
            gcgGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTtacttact                                                                                                   <  2:41275/136‑1 (MQ=255)
                   gCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAGCGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAAt                                                                                            >  2:159965/1‑136 (MQ=255)
                   gCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAGCGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAAt                                                                                            >  2:158940/1‑136 (MQ=255)
                        tGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACAt                                                                                       <  2:55312/136‑1 (MQ=255)
                             tACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATggcgg                                                                                  <  1:127711/136‑1 (MQ=255)
                                   tttGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAGCGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAAc                                                                            <  1:158940/136‑1 (MQ=255)
                                   tttGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAGCGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAAc                                                                            <  1:159965/136‑1 (MQ=255)
                                       gcgcTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCt                                                                        >  1:32202/1‑136 (MQ=255)
                                         gcTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTgg                                                                      >  2:66909/1‑136 (MQ=255)
                                                          gCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGAc                                                     >  1:239906/1‑136 (MQ=255)
                                                          gCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGAc                                                     >  2:235276/1‑136 (MQ=255)
                                                                     aCCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGc                                          >  1:49435/1‑136 (MQ=255)
                                                                       cAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCg                                        >  2:50519/1‑136 (MQ=255)
                                                                             acCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTg                                  <  1:39855/136‑1 (MQ=255)
                                                                             acCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTg                                  <  1:253406/136‑1 (MQ=255)
                                                                               cTGATTAAAGGTAAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGAt                                <  2:32202/136‑1 (MQ=255)
                                                                                  cTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGcc                             >  2:222858/1‑136 (MQ=255)
                                                                                   ttCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGccc                            >  2:57718/1‑136 (MQ=255)
                                                                                    tCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCg                           >  2:227457/1‑136 (MQ=255)
                                                                                     cAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCATTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGc                          >  1:200899/1‑136 (MQ=255)
                                                                                       aCGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCt                        >  1:204316/1‑136 (MQ=255)
                                                                                        cGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTg                       <  2:220145/136‑1 (MQ=255)
                                                                                            aaTGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAgg                   >  1:173220/1‑136 (MQ=255)
                                                                                               gCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAGGgcg                >  2:209830/1‑136 (MQ=255)
                                                                                               gCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAGGgcg                >  2:159689/1‑136 (MQ=255)
                                                                                               gCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAGGgcg                >  2:143694/1‑136 (MQ=255)
                                                                                               gCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAGGgcg                >  1:44502/1‑136 (MQ=255)
                                                                                                   tGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAGGGCGCCga            <  1:220595/136‑1 (MQ=255)
                                                                                                             cGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAGGGCGCCGAGAGCGCCCAg  >  1:63172/1‑136 (MQ=255)
                                                                                                                      |                                                                                                                              
GGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGCATTGCTCGGCGCCGACGGCGAACGCGCCGAGCGCCAGTTGCGTTACTTACTGGGCAATGACATGGCGGCCGAACGCCTGGTCGATGCCCTGGCCGACTGGCTTGATGCCGACAGTGATGCCCGCCTGCAGGGCGCCGAGAGCGCCCAG  >  NC_002947/1190708‑1190952

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