Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,018,201 | G→A | 100% | intergenic (+144/‑36) | mvaB → / → PP_mr44 | hydroxymethylglutaryl‑CoA lyase/CrcZ ncRNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,018,201 | 0 | G | A | 95.5% | 58.7 / ‑3.8 | 22 | intergenic (+144/‑36) | mvaB/PP_mr44 | hydroxymethylglutaryl‑CoA lyase/CrcZ ncRNA |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (10/11); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TGTCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAG > NC_002947/4018068‑4018326 | tgtCACATTGGTGTGGGTGAAGTGTTACCTAACGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTATAAAGGAtt > 2:491669/1‑136 (MQ=255) tGGTGTGGGTGAAGTGTTACCTCCAGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAATGGATTTTCAAAAg > 1:362728/1‑136 (MQ=255) gTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacg > 2:516682/1‑136 (MQ=255) tGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgc < 1:195423/136‑1 (MQ=255) gAAGTTTTACCTAGCGGACTTTTCAGAGCAAAAAATCGGGTAACAGGGAAACATATCGGCAGATTTTCATGCTGGGTGGTTTCGACCAAAAATCTTAACTCATTATTTTTAAAGGATTTTCAAAAGTTGgcacgca < 1:491669/136‑1 (MQ=255) gAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgca < 2:121199/136‑1 (MQ=255) aGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCAcc < 2:314281/136‑1 (MQ=255) cTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCtata > 1:416670/1‑136 (MQ=255) acaGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTacaaca > 2:130249/1‑136 (MQ=255) acaGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTacaaca > 2:343928/1‑136 (MQ=255) aaaTCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACaaaaaa > 2:137155/1‑136 (MQ=255) cggAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCATAATTCTTATCTCATTTATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTTCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCa < 2:362728/135‑1 (MQ=255) gAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAAt < 2:519873/136‑1 (MQ=255) aaCAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATaa < 2:98692/136‑1 (MQ=255) cAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAgg < 2:308408/136‑1 (MQ=255) cAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAgg > 2:300924/1‑136 (MQ=255) ggtggtGTCGACCAAAAATCTTAACTCATTAATTTTAAAGTATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACAt < 1:259736/136‑1 (MQ=255) gtggtgTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATa > 2:510111/1‑136 (MQ=255) aCCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGaaca < 2:386851/136‑1 (MQ=255) aCCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGATATATCTCTGTTACAACAACAAAAAAAACGCTATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGaaca < 2:456227/136‑1 (MQ=255) aCCAAAAATATTAACTAATTAATTTTAAAGGATTTTCAAAATTTGGAACGCACCCTGATATATCTCTGGTAAAACAAAAAAAAAAACGCTATACCCAATAAAAAAAACAGGTAACGGCTCTGACATAACAAGaaca < 2:180690/136‑1 (MQ=255) aTTTTAAAGGATTTTCAAAAGTTCGCACGCACCCTGCTATATCTCAGGTACAACAACAAAAAAAACGCGATACCCAAAAAAAACAACAGGAAACGGCTCTGACATAACAACAACCACACGGCAGACGCGAaccaac > 1:59529/1‑135 (MQ=255) | TGTCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAG > NC_002947/4018068‑4018326 |
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: ‑ |