| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_bme_pgi | 1,377,961 | 0 | G | . | 21.1% | 46.0 / 10.3 | 19 | coding (78/2268 nt) | ycjT | putative family 65 glycosyl hydrolase |
| Reads supporting (aligned to +/- strand): ref base G (7/8); new base . (2/2); total (9/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCAGCAACTGGTGGCGGCGTTGTATCAATCCGCAGAAACAGGGACGCGTGTGGAATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATCGGGCGGGAAAAGGTGAAATCAACGAACTGGTGAACCTGCCTGATGTCGTGGGGATGGAGA > NC_000913_3_bme_pgi/1377827‑1378109 | ccAGCAACTGGTGGCGGCGTTGTATCAATCCGCAGAAACAGGGACGCGTGTGGAATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCtt > 2:94161/1‑150 (MQ=255) ggtggtGGCGTTGTATCAATCCGCAGAAACAGGGACGCGTGTGGAATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCA‑GGGAACGGCTATCTTGGGCTTcgcg > 4:21018/1‑150 (MQ=255) ggtggCGGCGTTGTATCAATCCGCAGAAACAGGGACGCGTGTGGAATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCA‑GGGAACGGCTATCTTGGGCTTcgcg > 4:21016/1‑150 (MQ=255) cgTTGTATCAATCCGCGGAAACAGGGACGCGTGTGGGATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCgccagcca > 2:86279/1‑150 (MQ=255) gTATCAATCCGCAGAAACAGGGACGCGTGTGGAATTATGACCAGGCCAGTAACTTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGGTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATgaa < 2:104824/150‑1 (MQ=255) tCAATCCGCAGAAACAGGGACGCGTGTGGAATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATgaagaa < 1:94161/150‑1 (MQ=255) gtgGAATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGa < 1:54217/150‑1 (MQ=255) gAATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGt > 3:18314/1‑150 (MQ=255) gacggggCCAGTAAGGTTATCAGAACCCCGTTTCAGCAAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGACCGGCTATCTCGGGCTTCGCGCCAGCCCTGAAGAAGATTACACGCGCCAGACGCTAGGGGTGTATCTGGc < 4:86424/145‑1 (MQ=255) aGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATcgggcggg < 4:18314/150‑1 (MQ=255) tGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGa < 2:111305/89‑1 (MQ=255) tGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGa > 1:111305/1‑89 (MQ=255) tGAACAAGTATGCATCGCTGATGGCGCA‑GGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATCGGGCGGGAAAAGGTGAAATCAACGAACTGGTGaa < 3:21016/150‑1 (MQ=255) tGAACAAGTATGCATCGCTGATGGCGCA‑GGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATCGGGCGGGAAAAGGTGAAATCAACGAACTGGTGaa < 3:21018/150‑1 (MQ=255) aTCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATCGGGCGGGAAAAGGTGaa > 2:47415/1‑121 (MQ=255) aTCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATCGGGCGGGAAAAGGTGaa < 1:47415/121‑1 (MQ=255) cGCTGATGGCGCAGGGGAACGGCTCTCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATCGGGCGGGAAAAGGTGAAATCAACGAACTGGTGAACCTGCCTGATGTCg < 1:86279/150‑1 (MQ=255) cTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCCGACGCGAGGGACGTATCTGGCGGGGCCGTAACACCTGGCGGGGAAAGGAGAAATCAACGAACTGCGTAACCTGCCTGaggaggtc > 4:22245/1‑143 (MQ=255) aGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATCGGGCGGGAAAAGGTGAAATCAACGAACTGGTGAACCTGCCTGATGTCGTGGGGATGgaga > 4:28059/1‑150 (MQ=255) | CCAGCAACTGGTGGCGGCGTTGTATCAATCCGCAGAAACAGGGACGCGTGTGGAATTATGACCAGGCCAGTAACGTTATCAGAACCCCATTTCAGCCAGCATACCCTGAACAAGTATGCATCGCTGATGGCGCAGGGGAACGGCTATCTTGGGCTTCGCGCCAGCCATGAAGAAGATTACACGCGCCAGACGCGAGGGATGTATCTGGCGGGGCTGTATCATCGGGCGGGAAAAGGTGAAATCAACGAACTGGTGAACCTGCCTGATGTCGTGGGGATGGAGA > NC_000913_3_bme_pgi/1377827‑1378109 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |