| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,634,949 | 0 | G | C | 29.2% | 52.9 / 18.2 | 24 | A300G (GCT→GGT) | stfQ | Qin prophage; putative side tail fibre assembly protein |
| Reads supporting (aligned to +/- strand): ref base G (9/8); new base C (0/7); total (9/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.23e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTCATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAG > NC_000913/1634869‑1635048 | cTCATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaac < 2:1119043/100‑4 (MQ=1)cTCATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaac < 2:274128/100‑4 (MQ=1) tCATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacg < 2:2121181/100‑5 (MQ=1) tCATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacg < 2:1158610/100‑5 (MQ=1) tCATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacg < 1:1599510/100‑5 (MQ=1) tCATTCTGAATGCCATTATGCAAGCCTCACAATATAGGTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacn < 2:2638520/100‑5 (MQ=1) cATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgc < 1:1838878/100‑6 (MQ=1) aCGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTAtgtg < 1:2188381/100‑1 (MQ=32) aCGGTGTTTACCGCGTTACCAGCAGCGTTAACGGTGATGGTTTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTAtgtg < 1:132034/100‑1 (MQ=12) cGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGc > 2:35515/1‑100 (MQ=33) cGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGc > 1:2786825/1‑100 (MQ=33) cGGTGTTTTCCGCATTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGc < 1:1575927/100‑1 (MQ=18) gtgtTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCAt > 2:2493148/1‑100 (MQ=35) gtTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATgc > 1:1154687/1‑100 (MQ=38) gtTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATgc > 1:1154733/1‑100 (MQ=38) tCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATgcgcgc > 1:135238/1‑100 (MQ=255) tCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATgcgcgc > 2:2203773/1‑100 (MQ=255) ccGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGcc < 2:663569/100‑1 (MQ=255) ccGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGcc < 2:2545792/100‑1 (MQ=255) ccGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGcc < 2:1486701/100‑1 (MQ=255) cgcgTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCt > 2:2003630/1‑100 (MQ=255) cgcgTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCt < 2:2914527/100‑1 (MQ=255) cgTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTgc < 1:2006560/100‑1 (MQ=255) ccAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTg > 2:1763180/1‑100 (MQ=255) cagcagCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgc < 1:906372/100‑1 (MQ=255) cagcagCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgc < 1:2147732/100‑1 (MQ=255) cagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcagc < 2:2664749/100‑1 (MQ=255) agcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcagca > 2:2427401/1‑100 (MQ=255) agcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcagca > 2:552864/1‑100 (MQ=255) gcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGcagcag < 1:2343425/100‑1 (MQ=255) gcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGcagcag > 2:498246/1‑100 (MQ=255) gcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCAACAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGcagcag < 2:2641048/100‑1 (MQ=255) | CTCATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAG > NC_000913/1634869‑1635048 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |