Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 646,747 | 0 | T | G | 34.6% | 33.2 / 13.2 | 26 | T255P (ACC→CCC) | citG | triphosphoribosyl‑dephospho‑CoA synthase |
Reads supporting (aligned to +/- strand): ref base T (1/16); new base G (9/0); total (10/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.20e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCC > NC_000913/646619‑646880 | gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑cccc < 1:502804/140‑1 (MQ=255) cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCg > 2:317577/1‑129 (MQ=255) cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCg < 1:317577/129‑1 (MQ=255) cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑cccccct > 1:626658/1‑138 (MQ=255) atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCCTTTTTCaaa > 1:333919/1‑139 (MQ=255) tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTTGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCaataa < 1:653435/140‑1 (MQ=255) tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCAATAATGttttcc > 1:617181/1‑137 (MQ=255) tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCAATAATGTTTtcg > 2:18453/1‑140 (MQ=255) gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCCTTTTTCAATAAAGtttttgc > 2:690903/1‑136 (MQ=255) aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCa < 2:78276/140‑1 (MQ=255) aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg < 1:97246/140‑1 (MQ=255) tCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCCGATT‑TCCCCTTTTTTTAAAAAAATTTTCGCCCCCCGCCGTAGGCAgcgc > 1:751436/1‑140 (MQ=255) aggTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCa < 1:512235/138‑1 (MQ=255) aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAG‑CCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCCGTAGCCAGCGCAg > 1:563721/1‑140 (MQ=255) gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCCTTTTTCAATAATATTTTCGCCTCCCGCTGTAGCCAGCGCAgg > 1:639141/1‑140 (MQ=255) gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑GCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTAGCCAGCGCAGGcccc > 1:485728/1‑140 (MQ=255) cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCAATAAATTTTTCGCCCCCCGCTGTAGCCCGCGCAGGccccccccc > 1:30109/1‑137 (MQ=255) cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATG‑CCCCCCTTTTTGAATAATTTTTTCGCCCCGCGCTGTAGCCAGCGCACGccccccccccc > 2:136185/1‑135 (MQ=255) cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATG‑GCCCCCTTTTTTTATTATGTTTTAGGCCCCCCCTGTTGGCAGCGGCGGCCCCCCCCGcc > 2:512235/1‑140 (MQ=255) ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGa < 2:766757/140‑1 (MQ=255) cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc < 2:161485/140‑1 (MQ=255) tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCAAG‑CCCCCCTTTTTCAATAATATTTTCGCCCCCCCCCGTAGCCAACGCAGGGCCCCCCCGCCACCCCGCGATGCCACCTTGGGATTGCCGTTGAt > 2:55712/1‑140 (MQ=255) tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGcc < 1:257818/140‑1 (MQ=255) gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCCTTTTTTAATAAAATTTTCGCCCCCCCCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCCACGGTTGTATCGCCGCTGATTGccc > 2:257818/1‑139 (MQ=255) tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt < 1:136185/140‑1 (MQ=255) aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt < 1:589898/140‑1 (MQ=255) ccGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGa < 1:95881/140‑1 (MQ=255) gagaTAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa < 2:626658/140‑1 (MQ=255) gagaTAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa < 1:18453/140‑1 (MQ=255) gagaTAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCACCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa < 1:676835/140‑1 (MQ=255) aaTCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTa < 2:333919/140‑1 (MQ=255) agaTCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAGCAGTAGGAGCAAGGTATCGAg < 2:751436/140‑1 (MQ=255) aTCGGCGGGGGGTCGAATG‑CCCCCCTTTTTTAATAAAGTTTTTGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGGTGGTATCGCCGTTGGTCGCCCTCAGTAGGAGCAAGGGATCGAGCa > 2:397304/1‑140 (MQ=255) cggcggGGGGTCGGATGCCCCCCTTTTTCCAAAAAT‑TTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCACCCCCCGGTGCAACGGTGGTGTTGCCGGTTGTCGCCCCCCGTAGGAGCAAGGGATCGAGCaaa > 2:296472/1‑139 (MQ=255) cggGGGGTCGAAGG‑CCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGTGCAACGTTGGTATCGCCGTTGATCGCCCTCAGTAGGAGCAAGGGATCGAGCAAAGcc > 2:234784/1‑140 (MQ=255) | GTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCC > NC_000913/646619‑646880 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |