Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,190,658 T→G 19.1% Y58* (TAT→TAG 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,6580TG19.1% 38.7 / 4.1 20Y58* (TAT→TAGxcpXtype II secretion system protein K
Reads supporting (aligned to +/- strand):  ref base T (9/7);  new base G (2/2);  total (11/9)
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

GCGGCCCTTTTGATGGTGATGGTGGTGCTGGCAATGCTGGCGGCGGGCATGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGC  >  NC_002947/1190509‑1190804
                                                                                                                                                     |                                                                                                                                                  
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)
                                                                                                                                                 cctATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATg   <  2:489499/150‑1 (MQ=255)
                                                                                                                                                  ctATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGc  >  2:426177/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGGCCCTTTTGATGGTGATGGTGGTGCTGGCAATGCTGGCGGCGGGCATGGCCTGGCTGGTGGAAGATGGCCGGCGTCAGGTCGACGAAGTGCGCCTGCTGCACCAGCGAGTGCAGGTGCGGGCCATGGAGCAGGCCGGCCTGGCCTATGCCGGGCAGGCCCTGCGTGACCCCGCCTGGCGGCTCAGTCCGCTGTTCTGGCAAGCCTTGCGCGGGCAGCCGCTGAACTACGACTTTGGCGCTGGCCAGGCGCAGTTGCGGGTACGCGACCAGCACACCTGCTTCAACGTCAATGC  >  NC_002947/1190509‑1190804

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: