Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 555,771 | C→A | 16.7% | H29N (CAC→AAC) | rplN → | 50S ribosomal protein L14 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 555,771 | 0 | C | A | 16.7% | 54.7 / 5.0 | 24 | H29N (CAC→AAC) | rplN | 50S ribosomal protein L14 |
Reads supporting (aligned to +/- strand): ref base C (12/8); new base A (2/2); total (14/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01 |
TGGTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTCCATCATTCG > NC_002947/555629‑555919 | tGGTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCAccgc > 1:335075/1‑148 (MQ=255) tGGTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCAccgc < 2:335075/148‑1 (MQ=255) gTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAg > 1:394768/1‑150 (MQ=255) taagggctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTaa > 2:415083/1‑150 (MQ=255) gctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTc > 1:550359/1‑150 (MQ=255) taagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTCCTCGGCGGTTCGCACCGCCGTTACTCGGGTAGCGGTGAGATCATCAAAGGAACCGGCaa < 2:315761/150‑1 (MQ=255) ggAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCaa < 1:575219/150‑1 (MQ=255) agaAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAAtt < 2:66002/150‑1 (MQ=255) acAATCCATGCTCGATGTGGCCGATAACATCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCTCCGTTACGCGCGTATGGGTGACATCATAAAAGTAACCGTCAAGGAAGAAATTCAGGTCGGTAAGGTaaaaaa > 2:33371/2‑150 (MQ=255) aaTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAg < 1:196636/123‑1 (MQ=255) aaTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAg > 2:196636/1‑123 (MQ=255) aTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAAcc > 1:10944/1‑112 (MQ=255) aTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAAcc < 2:10944/112‑1 (MQ=255) tgtgGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGAc > 1:106137/1‑150 (MQ=255) gCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCt < 2:394768/150‑1 (MQ=255) tcgCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCa > 2:211011/1‑150 (MQ=255) aTGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGtt > 1:466182/1‑150 (MQ=255) aTCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTcgc > 2:18251/1‑150 (MQ=255) aacggcggTGCGAACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGa < 1:449104/144‑1 (MQ=255) aacggcggTGCGAACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGa > 2:449104/3‑146 (MQ=255) acggcggTGCGAACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAg < 1:592294/79‑1 (MQ=255) acggcggTGCGAACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAg > 2:592294/2‑80 (MQ=255) cggcggTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTc < 2:466182/150‑1 (MQ=255) gCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTCCATCATTCg > 2:91419/1‑150 (MQ=255) | TGGTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTCCATCATTCG > NC_002947/555629‑555919 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |