Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,135,975 | 0 | G | C | 41.2% | 7.9 / 10.9 | 18 | A138P (GCC→CCC) | flgH | flagellar protein of basal‑body outer‑membrane L ring |
Reads supporting (aligned to +/- strand): ref base G (5/5); new base C (0/7); total (5/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.41e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACG > NC_000913/1135922‑1136038 | cTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAAt > 2:120941/1‑68 (MQ=255) cTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAAt > 2:308445/1‑68 (MQ=255) ccGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATa < 2:64145/67‑1 (MQ=255) cGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATAc < 1:100063/68‑1 (MQ=255) aGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGt > 1:354317/1‑68 (MQ=255) cGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTtgacgg < 2:687561/68‑1 (MQ=255) cGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTATCGGCACGTtgacgg < 2:98424/68‑1 (MQ=255) cGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTtgacgg < 2:367625/68‑1 (MQ=255) ggtggtAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTtgacggt < 1:230583/68‑1 (MQ=255) cggtCAACGGAAAGGGCGGGTCCAATGCCAGCAATACCTTTAGCGGCACGTTGACTGTGACGGTTGAc < 2:50310/65‑1 (MQ=255) cGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc < 1:430572/68‑1 (MQ=255) tCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAg < 2:289358/68‑1 (MQ=255) cAACGGAAAGGGCGGGGCCATGGCCAGCAATACCTTTAGCGGCACGTTGACGTTGACGGTTGACCAgg < 2:616246/68‑1 (MQ=255) aaaGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTgg < 2:267568/68‑1 (MQ=255) ggCGGGGCAAATGCCAGCAATACCTTTAGCGGCAAGTTGACGGTGACGGTTGACAAGGTACTGGTaaa > 2:174233/1‑68 (MQ=255) gCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGCCCAGGTACTGGTCAAc > 1:181129/1‑68 (MQ=255) gCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAAc < 1:321871/68‑1 (MQ=255) cGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACg < 1:621422/68‑1 (MQ=255) | CTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACG > NC_000913/1135922‑1136038 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |