Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 580,431 | T→A | I125I (ATT→ATA) | nohB → | bacteriophage DNA packaging protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 580,431 | 0 | T | A | 75.7% | 24.3 / 25.3 | 37 | I125I (ATT→ATA) | nohB | bacteriophage DNA packaging protein |
Reads supporting (aligned to +/- strand): ref base T (7/2); new base A (15/13); total (22/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.62e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01 |
GATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC > NC_000913/580339‑580491 | gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt > 1:1056763/1‑93 (MQ=33) gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt > 2:9768515/1‑93 (MQ=33) gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt > 2:9361483/1‑93 (MQ=33) gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt < 2:5515691/93‑1 (MQ=33) gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt < 2:487398/93‑1 (MQ=33) gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt < 2:2412398/93‑1 (MQ=33) gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt < 1:7974575/93‑1 (MQ=33) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg > 1:1783885/1‑93 (MQ=32) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg > 2:6911456/1‑93 (MQ=32) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg > 1:7614774/1‑93 (MQ=32) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg > 1:6405752/1‑93 (MQ=32) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg > 1:2346704/1‑93 (MQ=32) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg > 1:4950852/1‑93 (MQ=32) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 1:5557290/1‑93 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 2:8191741/1‑93 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc < 1:8822324/93‑1 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 1:9598103/1‑93 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc < 2:874959/93‑1 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 2:8191742/1‑93 (MQ=31) actaaaaaaTGCAAAAGATTCCGCTGCGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 1:9057116/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGATAAGTAGCCAGCATtct < 2:2686827/89‑1 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 1:3841349/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct < 1:8135479/89‑1 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 1:9118193/5‑93 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt < 1:6060112/90‑2 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt < 2:7012888/90‑2 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 2:6528071/4‑92 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 1:5028649/4‑92 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 2:1695411/4‑92 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 2:1191409/4‑92 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 2:1187430/4‑92 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGCGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt < 2:7012890/90‑2 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg < 2:3433566/90‑1 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg < 2:3756741/90‑1 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 2:4904781/4‑93 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 2:9321164/4‑93 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 1:7450646/4‑93 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg < 1:7348744/90‑1 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg < 2:7719151/90‑1 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa > 1:6398876/5‑93 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa < 1:2703825/89‑1 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa > 1:5198769/5‑93 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa > 1:5517753/5‑93 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCACCGGCGCTTCACGGAGCTGGaaaa > 2:8383555/5‑93 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc < 2:4028738/90‑1 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc < 2:860820/90‑1 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc < 1:9697001/90‑1 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc > 1:1886367/4‑93 (MQ=0) | GATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC > NC_000913/580339‑580491 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |