Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,981,203 | T→C | 21.1% | intergenic (+99/‑201) | rpsA → / → ihfB | 30S ribosomal protein S1/integration host factor subunit beta |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,981,203 | 0 | T | C | 21.1% | 27.7 / 4.7 | 19 | intergenic (+99/‑201) | rpsA/ihfB | 30S ribosomal protein S1/integration host factor subunit beta |
Reads supporting (aligned to +/- strand): ref base T (9/6); new base C (2/2); total (11/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01 |
ACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTC > NC_002947/1981057‑1981335 | accaccATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTGGATCGGTAAAAAGGGCGACCCTCGGGTCGCGCTTTTTTTTGCTCCACGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCa > 2:13379/1‑150 (MQ=255) caccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAgg > 1:150522/1‑150 (MQ=255) aaTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGt < 2:150522/150‑1 (MQ=255) gtgttaGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc > 1:238776/6‑150 (MQ=255) cAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGTCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc < 2:202936/150‑1 (MQ=255) cAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc > 2:358903/1‑150 (MQ=255) cAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc < 2:311212/150‑1 (MQ=255) aacAGAACTCAGTTCTGCTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACAAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcgc > 2:5600/3‑150 (MQ=255) aTCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTa > 1:280753/1‑150 (MQ=255) ggCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgtgtt > 2:88646/1‑88 (MQ=255) ggCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgtgtt < 1:88646/92‑5 (MQ=255) ggCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCg < 2:249430/88‑1 (MQ=255) ggCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCg > 1:249430/1‑88 (MQ=255) tttttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAAt < 1:358903/150‑1 (MQ=255) gtggtgCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTACATGCTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTaa < 2:238776/150‑1 (MQ=255) aCTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCTAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGt < 1:5600/150‑1 (MQ=255) aCTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGt > 1:50719/1‑150 (MQ=255) aCTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGt > 1:274020/1‑150 (MQ=255) tgctgGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTc > 1:443116/1‑150 (MQ=255) | ACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTC > NC_002947/1981057‑1981335 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |