Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,957,560 | +C | intergenic (‑129/‑732) | 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,957,558 | 1 | . | C | 94.2% | 175.0 / 6.5 | 69 | intergenic (‑127/‑734) | ppiC/rep | peptidyl‑prolyl cis‑trans isomerase C/rep helicase |
Reads supporting (aligned to +/- strand): ref base . (4/0); new base C (33/32); total (37/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.18e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
TAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑T‑‑GC‑TCGCA‑GG‑T‑CCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTC > NC_000913/3957494‑3957628 | tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAgg < 2:27679/76‑1 (MQ=255) gaaTTGTCGGTGATGTGGTGGGTGTTACATATATGGCATAAGCGCGGAAAATCC‑GGGGT‑TCGCA‑GG‑TCCCTCAAgg < 2:59471/75‑1 (MQ=255) aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAgggg < 1:4104952/76‑1 (MQ=255) aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAgggg < 2:5318377/76‑1 (MQ=255) gtcagtcaGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGga‑ga > 1:3954680/1‑76 (MQ=255) caggcaggGAGTGGCTGGTACATATAGGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑Gaaa < 2:5303664/76‑1 (MQ=255) agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑Gaaaa < 1:4099186/76‑1 (MQ=255) agtcagGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCC‑GG‑GCTTCGCAGGG‑C‑CCTCCAGGGGA‑Gaaaa > 1:845718/1‑76 (MQ=255) gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAAc < 2:2574174/76‑1 (MQ=255) tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCCCA‑GG‑TCCCCCCAGGGGA‑GAAAACt > 2:256689/1‑76 (MQ=255) tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑GG‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACt > 1:2119854/1‑76 (MQ=255) agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTa < 2:1938691/76‑1 (MQ=255) agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTa < 1:3345907/76‑1 (MQ=255) gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAg < 1:1653847/76‑1 (MQ=255) gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAg < 1:4046882/76‑1 (MQ=255) gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑GG‑GCTTCGCAGGG‑C‑CCTCAAGGGGA‑GAAAAATTAg > 1:189100/1‑76 (MQ=255) ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAgg < 2:4395766/76‑1 (MQ=255) gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTa > 1:4297471/1‑73 (MQ=255) gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTa < 2:4297470/73‑1 (MQ=255) gAGTGGCTGGGACAGATATGGCAAAAGCCCGGAAAATCC‑CG‑GCTTCGCA‑GGCC‑CCTCCAGGGGA‑GAAAACTTAggg > 2:881598/1‑76 (MQ=255) aGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCCAGGGGGGGAAAACTAAggc > 1:850930/1‑75 (MQ=255) aGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCC‑GGGGT‑TCGCAGGG‑C‑CCTCCAGGGGA‑GAAAAATTAgggg > 2:3322551/1‑75 (MQ=255) tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCct > 2:5325025/1‑76 (MQ=255) cTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTa < 2:1672678/76‑1 (MQ=255) cTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGAGCCTCTa > 1:2017151/1‑76 (MQ=255) tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCACGGGCA‑GAAAACTTAGGGCCTCTag > 2:5584087/1‑75 (MQ=255) tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTAt > 1:3977248/1‑76 (MQ=255) tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTAt > 1:5084326/1‑76 (MQ=255) tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTAt < 2:2501805/76‑1 (MQ=255) tGGCACATAGAGGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCGAt > 1:683832/1‑76 (MQ=255) ccATATCTGGCAAACGCCCGGGAACTCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTATGGCCTCTATGGGc > 1:3771546/2‑76 (MQ=255) aCATATATGGCAAAAGCCCGGGAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGc > 1:3739029/1‑76 (MQ=255) aCATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGc > 2:2977680/1‑76 (MQ=255) aCATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGc > 1:1973741/1‑76 (MQ=255) cATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCt < 1:1325010/76‑1 (MQ=255) atatatGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACt > 1:5355472/1‑58 (MQ=255) atatatGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACt < 2:5355471/58‑1 (MQ=255) atatatGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTa < 2:1624104/76‑1 (MQ=255) atatGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCtcca < 2:4230286/76‑4 (MQ=255) atGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACaaa < 1:3927636/76‑1 (MQ=255) tGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAg > 2:4025058/1‑76 (MQ=255) ggCAAAATCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGt < 2:4644780/76‑1 (MQ=255) ggCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACCAAGt > 2:5300317/1‑76 (MQ=255) gCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTc < 1:4826092/76‑1 (MQ=255) cAAAAGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCa > 1:3013893/1‑76 (MQ=255) aaaaGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt > 1:5347528/1‑76 (MQ=255) aaaGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt > 2:3810144/1‑76 (MQ=255) aaGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa > 1:1262481/1‑76 (MQ=255) aGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt > 2:4598516/1‑76 (MQ=255) aGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt < 2:147653/76‑1 (MQ=255) aGCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt > 2:2443414/1‑76 (MQ=255) gCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc < 1:4189929/76‑1 (MQ=255) gCCCGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc > 1:5625757/1‑76 (MQ=255) cccGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg < 2:4927619/76‑1 (MQ=255) ccGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg > 2:934393/1‑76 (MQ=255) cGGAAAATCCGG‑‑GCTTCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg < 2:5596078/76‑1 (MQ=255) cttCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta < 2:4041776/74‑1 (MQ=255) ttCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat > 1:675557/2‑76 (MQ=255) ttCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat > 2:3233589/2‑76 (MQ=255) ttCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat > 1:2473354/2‑76 (MQ=255) ttCGCA‑GG‑ACCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCAACAAAGTCACTATCGCGGCTGGttatttcat > 1:5314072/2‑76 (MQ=255) tCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCa > 2:5066697/1‑53 (MQ=255) tCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCa < 1:5066698/53‑1 (MQ=255) tCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt < 1:256689/76‑1 (MQ=255) tCGCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt < 2:890388/76‑1 (MQ=255) gCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt > 1:5280272/1‑76 (MQ=255) gCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt < 1:3163554/76‑1 (MQ=255) gCA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt < 2:396928/76‑1 (MQ=255) cA‑GG‑TCCCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc > 2:4211582/1‑76 (MQ=255) | TAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑T‑‑GC‑TCGCA‑GG‑T‑CCTCAAGGGGA‑GAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTC > NC_000913/3957494‑3957628 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |