| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,433,718 | T→C | 100% | 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 | 52.1% | ‑0.4 / 87.8 | 71 | intergenic (‑44/+273) | pinR/ynaE | Rac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage |
| Reads supporting (aligned to +/- strand): ref base T (16/18); new base C (21/16); total (37/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.80e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01 | |||||||||||
TTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCAACCCCGTAATTTTGGTGAG > NC_000913/1433621‑1433808 | ttCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCTAACATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCa < 1:2380591/100‑1 (MQ=255)ttCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCa < 2:1595001/100‑1 (MQ=18) tCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCa < 2:563742/99‑1 (MQ=18) ggtggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGaa > 1:2409655/1‑100 (MQ=33) ggtggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGaa < 1:1301520/100‑1 (MQ=32) tggtCTGATCCAGCGTTGATATCCGACATTACGCAAAAATTCGAGACATAGCGAGACTCTATACGAAATACTTGTTCATATGATAATGCATCCCAGaaaa > 1:3051706/1‑100 (MQ=11) ggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAaat < 1:3119245/100‑1 (MQ=32) ggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAaat < 2:2713209/100‑1 (MQ=33) ggtCAGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaat < 1:2668533/100‑1 (MQ=12) gtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAaata > 1:999937/1‑100 (MQ=33) gtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAaata > 2:2549719/1‑100 (MQ=33) gtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAaata > 1:2507597/1‑100 (MQ=33) tCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAaataa > 2:955456/1‑100 (MQ=33) cTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaataat > 2:2811230/1‑100 (MQ=17) cTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaataat > 2:1163452/1‑100 (MQ=18) cTGATCCAGCGTTGACATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAaataat < 2:1490559/100‑1 (MQ=18) tGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGATAAATAGTGAGACTATATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAAtt < 1:1397982/100‑1 (MQ=255) gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAAttn < 2:187610/100‑2 (MQ=18) gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa > 2:831964/1‑100 (MQ=18) gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa > 2:3135796/1‑100 (MQ=18) gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa > 1:1135063/1‑100 (MQ=18) aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTAt > 2:895999/1‑100 (MQ=33) aTCCAGCATTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTTTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAaataat < 1:3004335/97‑1 (MQ=12) tCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTAt > 1:1978602/1‑99 (MQ=17) tCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATg > 2:461569/1‑100 (MQ=18) ccAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGa > 2:1316721/1‑100 (MQ=18) naGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGat > 1:1063514/2‑100 (MQ=18) aGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGatt > 2:1889411/1‑100 (MQ=33) gtgatGATATACGACAGTAAGCAAAAATGCTATACATAGTGAGACTTTATACGAAATTGGTGTACATATCAAAATGCAGCCCAGAAAATAATTATGatt < 1:209138/95‑1 (MQ=11) gTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGattatt > 1:19140/1‑100 (MQ=33) ttGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAttt > 1:2760843/1‑100 (MQ=18) tGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATGTCATAATGCATCTCAGAAAATAATTATGATTAtttt < 2:2387556/100‑1 (MQ=18) tGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAtttt < 2:3323120/100‑1 (MQ=18) tGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTAATATCATAATGCATCTCAGAAAATAATTATGATTAtttt > 1:391586/1‑100 (MQ=17) natatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTAttttt > 1:2002999/2‑100 (MQ=18) gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTAGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTAttttt < 1:2518657/100‑1 (MQ=17) atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtg < 1:3259272/100‑1 (MQ=32) tatCCGACAGTAAGCAAAAATTCGATACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtgt < 1:3321852/100‑1 (MQ=17) tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtgt > 1:349356/1‑100 (MQ=33) tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgt < 2:367322/100‑1 (MQ=18) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtgtg > 1:2568285/1‑100 (MQ=33) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg > 1:733331/1‑100 (MQ=18) tCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGc > 2:2742565/1‑100 (MQ=32) tCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGc > 1:754996/1‑100 (MQ=18) tCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGc > 1:1092221/1‑100 (MQ=18) ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATCATAATTAttttctgaga < 2:206447/100‑7 (MQ=11) cGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCat > 1:1315060/1‑100 (MQ=17) gaaGCAAAAATGCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCATAAAATAATTATGATTATTTTTGTGCATATTtg < 2:989422/98‑1 (MQ=255) naGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta > 2:2544910/2‑100 (MQ=18) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 2:3308283/100‑1 (MQ=18) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 1:331519/100‑1 (MQ=17) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 2:2977338/100‑1 (MQ=18) aaGAAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 1:1070396/100‑1 (MQ=16) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 2:3133371/1‑100 (MQ=33) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat < 2:397878/100‑1 (MQ=33) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 2:1864517/1‑100 (MQ=17) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 1:2182845/1‑100 (MQ=18) gCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatg < 2:266991/100‑1 (MQ=18) gCAAAAATTCAAGACATATTGAGACTCTATATGAAATTGGTGTTCATATCATAATGAATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatg < 1:451354/100‑1 (MQ=11) aaaaaTTCGAGACATAGTGAGACTCTATAAGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatgta < 2:2368868/100‑1 (MQ=12) aaaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTac > 1:147903/1‑100 (MQ=33) aaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTaca > 1:2771871/1‑100 (MQ=18) aTTCGATACATAGTGAGACTCTATACGAAATTGGTGTTAATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACg < 2:108451/100‑1 (MQ=255) aTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACg < 1:228723/100‑1 (MQ=32) tCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGtt < 2:3226010/100‑1 (MQ=31) tCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGt > 1:2977191/1‑99 (MQ=14) cGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTc > 2:93837/1‑100 (MQ=11) cGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTc > 2:69532/1‑100 (MQ=11) gagtccctATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATTATTATTTTTGTTTATATTTTTATGTACACGTTCGAAAATAAACg < 1:415077/94‑1 (MQ=1) gtattCATATCATAATGCATCTCAGAAAATAATTATTATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCAAccc < 1:2635786/94‑1 (MQ=0) tgtaTCTCGGAAATTAATTGTTTTTTTTTTTGTGCTTATTTTTTTGTACACGTTTGTAAATAATCGAATGTTTATGCAACGCCGTAATTTTGGTgag < 1:950329/94‑1 (MQ=0) | TTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCAACCCCGTAATTTTGGTGAG > NC_000913/1433621‑1433808 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |