Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,793,795 | G→T | 12.9% | intergenic (+61/+66) | PP_2445 → / ← PP_2446 | membrane protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,793,795 | 0 | G | T | 12.9% | 83.5 / 3.5 | 31 | intergenic (+61/+66) | PP_2445/PP_2446 | membrane protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (12/15); new base T (2/2); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01 |
CATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGAC > NC_002947/2793653‑2793943 | cATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTg‑g < 2:1255968/150‑1 (MQ=255) gTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTg‑g < 2:712875/144‑1 (MQ=255) gTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTg‑g > 1:712875/1‑144 (MQ=255) ctcGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATGCGCG‑TGGCACTG‑GCAACGCGaa < 1:807164/150‑1 (MQ=255) tcGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCAc > 1:296126/1‑137 (MQ=255) tcGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCAc < 2:296126/137‑1 (MQ=255) gACGGCGATCGAGGTCCTGTGCATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TTGCCAGT‑GCCACGCGAATAc < 2:1036828/150‑1 (MQ=37) gTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGc < 1:1020009/150‑1 (MQ=255) tCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCTGACAAGCCCGCTGCACAGGTATTCGCG‑TTGCCAGT‑GCCACGCGAATACCTGTGCAGCGGGCt > 1:389024/1‑150 (MQ=25) ccTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCtt < 2:1069583/150‑1 (MQ=255) cTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAAc > 2:423245/1‑117 (MQ=255) cTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAAc < 1:423245/117‑1 (MQ=255) gTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑Tgg < 1:145734/100‑1 (MQ=255) gTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑Tgg > 2:145734/1‑100 (MQ=255) cgcgTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGaaa < 1:734402/150‑1 (MQ=255) tttGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTG‑CC‑ACGCGAATACCtgt > 1:1044562/1‑112 (MQ=21) tttGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TTGC‑CAGTGCCACGCGAATACCtgt < 2:1044562/112‑1 (MQ=21) gAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAg < 2:1272709/150‑1 (MQ=255) gTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCtg < 1:364750/82‑1 (MQ=255) gTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCtg > 2:364750/1‑82 (MQ=255) gcaagaGGGGCGCGAGGCGCCCCTTTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAg < 1:632091/118‑1 (MQ=37) aaaGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TTGC‑CAGTGCCACGCGAATACCt > 1:942150/1‑77 (MQ=11) cccTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATgccgc > 1:471799/1‑126 (MQ=255) cccTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATgccgc < 2:471799/126‑1 (MQ=255) aaGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCGACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCgg > 1:300027/1‑144 (MQ=255) aaGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCGACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCgg < 2:300027/144‑1 (MQ=255) cccGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTTTCTGTGTCTTGGCGGGTGCCTGCg > 2:298541/1‑150 (MQ=255) acaGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGGTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAgg < 2:608261/150‑1 (MQ=255) acaGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAgg > 1:384822/1‑150 (MQ=255) gTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAggtggt > 2:307308/1‑150 (MQ=255) gcg‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCag > 1:18390/1‑149 (MQ=255) gcg‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCAt > 1:220741/1‑150 (MQ=255) tGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGAc < 2:18390/150‑1 (MQ=255) | CATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCG‑TGGCACTG‑GCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGAC > NC_002947/2793653‑2793943 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |