| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 923,656 | 0 | C | A | 55.6% | 16.1 / 25.8 | 20 | D514Y (GAT→TAT) | PP_0804 | protein secretion ABC transporter permease/ATP‑binding protein |
| Reads supporting (aligned to +/- strand): ref base C (0/8); new base A (10/0); total (12/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.29e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGCGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCGGCGACT > NC_002947/923530‑923787 | tGGCGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCGCGGCGGGGGATACAGCCGGcc > 2:312979/1‑136 (MQ=255)tGGCGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCGCCGCGGGTTATACCGCCGGcg > 1:395988/1‑135 (MQ=255) cGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAAt < 2:287112/136‑1 (MQ=255) gATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATACCGCCGGCCAATg > 2:331501/1‑136 (MQ=255) aGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCGGGGTATACAGCCGGCCCATGTTGAGCACTg > 2:339673/1‑136 (MQ=255) tcGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCGCCGCGGGGTATGCAGCCCGCCAAAGTTGAGCACAGGGGGGTTGt > 1:121486/1‑136 (MQ=255) tcGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGGGGGGTATGCAGCCCGCCAAACTTCAGCACAGGCGGGGAGt > 1:349574/1‑136 (MQ=255) cccGGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCtt < 2:210135/136‑1 (MQ=255) tGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCCGCGGGGTATAAACCCGGCCCATGTTAAGCACAGGCGGGTAGTCTTAGCTGTATCGGAAGtt > 2:381143/1‑136 (MQ=255) gCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCCGCAGGATATAAACCCCCCCAAAGTTAAACACAGGCCGGGAGTCTTAGCAGTATCGGAAAtttt > 1:59869/1‑134 (MQ=255) gCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCCCCGGGGTATACAGCCCGCCCAATATGAGCACAGGCGGGGAGTAATACGTGTTGAGGAAGTTgg > 2:374674/1‑136 (MQ=255) ccccACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCCGCGGGATATACACCCCGCCAAAGTTAAGCACAGGCGGGGAGTATTCGCGGTAGCGGGAGtttgg > 2:102270/1‑133 (MQ=255) ccTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCg < 1:310721/136‑1 (MQ=255) ccTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCg < 1:312979/136‑1 (MQ=255) agcaAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCGCGGCGGGGATTACAGCCCGCCCATGTTCAGCACTGGCGGGGATTAATAGCTGTAGCGGGAGTTTGACAGGCGCCGGCGGGAca > 2:286422/1‑135 (MQ=255) gcTCGCCCCGGGGGGTATACAGCCGGCCCATGTTGAGCACTGGGGGGGTGTGTTGGGTGTAGCGGGAGTTGGGCTGGCGCAGGCGGTTCTCGCCGTGAATGGCCGGCAGGGGGACGCGGGCGTCGCCGTCGGGGTg > 2:372789/1‑136 (MQ=255) ggCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCg < 2:138321/136‑1 (MQ=255) ggCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCg < 2:135600/136‑1 (MQ=255) ggtGGATATCCAGCCGGCCAATGTTCAGCACTGGCTGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCgg < 1:58374/134‑1 (MQ=255) atatCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCGGCGACt < 1:381143/136‑1 (MQ=255) | TGGCGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCGGCGACT > NC_002947/923530‑923787 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |