| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,633,966 | 0 | T | C | 53.8% | ‑5.7 / 10.7 | 13 | pseudogene (426/528 nt) | nohQ | pseudogene, Qin prophage; predicted packaging protein;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin |
| Reads supporting (aligned to +/- strand): ref base T (0/6); new base C (7/0); total (7/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAA > NC_000913/1633876‑1634055 | cacccccGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAAc < 1:3797682/98‑1 (MQ=1) gggaTCGGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTCTTCCAGTTCCGGAAAACGCCGCTg < 2:3493986/98‑1 (MQ=0) gggaTCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGATTCAGGAAATCAACATGTCGGGTTTCCAGTTCCGGGAAACGCCGCtc > 1:3368044/4‑100 (MQ=1) cacccAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAgg < 2:4442134/98‑1 (MQ=1) attgCTTTGATTATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGc < 2:6690353/98‑1 (MQ=0) tattatatCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCAc < 1:8494260/97‑1 (MQ=1) tgctatCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCt < 2:5477537/98‑1 (MQ=1) tcacGTTTCATGAAATCAACATGTCGGTATTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGAt < 2:3584923/98‑1 (MQ=0) cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCcgcg > 1:12761905/5‑101 (MQ=0) cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCcgcg > 1:5658642/5‑101 (MQ=0) cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCcgcg > 1:4952798/5‑101 (MQ=0) cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCcgcg > 1:2986271/5‑101 (MQ=0) cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGCCAATCcgcg > 1:12842701/5‑101 (MQ=0) ttaaGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACACAGCTGCACCCACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGa < 2:4507280/98‑1 (MQ=1) aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGaacg > 1:12358707/4‑99 (MQ=0) aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTATTCGGCCATCCGCGACAGCACGAAAg > 1:13918625/4‑101 (MQ=0) aataTGTCGGTTTGCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGaacg < 1:5003389/98‑3 (MQ=0) tctgTTTTCCAGTTCCGTAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGa < 2:6894158/98‑1 (MQ=0) cgattttCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGaa < 1:6045292/98‑1 (MQ=1) | CAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAA > NC_000913/1633876‑1634055 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |