| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 580,422 | 0 | A | C | 54.5% | 73.9 / 76.5 | 55 | A122A (GCA→GCC) | nohB | bacteriophage DNA packaging protein |
| Reads supporting (aligned to +/- strand): ref base A (14/11); new base C (20/10); total (34/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.78e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC > NC_000913/580330‑580491 | gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa > 2:2797544/1‑93 (MQ=33)gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa < 1:5953495/93‑1 (MQ=33) cgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCTGATCGCAg < 1:2932367/93‑1 (MQ=18) cgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAg < 1:914569/93‑1 (MQ=33) cgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAg < 1:2113296/93‑1 (MQ=33) gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg > 1:6365109/1‑93 (MQ=2) gcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAgg < 2:480141/93‑1 (MQ=33) gcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAgg > 2:6513731/1‑93 (MQ=33) gcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAgg < 1:2711374/93‑1 (MQ=33) gcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAgg > 2:3657582/1‑93 (MQ=33) cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt < 2:4046005/93‑1 (MQ=33) aGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTg > 1:2190377/1‑93 (MQ=34) aGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTg > 1:5727147/1‑93 (MQ=34) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag > 2:730446/1‑91 (MQ=2) aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag > 1:3945318/1‑91 (MQ=2) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCCGGTGa > 1:5627012/1‑93 (MQ=21) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa > 1:914279/1‑93 (MQ=34) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa > 2:429108/1‑93 (MQ=34) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa < 2:4841879/93‑1 (MQ=33) ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa > 1:2569000/1‑93 (MQ=34) gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa > 2:1895694/4‑93 (MQ=2) gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa > 1:4646876/4‑93 (MQ=2) gCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGaa > 1:1765562/1‑93 (MQ=34) ccGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGaaa > 2:4512280/1‑93 (MQ=33) cagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaag < 2:6298634/91‑2 (MQ=2) cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt > 2:2446193/1‑93 (MQ=33) cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt < 1:3445000/93‑1 (MQ=33) cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt < 2:6290417/93‑1 (MQ=33) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt > 2:3088287/2‑91 (MQ=2) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt > 2:2315201/2‑91 (MQ=2) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt > 1:2393910/2‑91 (MQ=2) agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaaga > 2:1361427/2‑91 (MQ=2) gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta > 1:197416/1‑90 (MQ=2) gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt > 1:561876/1‑93 (MQ=33) gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt > 1:4521327/1‑93 (MQ=33) aTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagtag > 1:6424416/1‑89 (MQ=2) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg > 2:4699916/1‑93 (MQ=32) aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg < 2:82582/93‑1 (MQ=32) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 2:5451048/1‑93 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 2:1114802/1‑93 (MQ=31) tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc > 1:4146739/1‑93 (MQ=31) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct > 1:1164190/5‑93 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct < 1:6169472/89‑1 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct < 2:4832437/89‑1 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct < 2:3922553/89‑1 (MQ=0) actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGGAGCCAGCATtct > 1:3644108/5‑93 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 2:1537497/4‑92 (MQ=0) ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt > 1:3982570/4‑92 (MQ=0) taaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtcttg < 1:2090323/91‑3 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 1:4012051/4‑93 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 1:3810454/4‑93 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 2:617944/4‑93 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg < 1:2354154/90‑1 (MQ=0) gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg > 1:1279293/4‑93 (MQ=0) gtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcgg < 1:2315521/93‑1 (MQ=0) gtggtgGAAACCGCACTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcgg < 2:2569660/93‑1 (MQ=0) gtgcTTTCGTGCTGTCGAGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCGCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGaaa < 1:6028783/90‑1 (MQ=0) cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa < 2:5164712/89‑1 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc > 1:4710213/4‑93 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc > 2:1207626/4‑93 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc < 2:6296423/90‑1 (MQ=0) acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc > 2:2131119/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: ‑ |