Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 591,178 | C→T | 42.4% | A465V (GCA→GTA) | rpoC → | DNA‑directed RNA polymerase subunit beta' |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 591,178 | 0 | C | T | 42.4% | 11.6 / 25.0 | 28 | A465V (GCA→GTA) | rpoC | DNA‑directed RNA polymerase subunit beta' |
Reads supporting (aligned to +/- strand): ref base C (8/8); new base T (5/7); total (13/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.18e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 |
TAGACTTGGTATTCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACTTTAGAAAGAAAAGAT > USA300TCH1516_ALE/591041‑591308 | tAGACTCGGTATTCAAGCATTTACACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCCCTTGTAACAACATCCTACAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGcaca < 2:264275/141‑1 (MQ=255) gTATTCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAg < 2:338593/141‑1 (MQ=255) ttCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCaa < 1:101106/141‑1 (MQ=255) tCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCACGAGAGGTACAACCTGACGCCAg > 2:169421/1‑141 (MQ=255) gCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATg < 2:281656/141‑1 (MQ=255) aCCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCGCAAGCTGAAGCAAGAATGTTGATGtt > 2:175235/1‑141 (MQ=255) ttAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTAAAGCAAGAATGTTGATGTTAgccgca > 2:290172/1‑141 (MQ=255) aGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGcaca < 2:319994/141‑1 (MQ=255) tGAAGGTCGTGCGATTCGTTTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACaaaa < 2:29333/141‑1 (MQ=255) tGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACaaaa < 2:194706/141‑1 (MQ=255) aGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACAt > 1:85865/1‑141 (MQ=255) aGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACAt < 1:331519/141‑1 (MQ=255) ggTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATc < 1:124280/141‑1 (MQ=255) gTGCGATTCGTCTACATCCACTTGTAACAACCGATTATAAAGCTGACTTTGACGGTGCCCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCACGAATGTTGATGTTAGCAGCACAAAACATCTTGa > 1:74323/1‑141 (MQ=255) gTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGa > 1:68762/1‑141 (MQ=255) gTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGa > 1:150780/1‑141 (MQ=255) gATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAAccc < 1:238172/141‑1 (MQ=255) ttCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACACAACATCTTGAACCCTa > 1:101271/1‑141 (MQ=255) tGTAAAAACAGCGTATAACGCTTACTTTGACTGTGACCAAATGGCTGTTCACGTTCCTTTATCAATAGGAGTACACGCTGAAGCAAGAATGTTGATGTTAGCAGCATAACACATCTTGATCCCTAAAGATGGTAAACCTGt < 2:310781/141‑1 (MQ=255) gTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTa < 1:173662/141‑1 (MQ=255) acaGCTTATAACGCTCACTTTGACGGTGACCAAATGGCTGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTGAACCTGTAGTTaca < 2:270112/141‑1 (MQ=255) aGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACAcc < 1:212223/141‑1 (MQ=255) ataACGCTGAGTTTGACGATGACCAAATGGCGGGTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCACGAATGTTCATGTTAGCAGCACAAACCAGCTTGAACCCTAAATATGGCAAACCTGTAGTTACACCATcac > 2:132155/1‑141 (MQ=255) cTGACTTTGAGGGTGAAAAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAACTATGTTGATGTTAGTAGCACAAAACAGCTTGAACCCTAAAGATGGTAAACCTGTAGTAACACCATCACAAGata < 2:184867/141‑1 (MQ=255) cTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTCATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGt < 1:182878/141‑1 (MQ=255) gTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACt > 2:20588/1‑141 (MQ=255) ttCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACtt > 1:181500/1‑141 (MQ=255) ccTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTGGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACTTTagaaaga > 2:50873/1‑141 (MQ=255) tCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACTTTAGAAAGAAAAGAt > 2:229100/1‑141 (MQ=255) | TAGACTTGGTATTCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACTTTAGAAAGAAAAGAT > USA300TCH1516_ALE/591041‑591308 |
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: ‑ |