Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 580,422 | A→C | A122A (GCA→GCC) | nohB → | bacteriophage DNA packaging protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 580,422 | 0 | A | C | 76.2% | 17.9 / 21.2 | 42 | A122A (GCA→GCC) | nohB | bacteriophage DNA packaging protein |
Reads supporting (aligned to +/- strand): ref base A (5/5); new base C (24/8); total (29/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01 |
GCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC > NC_000913/580330‑580491 | gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa < 1:4299638/93‑1 (MQ=33) gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa < 1:4134047/93‑1 (MQ=33) gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa > 2:3354427/1‑93 (MQ=34) cgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAg < 2:479552/93‑1 (MQ=33) gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg > 2:340795/1‑93 (MQ=2) cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt < 1:4192499/93‑1 (MQ=33) cAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgga > 1:2512023/1‑92 (MQ=2) cAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgga > 1:985303/1‑92 (MQ=2) aGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTg < 1:2971500/93‑1 (MQ=33) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAGCCGCATTCTGCCCGTTCGTGCTGTCGCGGATGGCCggag < 2:2467816/93‑3 (MQ=2) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACGGCATTCTGCACGTTCGTGCTGTCGCGGATGGCgggag > 1:3736710/1‑91 (MQ=2) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag > 1:4528104/1‑91 (MQ=2) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag > 1:2716779/1‑91 (MQ=2) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag > 1:2716778/1‑91 (MQ=2) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa > 2:2833688/1‑93 (MQ=33) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa > 2:2833563/1‑93 (MQ=33) ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGTGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGCCGCGGATGGCCGgaga > 2:197374/1‑93 (MQ=2) ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGgaga > 1:3004157/1‑93 (MQ=2) ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGgaga > 1:265162/1‑93 (MQ=2) ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTATGCACGTTCGTGCTGTCGCGGATGGCCGGCGa > 2:3593784/1‑93 (MQ=2) gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa > 1:3149965/4‑93 (MQ=2) gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa > 1:573751/4‑93 (MQ=2) gCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGaa > 1:3855357/1‑93 (MQ=33) cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt > 1:2481117/1‑93 (MQ=33) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt > 1:563521/2‑91 (MQ=2) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt > 1:3968793/2‑91 (MQ=2) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt < 1:2596499/92‑3 (MQ=2) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt > 1:1980229/2‑91 (MQ=2) gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta > 2:1211934/1‑90 (MQ=2) gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta > 2:1251750/1‑90 (MQ=2) aTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaaggag > 1:1608377/1‑89 (MQ=2) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg < 2:2106863/93‑1 (MQ=32) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc < 2:4433060/93‑1 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 2:449384/1‑93 (MQ=31) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 2:3746520/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 2:799262/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct < 2:534810/89‑1 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 2:1651226/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 1:399768/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 1:3807830/5‑93 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt < 1:4528589/90‑2 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 1:4466446/4‑92 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg < 1:365722/90‑1 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 1:2669311/4‑93 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa < 2:2390395/89‑1 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa < 1:2839249/89‑1 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc < 1:822790/90‑1 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc > 1:4014100/4‑93 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc > 1:1502421/4‑93 (MQ=0) | GCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC > NC_000913/580330‑580491 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |