Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,816,899 | C→T | 10.5% | intergenic (+29/+100) | 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,899 | 0 | C | T | 10.5% | 148.6 / 8.4 | 57 | intergenic (+29/+100) | PP_5092/pilT | NLP/P60 protein/twitching motility protein |
Reads supporting (aligned to +/- strand): ref base C (25/26); new base T (4/2); total (29/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.70e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01 |
ACGCCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCC > NC_002947/5816753‑5817046 | acgcCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTc > 2:1407333/1‑150 (MQ=255) gcCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATAGGCGCTATCCGGTACCGGCCCCTTcgc > 2:770972/1‑150 (MQ=255) cGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCgg < 1:145803/150‑1 (MQ=255) cGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACCACGTGGTCAGC‑GGCGTTGTTCTGCCAACAGAGCAGCCGCTGAGTACGGCAAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCgg < 1:1032966/150‑1 (MQ=255) gttcgggaggggcggttggcgcagttgtcttcagtgatgcagggcggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGACGCTtcgcggcccc < 2:1371059‑M1/17‑5 (MQ=255) aaGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGc > 2:838878/1‑150 (MQ=255) aaGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGc < 2:534134/150‑1 (MQ=255) aaGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGc > 1:864218/1‑150 (MQ=255) gggTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGc > 2:72284/1‑150 (MQ=255) aGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCca < 1:838878/150‑1 (MQ=255) aGTAGCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGcgaaag < 2:512939/150‑6 (MQ=255) tATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTcc > 2:982660/1‑146 (MQ=255) tATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTcc < 1:982660/146‑1 (MQ=255) tATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCcaca > 1:688967/1‑150 (MQ=255) cGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGt < 1:72284/150‑1 (MQ=255) cacCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAgc < 2:173575/150‑1 (MQ=255) cTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGcgcg > 1:458995/1‑150 (MQ=255) aaGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGctct < 2:688967/150‑1 (MQ=255) ttACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTcgcgc > 1:259083/1‑150 (MQ=255) atgcagggcggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCgcgcgc < 2:249379‑M1/53‑1 (MQ=255) cggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTcaca > 2:372815‑M1/90‑150 (MQ=255) cTACGTGGTGGTCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCt < 2:438318/150‑1 (MQ=255) cTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCt > 1:280247/1‑150 (MQ=255) cggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg > 1:249379‑M1/82‑150 (MQ=255) gCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgt > 2:1423363/1‑142 (MQ=255) gCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCcgt < 1:1423363/142‑1 (MQ=255) gatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCgg > 2:1347839‑M1/75‑150 (MQ=255) cgaaGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCggg < 2:1018270/148‑1 (MQ=255) cAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCgg < 2:681961/85‑1 (MQ=255) cAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCgg > 1:681961/1‑85 (MQ=255) cAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgc > 2:465096/1‑150 (MQ=255) aaGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg > 1:45655/1‑140 (MQ=255) aaGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg < 2:45655/140‑1 (MQ=255) tgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg < 1:394207‑M1/70‑1 (MQ=255) tgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg > 2:394207‑M1/70‑139 (MQ=255) tgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGaa > 2:832039/1‑150 (MQ=255) gcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGaaa < 1:563450/150‑1 (MQ=255) ccGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAg < 2:280247/150‑1 (MQ=255) cGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAgg > 1:523756/1‑150 (MQ=255) gAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAggg < 2:365810/150‑1 (MQ=255) gcTTAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAAACATCCCGGCGATCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGta > 1:602389/1‑150 (MQ=38) tGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGc < 1:465096/150‑1 (MQ=255) agggttgatcaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtg < 1:1007259‑M1/67‑1 (MQ=255) agggttgatcaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtg > 2:1007259‑M1/33‑99 (MQ=255) agggttgatcaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg > 1:1395030‑M1/33‑102 (MQ=255) agggttgatcaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg < 2:1395030‑M1/70‑1 (MQ=255) aaGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg < 1:747411/102‑1 (MQ=255) aaGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg > 2:747411/1‑102 (MQ=255) gTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAgg < 1:978194/120‑1 (MQ=255) gTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAgg > 2:978194/1‑120 (MQ=255) gCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCt > 1:1234842/1‑150 (MQ=255) cTGATCGGCGCTATCCGGTACCGGCCCCTGCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTc < 2:676692/150‑1 (MQ=255) tCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCt > 1:578376/1‑150 (MQ=255) ccGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAAc < 2:458995/150‑1 (MQ=255) GCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCCGCCTCCGCAACGcc < 2:66690‑M1/150‑1 (MQ=255) GCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGcc > 1:832039‑M1/1‑150 (MQ=255) CCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGcccc > 2:602389‑M1/1‑150 (MQ=255) | ACGCCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCG‑GGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCC > NC_002947/5816753‑5817046 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |