| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_hsa_tpiA | 1,642,370 | 0 | C | G | 20.0% | 36.2 / 11.4 | 25 | intergenic (+303/+119) | cspF/quuQ | Qin prophage; cold shock protein/Qin prophage; putative antitermination protein Q |
| Reads supporting (aligned to +/- strand): ref base C (12/8); new base G (2/3); total (14/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.23e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGC > NC_000913_3_hsa_tpiA/1642235‑1642515 | gACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt < 1:139817/149‑1 (MQ=255) cATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTCCCAGAAAATGTAACAGACATGAGTCCCGGGATGCGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt < 1:421140/149‑3 (MQ=255) ttGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGtttt < 2:404849/145‑1 (MQ=255) ttGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGtttt > 1:404849/1‑145 (MQ=255) tACTGGAAGAAAACATTTTAAACATCAGGGAAATAACCAAGTGACCAGCTAAATAATAAGTTAACAGACATGAGCCCCGGTATGAGATTAACCATTACCATTGCCCCCTCTAAAGCACAAAACCCTCTCATCAGCGGGTGTTCTACttt < 2:103054/149‑1 (MQ=255) tGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTCTCTACTTTTTc > 1:2351/1‑149 (MQ=255) ggAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAAGAATAATTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt < 1:228082/149‑1 (MQ=255) aaaaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACg > 1:263821/1‑149 (MQ=255) aaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTc > 1:428573/1‑149 (MQ=255) cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg < 1:5988/149‑1 (MQ=255) cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg < 1:97440/149‑1 (MQ=255) aGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGccc > 1:329893/1‑149 (MQ=255) aaCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa < 1:320583/149‑1 (MQ=255) taggtggccAGCTAAATAATAAGTAAACAGACATGAGTCCCGGGATGAGATTCAACATTACAATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCCGATTTTCTACTTTTTCTTAACGTCGAGTATACAAAGCCCATCGTTGaaaa > 2:269022/8‑149 (MQ=255) aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtctt > 1:213523/1‑149 (MQ=255) cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt < 2:270342/149‑1 (MQ=255) gAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcat > 1:356127/1‑149 (MQ=255) gATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCACACATCATGTCGcc > 2:357907/1‑149 (MQ=255) tGAGATTCATCATTACCATTGCCCCATTTAAAGCATAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATATAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAGATGCGAACCTCTTGGCGACAt < 2:203193/149‑1 (MQ=255) tGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCAt > 2:429349/1‑149 (MQ=255) ttACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAAATCATGTCGCCATCTTCAGCaaaaa > 1:72604/1‑149 (MQ=255) cATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATGAt > 2:229156/1‑149 (MQ=255) aaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAt < 1:104866/71‑1 (MQ=38) aaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAt > 2:104866/1‑71 (MQ=38) gCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTgcgc > 1:203989/1‑149 (MQ=255) | GACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGC > NC_000913_3_hsa_tpiA/1642235‑1642515 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |