Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,190,651 T→C 20.4% L56P (CTG→CCG)  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,6510TC20.4% 35.6 / 2.1 19L56P (CTG→CCG) xcpXtype II secretion system protein K
Reads supporting (aligned to +/- strand):  ref base T (9/6);  new base C (2/2);  total (11/8)
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

GCGGCCCTTTTGATGGTGATGGTGGTGCTGGCAATGCTGGCGGCGGGCATGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTT  >  NC_002947/1190509‑1190792
                                                                                                                                              |                                                                                                                                             
gCGGCCCTTTTGATGGTGATGGTGGTGCTGGCAATGCTGGCGGCGGGCATGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTAt                                                                                                                                        >  2:264482/1‑150 (MQ=255)
           gatggtgatggtGGTGCTGGCAATGCTGGCGGCGGGCATGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCTGCCCGGCATAGGCCAGGCCGGc                                                                                                                             <  2:460656/150‑1 (MQ=255)
                     gtggtgCTGGCAATGCTGGCGGCGGGCATGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGAc                                                                                                                   <  2:524465/150‑1 (MQ=255)
                                                aTGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCa                                                                                                                                >  2:581847/1‑110 (MQ=255)
                                                aTGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTc                                                                                        >  2:458512/1‑150 (MQ=255)
                                                aTGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGCTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCa                                                                                                                                <  1:581847/110‑1 (MQ=255)
                                                          tggtggAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCt                                                                              >  2:7146/1‑150 (MQ=255)
                                                             tggAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCTGCCCGGCATAGGCCAGGCCGGc                                                                                                                             <  2:590954/100‑1 (MQ=38)
                                                             tggAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCTGCCCGGCATAGGCCAGGCCGGc                                                                                                                             >  1:590954/1‑100 (MQ=38)
                                                                         ggCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGc                                                               <  2:277870/150‑1 (MQ=255)
                                                                                               ccTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTgg                                         >  2:492534/1‑150 (MQ=255)
                                                                                                  gctgcACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCa                                      <  1:264482/150‑1 (MQ=255)
                                                                                                      cACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGgcg                                  >  1:489499/1‑150 (MQ=255)
                                                                                                          aGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGt                              >  2:581673/1‑150 (MQ=255)
                                                                                                             gAGGGCAGGGGCGGGCCGTGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGc                           <  2:66476/150‑1 (MQ=255)
                                                                                                               gTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCgg                         <  1:492534/150‑1 (MQ=255)
                                                                                                                          ggCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCa              >  1:299366/1‑150 (MQ=255)
                                                                                                                                    agggcctgccCGGCATAGGCCAGGCCGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGc    >  1:112197/3‑150 (MQ=255)
                                                                                                                                      ggccggccTGGCCTATGCCGGGCAGGCCCTGCGTAACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCtt  >  1:292906/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GCGGCCCTTTTGATGGTGATGGTGGTGCTGGCAATGCTGGCGGCGGGCATGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTT  >  NC_002947/1190509‑1190792

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