Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 580,419 | 0 | C | G | 73.2% | 32.4 / 43.3 | 56 | I121M (ATC→ATG) | nohB | bacteriophage DNA packaging protein |
Reads supporting (aligned to +/- strand): ref base C (7/8); new base G (27/14); total (34/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC > NC_000913/580327‑580491 | cgcgcgcAGGCCGCTGCAGAGGAGCTGAAACCTGCCAGAGACTCCGCTGAAGCGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATc < 2:1360653/90‑1 (MQ=11) cgTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATc > 2:4133085/1‑93 (MQ=33) tGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGc > 2:1569019/1‑93 (MQ=33) tGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGc < 2:1233102/93‑1 (MQ=33) agcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGc < 1:4686719/92‑1 (MQ=2) agcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGc > 1:4627147/2‑93 (MQ=2) agcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGc < 1:319320/92‑1 (MQ=2) gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa > 2:3354427/1‑93 (MQ=34) gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa < 1:4299638/93‑1 (MQ=33) gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa < 1:4134047/93‑1 (MQ=33) gcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGcc < 1:1577885/93‑2 (MQ=11) gcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGcc > 2:4269305/1‑92 (MQ=2) gcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGcc < 2:3288658/93‑2 (MQ=2) acgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGcc < 2:1282009/92‑2 (MQ=2) cgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAg < 2:479552/93‑1 (MQ=33) cgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGccg > 1:4518645/1‑91 (MQ=2) cgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGccg < 2:2556189/93‑3 (MQ=2) gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg > 2:340795/1‑93 (MQ=2) cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt < 1:4192499/93‑1 (MQ=33) cAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgga > 1:985303/1‑92 (MQ=2) cAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgga > 1:2512023/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:2716779/1‑91 (MQ=2) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag > 1:4528104/1‑91 (MQ=2) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag > 1:2716778/1‑91 (MQ=2) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa > 2:2833563/1‑93 (MQ=33) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa > 2:2833688/1‑93 (MQ=33) ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGTGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGCCGCGGATGGCCGgaga > 2:197374/1‑93 (MQ=2) ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGgaga > 1:265162/1‑93 (MQ=2) ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGgaga > 1:3004157/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:2596499/92‑3 (MQ=2) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt > 1:1980229/2‑91 (MQ=2) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt > 1:3968793/2‑91 (MQ=2) gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta > 2:1251750/1‑90 (MQ=2) gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta > 2:1211934/1‑90 (MQ=2) aTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaaggag > 1:1608377/1‑89 (MQ=2) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg < 2:2106863/93‑1 (MQ=32) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 2:449384/1‑93 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc < 2:4433060/93‑1 (MQ=31) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct < 2:534810/89‑1 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 2:3746520/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 1:3807830/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 2:799262/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 1:399768/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 2:1651226/5‑93 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 1:4466446/4‑92 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt < 1:4528589/90‑2 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg < 1:365722/90‑1 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 1:2669311/4‑93 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa < 1:2839249/89‑1 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa < 2:2390395/89‑1 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc > 1:4014100/4‑93 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc > 1:1502421/4‑93 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc < 1:822790/90‑1 (MQ=0) | CGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC > NC_000913/580327‑580491 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |