| 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 | 97.0% | 79.5 / ‑4.3 | 44 | intergenic (‑44/+273) | pinR/ynaE | Rac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage |
| Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (23/19); total (23/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.22e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTA > NC_000913/1433621‑1433785 | ttCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCa > 1:1470704/1‑100 (MQ=18) tCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAg > 2:2483108/1‑100 (MQ=18) cGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGa < 1:530735/100‑1 (MQ=17) gtggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGaaa > 1:897915/1‑100 (MQ=18) ggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGCCTCTATACGAAATTGGTGTTCATATCATAATGCATCCCCGaaaa > 2:1508979/1‑99 (MQ=255) gtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACCTAGTGAGACTCTATACGAAATTGGTGTTCATAGCATAATGCACCCCAGAAaat > 2:919923/1‑99 (MQ=11) tCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaat < 1:472844/98‑1 (MQ=17) gtactCCGGCGTTGATATCTTACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaataat > 2:3000499/5‑100 (MQ=11) cTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaataat > 2:37172/1‑100 (MQ=18) cTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCACAATGCATCCCAGAAaataat < 2:1236560/100‑1 (MQ=12) gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa > 1:2351583/1‑100 (MQ=18) gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa > 1:1974082/1‑100 (MQ=18) gATCCAGCGCTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa > 1:3030457/1‑100 (MQ=255) aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTAt > 1:2894788/1‑100 (MQ=18) aGCGTTGATATCCGATAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGat < 2:1768050/99‑1 (MQ=12) aGCGTTGATATCCGACAGTAAGCAAAAATTCGAGATATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGat < 1:1198019/99‑1 (MQ=12) gCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATACATCCCAGAAAATAATTATGatta > 1:370808/1‑100 (MQ=16) cGTTGATATCCGACAGTGAGCAAAAATTCGAGACATAGCGAGACTCTAGACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATg < 2:2544321/95‑1 (MQ=11) gTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGattatt > 1:2766388/1‑100 (MQ=18) tGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCAAGAAAATAATTATGATTAtttt > 1:358274/1‑100 (MQ=16) atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg > 1:290471/1‑100 (MQ=18) atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg < 1:2391079/100‑1 (MQ=17) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGCCTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg < 2:2239088/100‑1 (MQ=12) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg < 2:332670/100‑1 (MQ=17) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg < 1:2137354/100‑1 (MQ=17) atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATCATAATTAttttctgg < 1:1883894/100‑5 (MQ=11) tCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGc > 2:1533310/1‑100 (MQ=18) tCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGc < 2:856682/100‑1 (MQ=17) ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATACATCCCAGAAAATAATTATGATTATTTTTGTGCa > 1:163821/1‑100 (MQ=16) gACAGTAATCAGAAAGTCGAGTCATGGTGGGGCTCTAGACGAAATTGGTGTTCATATCATAATGCATCCCAGATTATAATTATGATTATTTTTGTGCatg < 2:2123849/100‑2 (MQ=11) gACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCata > 2:2322170/1‑100 (MQ=18) aCAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCatat > 2:2909558/1‑100 (MQ=18) cAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATAtt < 1:1364884/100‑1 (MQ=17) gTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGTTGTTCACTTCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCatat < 2:1678050/97‑1 (MQ=255) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 2:1254612/1‑100 (MQ=18) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 2:1068893/1‑100 (MQ=18) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 1:14762/1‑100 (MQ=18) cAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatgt > 1:2022416/1‑100 (MQ=18) aaaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTac < 2:959230/100‑1 (MQ=18) aaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCTGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTaca < 2:1935206/100‑1 (MQ=12) aaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTaca < 1:1533328/100‑1 (MQ=18) aTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACg < 2:2351556/100‑1 (MQ=18) aatccataCGTAATTGGTGTTCATATCTTAATGCATACCATAATATAATTATGATTATTTTTGTTCATATTTGTATGTACACGTTCGAAAATAAAc < 2:1557914/91‑1 (MQ=1) tatgAAATTGGTGTTCATATCATAATGCGTCTCAGAAAATAATTATTATTATTTTTGTTCATATTTTTATGTTCACTTTCGAAAATAAACGAATGCgga < 2:709528/96‑3 (MQ=1) | TTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTA > NC_000913/1433621‑1433785 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |