Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,816,917 | C→T | 20.0% | intergenic (+103/+26) | PP_5092 → / ← pilT | NLP/P60 protein/twitching motility protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,816,917 | 0 | C | T | 20.0% | 68.4 / 11.0 | 35 | intergenic (+103/+26) | PP_5092/pilT | NLP/P60 protein/twitching motility protein |
Reads supporting (aligned to +/- strand): ref base C (8/20); new base T (3/4); total (11/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.52e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGCCCCTCC > NC_002947/5816769‑5817060 | gCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgcgc < 8:329105/149‑1 (MQ=255) gCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgcgc < 7:295000/149‑1 (MQ=255) tgctgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCGCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCGCCCATACCAGTGGGCGCGGGCGCgcccc < 6:308061‑M2/149‑6 (MQ=255) gAAAGAGCAGCAGCTGAGTCAGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAgg > 5:147032/1‑149 (MQ=255) gagCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCt < 7:219263/149‑1 (MQ=255) agcagcTGAGGCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACCGGTAGCGCGGCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGc < 5:234231/149‑1 (MQ=255) gCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGc < 3:303599/148‑1 (MQ=255) aaGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACCCCCCCGAACACACGTATCCCGTGGCTCCACCCCGCAGCGCACAAAACGGTGGGACCGGGCCCGCCCCCGCAACGGGCGGTGCAGGTACGACaccgcgtt > 4:324673/1‑142 (MQ=255) tccGTTACGCTGATCGGCGCTAGCCGGTCCCGGCCCCTTCGGGGCCACGCCCGCTCCCACAGGTCGCGCGTCGCTCTCGCGCGCAGGTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTg < 3:44113/147‑1 (MQ=38) gCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGa < 6:147032/149‑1 (MQ=255) cGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGAt < 5:308061/149‑1 (MQ=255) aTCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGGGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCg < 3:98624/148‑1 (MQ=255) ccGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCaa < 4:116327/149‑1 (MQ=255) cGGCCCCTGCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGtac < 1:283504/107‑1 (MQ=255) gCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCa < 7:361776/136‑1 (MQ=255) gCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGc < 7:216542/149‑1 (MQ=255) gggCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTg > 7:144898/1‑107 (MQ=255) gggCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCAAATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTg < 8:144898/107‑1 (MQ=255) gCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaagggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcaccc < 6:38614‑M2/149‑83 (MQ=255) gCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACCTACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGCTCCACCCTCGCTCCGTCAGCCTCCGCAACGCCCGCGCCAGCGGCAGc > 3:133402/1‑149 (MQ=38) cccGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttg < 2:358654‑M2/149‑88 (MQ=255) tCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttggccac < 2:298385‑M2/149‑93 (MQ=255) ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttggccacca < 3:336307‑M2/149‑95 (MQ=255) ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctg > 5:99054‑M2/1‑55 (MQ=255) ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctg < 6:99054‑M2/117‑63 (MQ=255) ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcag < 7:209433‑M2/132‑78 (MQ=255) ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcag > 8:209433‑M2/1‑55 (MQ=255) cacaGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgctt > 1:219933‑M2/1‑54 (MQ=255) cacaGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgctt < 2:219933‑M2/101‑48 (MQ=255) aGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttggccaccacgta > 4:336307‑M2/1‑51 (MQ=255) gTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttggccaccacgtagt < 2:6302‑M2/149‑101 (MQ=255) gcgctcGCGGGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCTGCCCCGCCCTGc < 6:11608/146‑1 (MQ=255) gCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaactcttggccagactcagctgctgctctttcggcagcacccgcttggccaccacgtagtgcttctgcca > 1:358654‑M2/1‑39 (MQ=255) cgcAGCTGACATAGCTGGGGGAGCGGGCGCGGCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGAGCAACCCTCGCTCCCTCTGCCTCCGCAACGCCCCCTCCTTCGTCAGCATCCCCTCCTCCGCCCCGCCTTGCATCACTGaa < 6:276210/149‑1 (MQ=255) acaTACCTGTGGTAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACaa < 3:324673/146‑1 (MQ=255) gAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcactccgcttggccaccacgtagtgcttctgccagtaaggcttgctcaaggtgtcgatactc > 1:298385‑M2/1‑10 (MQ=255) gggCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCAGCGTCAGCATCCCCTCCGCCGCCGCGGCCGGGATCACTGAAGACCACTGCGCCCACCGCCCCTcc > 4:165812/1‑149 (MQ=255) | GCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGCCCCTCC > NC_002947/5816769‑5817060 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |