Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,433,718 | T→C | intergenic (‑44/+273) | pinR ← / ← ynaE | Rac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,433,718 | 0 | T | C | 87.2% | 76.0 / 7.3 | 47 | intergenic (‑44/+273) | pinR/ynaE | Rac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage |
Reads supporting (aligned to +/- strand): ref base T (3/3); new base C (20/21); total (23/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTC > NC_000913/1433622‑1433766 | taggtggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAg > 2:1316859/3‑100 (MQ=16) tCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAAGTCTAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAg < 1:1839971/100‑1 (MQ=255) tggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGaaaa < 1:64845/100‑1 (MQ=17) tggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGaaaa < 2:1042324/100‑1 (MQ=18) ggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaat < 2:200311/100‑1 (MQ=18) gtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaata < 2:1166192/100‑1 (MQ=18) gtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaata > 2:632557/1‑100 (MQ=18) gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa > 2:1263008/1‑100 (MQ=18) aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATCATTa > 2:1722216/1‑99 (MQ=17) aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTAt > 2:826824/1‑100 (MQ=33) aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTAt > 1:434260/1‑100 (MQ=18) cAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGat > 2:1800250/1‑100 (MQ=18) cGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGattat > 1:1397205/1‑100 (MQ=18) tGATATCCGACAGTAAGCAAAAATTCGATACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAtttt < 1:870435/100‑1 (MQ=12) gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAttttt > 2:1772433/1‑100 (MQ=18) atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtg < 2:938673/100‑1 (MQ=32) atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg < 2:247256/100‑1 (MQ=17) atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg < 1:241678/100‑1 (MQ=17) atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg < 2:1637207/100‑1 (MQ=18) tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtgt > 2:591498/1‑100 (MQ=33) tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgt > 1:342466/1‑100 (MQ=17) tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgt < 1:177611/100‑1 (MQ=17) tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgt > 2:711565/1‑100 (MQ=18) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATTTCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg < 1:1201216/100‑1 (MQ=12) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg < 2:305320/100‑1 (MQ=18) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg < 1:849034/100‑1 (MQ=18) ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCa < 2:552724/100‑1 (MQ=18) ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCa < 2:299938/100‑1 (MQ=18) ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAAAAAATAATTATGATTATTTTTGTGCa > 2:1215924/1‑100 (MQ=255) cGACAGTAAGCAAAAATTCGAGACATAGTGCGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGc > 2:1234100/1‑98 (MQ=255) gTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTtg < 1:1459457/100‑1 (MQ=32) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta > 2:161517/1‑100 (MQ=18) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 2:373635/100‑1 (MQ=18) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta > 1:1428576/1‑100 (MQ=18) aaGCAAAAATTCGAGAAATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 2:536498/100‑1 (MQ=12) aaGCAAAAATGCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 2:1133497/100‑1 (MQ=12) aaGAAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 1:296051/100‑1 (MQ=16) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat < 2:1467493/100‑1 (MQ=18) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 2:43264/1‑100 (MQ=18) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGATGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 1:83278/99‑1 (MQ=17) gCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatg > 2:31438/1‑100 (MQ=18) gCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatg > 2:1290776/1‑100 (MQ=18) gCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgt > 1:890923/1‑97 (MQ=18) cAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatgt > 2:637601/1‑100 (MQ=18) aaaaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatgta > 2:1450363/1‑100 (MQ=18) aaaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTac > 1:1145187/1‑100 (MQ=18) cGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTc < 1:727152/100‑1 (MQ=11) | TCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTC > NC_000913/1433622‑1433766 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |