Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,958,261 | T→C | intergenic (‑830/‑31) | ppiC ← / → rep | peptidyl‑prolyl cis‑trans isomerase C/rep helicase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,958,261 | 0 | T | C | 90.5% | 48.7 / ‑1.2 | 21 | intergenic (‑830/‑31) | ppiC/rep | peptidyl‑prolyl cis‑trans isomerase C/rep helicase |
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (7/12); total (8/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.61e-01 |
TTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑CT‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTAC > NC_000913/3958171‑3958335 | ttCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑t‑‑‑c‑‑‑‑ct‑cc < 2:2731488/93‑1 (MQ=255) cGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑CT‑cccc < 2:1644757/93‑1 (MQ=255) ggTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑c‑‑‑‑cc‑ccccc < 2:955977/93‑1 (MQ=255) ggTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑c‑‑‑‑cc‑ccccc > 2:1094200/1‑93 (MQ=255) gTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑T‑‑‑‑CC‑cccgcc > 1:1900866/1‑93 (MQ=255) gTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑cc‑cccccc > 1:1877640/1‑93 (MQ=255) tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑cc‑ccccccc > 1:3323309/1‑92 (MQ=255) cttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑cc‑cccccccc > 1:783096/1‑91 (MQ=255) ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑cc‑cccccccgt > 2:1124731/1‑90 (MQ=255) ctGATTTAAACATCCGCAGCCAACCGGTTAGCTGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑CT‑CCCCCCGTTCGAAGAtt < 2:2364783/93‑1 (MQ=255) ctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TCC‑C‑‑‑‑CC‑CCCCCCGTTCGAAGa < 2:4618737/93‑1 (MQ=255) aTTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TC‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGa < 1:875961/93‑1 (MQ=255) tttAAACACCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TC‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAg < 1:3131098/93‑1 (MQ=255) ttAAACAACCCCCGGCAACCGGTTTGGGGGTTAACCCGGGTCCAATTCAAAAGGGATTTT‑GCTACAAACCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGa < 2:2878456/93‑1 (MQ=255) aaaCACCCCCAGCAACCCGGTTAGGGGTTTAACCCGGGTCAAATTAAAATGGAATTTG‑GCTACAACCCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGc < 2:2813339/93‑1 (MQ=255) aaCATCCGCAGCCAACCGGGTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCT‑GCTACAAACCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCa < 2:4366784/93‑1 (MQ=255) cgtCCGCCGCCAACCTGTTAGCGGCTTACACGCGGTCCCATTCAAATGCGATTCT‑GCTACAAACCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAAt < 2:148386/91‑1 (MQ=255) aTCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TCCCC‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATa < 1:4469294/93‑1 (MQ=255) aTCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TC‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATacc > 1:4131163/1‑92 (MQ=255) ccGCAGCCAACCGGTTAGGGGCTTACACGGGGGCCCATTCAAATGCGATTCT‑GCTACAACCCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATaca < 2:2108897/93‑1 (MQ=255) cAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TC‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTAt < 2:460553/93‑1 (MQ=255) cgCCACCCGGTTAGCGGCTTACCGGGGGTCCAATTAAAATG‑GGTTTT‑GCTTCAAACCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTAt < 2:1346246/92‑1 (MQ=255) cgCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGGGATTCTGGCAACAAC‑‑C‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTAt < 2:4716459/92‑1 (MQ=255) cgCCAACCCGGTAGCGGGTTACACGCGGTCCCATTCAAATGCCATTTT‑CCTACAAACCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTa < 2:3418945/92‑1 (MQ=255) aCCGGTTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TCC‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTc < 1:3221855/93‑1 (MQ=255) aCCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑CCCCACC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑Cg > 2:579885/1‑93 (MQ=255) aCCGGTTAGCGCGTTAACCGCGGTCACATCGAAAGGCGATTTC‑GCCAAAAACCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGt < 2:1232312/93‑1 (MQ=255) ggTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TCC‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCAATG‑CGTCTaa > 1:792155/1‑93 (MQ=255) ggTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TC‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACAACTATG‑CGTATaaa > 2:3747202/1‑93 (MQ=255) gTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TC‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAAc < 1:2704255/93‑1 (MQ=255) gTTAACGGCTTACCCCCCGTACAATTAAAATTCGATTTG‑GCTACAACC‑C‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTaaa < 2:1545963/93‑1 (MQ=255) tAGCGGCTTACACGCGGTAACAATCAAAATCCATTCT‑GCTACAA‑‑‑TCC‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAAcc < 2:622150/93‑1 (MQ=255) cgccTTACACGCGGACACATTAAAATGGGATTCT‑GCAACAAAC‑C‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCg < 1:4368907/90‑1 (MQ=255) ggCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TC‑‑C‑‑‑‑CC‑CCCCCCGTTTGAAGATTGAGCAATACACCTATG‑CGTTTCAACCCCgtc > 2:288865/1‑91 (MQ=255) ggCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑CT‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGcc > 2:2009112/1‑93 (MQ=255) gCTTACCCCGGGTCCCATTTAAATGCCATTCT‑GGCACAACC‑C‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGc < 2:2217266/93‑1 (MQ=255) ttACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TC‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCcaa < 1:1055276/93‑1 (MQ=255) acacGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑TCC‑C‑‑‑‑CC‑CCCCCCCTTCGAAGATTGAGCAATACCCCTATG‑CGCCTAAACCCCCGCcaac > 1:4077643/1‑93 (MQ=255) cccccgGTCACATTCAAAAGCGATTCT‑GCTAAAAACCC‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCcaac < 2:2434500/89‑1 (MQ=255) acGCGGCCACATTCAAAGGCGATGCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑CT‑ACCCTCGTTCGAAGATTGAGCAATACACCTATG‑CGGCTAAGCCCCTGCcaacaaca < 2:2223663/93‑1 (MQ=255) ggTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑CCCCCCCCCGGTTGGAGATTTAGCAAAAACACTATGCCGGTTAAACCCCGCCaaacaaaaaa > 2:1275527/1‑92 (MQ=255) atCAAATGCGATTCT‑GCTACAA‑‑‑TCC‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCCAACAACAAGCTGTCGa < 2:1250136/92‑1 (MQ=255) tCAAATGCGATTTT‑GCTACAACC‑C‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCCAACAACAAGCTGTCGaa < 2:1139802/93‑1 (MQ=255) cAAAGGCGTTTCT‑GCTACAA‑‑‑C‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTc < 1:2580970/93‑1 (MQ=255) gCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTAc < 1:4558233/93‑1 (MQ=255) gCGATTCT‑GCTACAA‑‑‑G‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTAc < 1:4558232/93‑1 (MQ=255) cGATTCT‑GCTACAACC‑C‑‑‑C‑‑‑‑CC‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTa < 2:155745/93‑1 (MQ=255) | TTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCT‑GCTACAA‑‑‑T‑‑‑C‑‑‑‑CT‑CCCCCCGTTCGAAGATTGAGCAATACACCTATG‑CGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTAC > NC_000913/3958171‑3958335 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |