Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,159,143 | C→T | 17.2% | intergenic (+91/‑118) | PP_3661 → / → PP_3662 | membrane protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,159,143 | 0 | C | T | 17.2% | 62.9 / 5.1 | 29 | intergenic (+91/‑118) | PP_3661/PP_3662 | membrane protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (13/11); new base T (3/2); total (16/13) | |||||||||||
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.80e-01 |
TGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTC > NC_002947/4159070‑4159278 | tgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtggt > 6:106902/1‑149 (MQ=255) ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt > 6:87728/1‑149 (MQ=255) gggTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCaa < 5:106902/149‑1 (MQ=255) ggTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCa > 1:316636/1‑147 (MQ=255) gTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACAATTCTCGGTATGGTGGTGTCa > 5:188577/1‑146 (MQ=255) cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg < 6:11332‑M2/69‑1 (MQ=255) cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg > 5:11332‑M2/59‑127 (MQ=255) CCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGtt < 3:6237‑M1/110‑1 (MQ=255) CCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGtt > 4:6237‑M1/1‑110 (MQ=255) CCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACGTTTCTCGGTATGGTGGTGTCAAAGAAGGa < 6:120399‑M1/149‑1 (MQ=255) ccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGAt > 5:304222‑M2/58‑149 (MQ=255) CCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAgg < 8:173743‑M1/147‑1 (MQ=255) gctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGt > 3:119810‑M2/56‑149 (MQ=255) tcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGCCGACTCGGTCTCCAG‑‑‑‑‑T‑‑CAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCg > 7:15341‑M2/54‑148 (MQ=255) CCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTAGGTATGGTGGTGTCAATGacggagttt < 4:305903‑M1/147‑5 (MQ=255) ctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTa > 4:21869‑M2/52‑149 (MQ=255) CTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTa < 3:121914‑M1/149‑1 (MQ=255) TACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTAc > 5:87728‑M1/1‑149 (MQ=255) AGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCg < 8:36449‑M1/149‑1 (MQ=255) gtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAg > 2:214545/1‑148 (MQ=255) ggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGc < 3:21869‑M2/122‑1 (MQ=255) aCGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGc < 7:173743/149‑1 (MQ=255) aaGGCCCCCGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCaa > 8:97279/1‑149 (MQ=255) gCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGcaca > 4:178381/1‑125 (MQ=255) cctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCAt < 6:304222‑M2/133‑1 (MQ=255) cctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCAt > 1:195328‑M2/17‑149 (MQ=255) gtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGc < 1:65892‑M2/122‑1 (MQ=255) gtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGc > 2:65892‑M2/14‑135 (MQ=255) taggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGtc < 2:195328‑M2/136‑1 (MQ=255) | TGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTC > NC_002947/4159070‑4159278 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |