Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 4,117,071 | T→C | F44S (TTC→TCC) | rpmE → | 50S ribosomal protein L31 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,117,071 | 0 | T | C | 82.8% | 54.9 / 9.6 | 29 | F44S (TTC→TCC) | rpmE | 50S ribosomal protein L31 |
Reads supporting (aligned to +/- strand): ref base T (4/1); new base C (12/12); total (16/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.43e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01 |
AACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATA > NZ_CP009273/4116998‑4117152 | aaCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGt < 2:50376/90‑1 (MQ=255) aCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTg > 1:229426/1‑90 (MQ=255) aTGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGtt < 2:324101/90‑1 (MQ=255) gAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGc < 1:179156/90‑1 (MQ=255) gCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTg > 1:783883/1‑90 (MQ=255) cTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTgg < 2:50748/90‑1 (MQ=255) gTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGtt > 1:766775/1‑90 (MQ=255) tCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCg > 2:619907/1‑90 (MQ=255) tCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCg < 1:216942/90‑1 (MQ=255) tCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCg < 2:447495/90‑1 (MQ=255) aTGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCt > 2:415167/1‑90 (MQ=255) tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt > 2:665322/1‑90 (MQ=255) tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt > 2:526668/1‑90 (MQ=255) tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt > 1:30473/1‑90 (MQ=255) gACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTc > 1:386307/1‑90 (MQ=255) gACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGc < 2:288141/63‑1 (MQ=255) gACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGc > 1:288141/1‑63 (MQ=255) cGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCaa < 1:415167/90‑1 (MQ=255) gtgtGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAAc > 1:99990/1‑90 (MQ=255) gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa < 2:23012/78‑1 (MQ=255) gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa > 1:23012/1‑78 (MQ=255) gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACAt > 1:321409/1‑90 (MQ=255) aaGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgc < 1:748209/90‑1 (MQ=255) aaGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgc < 1:150230/90‑1 (MQ=255) aGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgca > 2:100874/1‑90 (MQ=255) gTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCTGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcag < 1:101167/90‑1 (MQ=255) tGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcagc > 1:525889/1‑90 (MQ=255) tGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcagc > 1:125059/1‑90 (MQ=255) cccGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATa < 2:481197/90‑1 (MQ=255) | AACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATA > NZ_CP009273/4116998‑4117152 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |