Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_pae_tpiA | 4,158,703 | 0 | A | G | 33.3% | 75.9 / 55.5 | 64 | T100A (ACA→GCA) | oxyR | oxidative and nitrosative stress transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base A (23/19); new base G (10/11); total (33/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGT > NC_000913_3_pae_tpiA/4158558‑4158851 | cAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGt > 4:51943/1‑150 (MQ=255) gCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTg > 4:253344/1‑150 (MQ=255) tgttgtTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTAcc < 3:51943/150‑1 (MQ=255) gttgttCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCt > 2:722869/1‑150 (MQ=255) aCCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCg < 4:149705/138‑1 (MQ=255) aCCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCg > 3:149705/1‑138 (MQ=255) ccAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCa > 3:284738/1‑150 (MQ=255) gAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATcc < 4:324025/150‑1 (MQ=255) gtggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCat < 1:4226/132‑1 (MQ=255) gtggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCat > 2:4226/1‑132 (MQ=255) gATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAg > 4:489843/1‑150 (MQ=255) aTCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATcc < 3:541596/135‑1 (MQ=255) aTCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATcc > 4:541596/1‑135 (MQ=255) aGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTAt < 3:524771/136‑1 (MQ=255) aGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTAt > 4:524771/1‑136 (MQ=255) gcgcGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCttt > 2:869229/1‑150 (MQ=255) cgcgTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTc < 2:457166/150‑1 (MQ=255) cgTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCa > 1:973128/1‑150 (MQ=255) aCCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAg > 1:61322/1‑150 (MQ=255) cGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCt < 1:927542/150‑1 (MQ=255) tgcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGaaa < 4:173920/150‑1 (MQ=255) gcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAAt > 2:666147/1‑150 (MQ=255) gcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCCCAGTCGGACCGTTCCTGCTACCGCATATTATCCCTCTGCTGGACCAGAAATTTCAAACGCTGGAAAt < 2:441526/150‑1 (MQ=255) gTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGt < 4:478894/150‑1 (MQ=255) gTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGt < 3:489843/150‑1 (MQ=255) aGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGa < 1:566540/126‑1 (MQ=255) aGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGa > 2:566540/1‑126 (MQ=255) aGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATc < 3:253344/150‑1 (MQ=255) aGGTGAAAGTCCTGAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATc > 2:461117/1‑150 (MQ=255) ggTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCTGACCTTTCCAAAGCTGGAAATGTATCt < 1:722869/150‑1 (MQ=255) ttAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCatat < 1:320532/94‑1 (MQ=255) ttAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCatat > 2:320532/1‑94 (MQ=255) aaaGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGc > 3:162909/1‑124 (MQ=255) aaaGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGc < 4:162909/124‑1 (MQ=255) agagATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGAcc < 2:731425/112‑1 (MQ=255) agagATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGAcc > 1:731425/1‑112 (MQ=255) agagATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGAc > 2:611121/1‑150 (MQ=255) aTGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAACGCTGGAAATGTATCTGCATGAAGCACAGACccac > 1:457266/1‑150 (MQ=255) aTGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccac > 3:439388/1‑150 (MQ=255) ggCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccacca < 2:61322/150‑1 (MQ=255) aaGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGtt > 2:825526/1‑150 (MQ=255) cagcagGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTg < 2:973128/150‑1 (MQ=255) agcagGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTgg > 2:576836/1‑150 (MQ=255) gCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAAc < 2:457266/150‑1 (MQ=255) cGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCttt > 3:565164/1‑94 (MQ=255) cGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCttt < 4:565164/94‑1 (MQ=255) gagaCGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCata > 1:116829/1‑67 (MQ=255) gagaCGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCata < 2:116829/67‑1 (MQ=255) tGTCCGGACCGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCg < 1:666147/150‑1 (MQ=255) tGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCg < 4:439388/150‑1 (MQ=255) ccGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCa > 3:329437/1‑150 (MQ=255) ggACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCaaa < 3:29131/150‑1 (MQ=255) cGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCaa > 1:850916/1‑81 (MQ=255) cGCTGCACATTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCaa < 2:850916/81‑1 (MQ=255) cGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCg > 2:501307/1‑150 (MQ=255) tGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCt > 4:331776/1‑94 (MQ=255) tGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCt < 3:331776/94‑1 (MQ=255) gCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACcca > 2:49878/1‑111 (MQ=255) gCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACcca < 1:49878/111‑1 (MQ=255) acaTTGGTTTGATTCCCGCAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCg < 1:576836/150‑1 (MQ=255) tttGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCa > 2:595560/1‑54 (MQ=255) tttGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCa < 1:595560/54‑1 (MQ=255) cccACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTg > 1:9742/1‑150 (MQ=255) cccACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTg > 1:10057/1‑150 (MQ=255) cacaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCAc < 1:624928/46‑1 (MQ=255) cacaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCAc > 2:624928/1‑46 (MQ=255) cacaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCa < 4:200871/48‑1 (MQ=255) cacaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCa < 4:200860/48‑1 (MQ=255) cacaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCa > 3:200871/1‑48 (MQ=255) cacaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCa > 3:200860/1‑48 (MQ=255) cacaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGt < 2:10057/150‑1 (MQ=255) cacaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGt < 2:9742/150‑1 (MQ=255) | CAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGT > NC_000913_3_pae_tpiA/4158558‑4158851 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |