Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 932,321 | C→T | intergenic (‑271/‑274) | trxB ← / → lrp | thioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 932,321 | 0 | C | T | 95.0% | 68.5 / ‑4.0 | 20 | intergenic (‑271/‑274) | trxB/lrp | thioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (6/13); total (7/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.50e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01 |
TCACCTCTGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCTCAC > NC_000913/932187‑932468 | tCACCTCTGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaa > 1:161954/1‑151 (MQ=255) ccTCTGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGaaaaaaaa > 1:18001/1‑151 (MQ=255) cTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCt < 1:235488/151‑1 (MQ=255) tttCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATg < 1:226476/151‑1 (MQ=255) tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAAc < 1:153711/151‑1 (MQ=255) tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAAc < 1:201192/151‑1 (MQ=255) tataAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACca < 2:3761/151‑1 (MQ=255) aaGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACcatgc < 2:162310/151‑1 (MQ=255) aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc < 1:219840/151‑1 (MQ=255) aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc < 2:138685/151‑1 (MQ=255) aaTGATGAATAAACGCCCCTGTTAATGAACATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc < 2:17182/151‑1 (MQ=255) tAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa > 1:250704/1‑150 (MQ=255) gTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa > 2:64091/1‑151 (MQ=255) gTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa > 1:48364/1‑151 (MQ=255) ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGAc < 1:164928/151‑1 (MQ=255) aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTg < 1:147737/150‑1 (MQ=255) gTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATg > 1:242523/1‑151 (MQ=255) tttGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAg < 2:83171/151‑1 (MQ=255) aTCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCca > 2:220686/1‑150 (MQ=255) aaaTATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCTCAc < 1:177454/151‑1 (MQ=255) | TCACCTCTGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCTCAC > NC_000913/932187‑932468 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |