Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 27,581 | C→T | 15.4% | intergenic (‑428/+18) | PP_0022 ← / ← PP_0023 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 27,581 | 0 | C | T | 15.4% | 84.6 / 9.2 | 39 | intergenic (‑428/+18) | PP_0022/PP_0023 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (17/16); new base T (4/2); total (21/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01 |
TCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATG > NC_002947/27440‑27728 | tcgtcgTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAATCAGCTCCCGAGAgct > 1:262915/1‑150 (MQ=255) acggAAGGGATAACCCATCCGGAGGTAGGTTCAGAATCCACCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGGGAGACGGTCGTTAGGGTTTGCAAACAGCCCCAAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTc < 1:343003/147‑1 (MQ=255) aaGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGTCTCACGACGCGGTCCGATTGATCAGCGAGAGGGTCGTGAGGGTTTGCACAAAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTc < 1:436755/150‑1 (MQ=255) ccGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCA‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCtttt > 1:40361/1‑150 (MQ=255) cGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTtg > 1:171723/1‑150 (MQ=255) gAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTtgt < 1:532571/150‑1 (MQ=255) ttCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCttttt > 1:482376/1‑134 (MQ=255) ttCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCttttt < 2:482376/134‑1 (MQ=255) tCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCA‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTc > 1:126836/1‑150 (MQ=255) aaTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTa < 2:403589/117‑1 (MQ=255) aaTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTa > 1:403589/1‑117 (MQ=255) tccagccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGCCGGTCGTCAGGGATTGCACACAGCCCC‑‑AAAAAAAAAGCATCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATg < 1:169808/147‑1 (MQ=255) aagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc > 2:422592/1‑150 (MQ=255) gAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCtt < 1:15267/150‑1 (MQ=255) cTCACGACGCGGTCCTATTGAGCAGCGAGACGGTCGTGAGGGTTTGCACACAGCGCC‑‑GAAAAGAAAGTAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAg < 2:552714/150‑1 (MQ=255) acgacgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCA‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg < 2:126836/150‑1 (MQ=255) acgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgg > 2:220790/1‑145 (MQ=255) acgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgg < 1:220790/145‑1 (MQ=255) gTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATg > 1:594531/1‑122 (MQ=255) gTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATg < 2:594531/122‑1 (MQ=255) aGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAg < 1:519927/101‑1 (MQ=255) aGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAg > 2:519927/1‑101 (MQ=255) cGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATc < 1:422592/150‑1 (MQ=255) tttGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc > 1:561490/1‑98 (MQ=255) tttGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc < 2:561490/98‑1 (MQ=255) acacaGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTgcgc > 1:360449/1‑150 (MQ=255) acaGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCCTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCtt > 2:351036/1‑150 (MQ=255) ccca‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTggg > 1:452759/1‑150 (MQ=255) taagaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTTTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTgg < 2:708961/142‑1 (MQ=39) taagaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTTTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTgg > 1:708961/5‑146 (MQ=39) taagaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTTTTTTTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCt > 2:505986/5‑150 (MQ=255) agaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTgg > 1:112647/3‑150 (MQ=255) aaaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTgg > 1:158843/1‑150 (MQ=255) aaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGTCTGGc > 1:555197/1‑150 (MQ=255) aaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTTCGCTTTGGGCCTGGCTTGa > 2:282587/1‑150 (MQ=255) aaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGAt < 2:171723/150‑1 (MQ=255) agcTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTc < 2:112647/150‑1 (MQ=255) cTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCAt < 2:360449/150‑1 (MQ=255) tCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGGGGGGTTTCAGGTCCGGATGGGCCTTCTGCGTCTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCGTGATGCTCGTg > 1:100969/1‑150 (MQ=255) ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGTTGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg > 1:548968/1‑79 (MQ=255) ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg < 2:548968/79‑1 (MQ=255) ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGc > 1:674831/1‑91 (MQ=255) ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGc < 2:674831/91‑1 (MQ=255) | TCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATG > NC_002947/27440‑27728 |
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: ‑ |