Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 623,110 | C→A | 23.6% | intergenic (+59/+103) | PP_0535 → / ← PP_0536 | ferric siderophore receptor/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 623,110 | 0 | C | A | 23.6% | 25.9 / 5.1 | 17 | intergenic (+59/+103) | PP_0535/PP_0536 | ferric siderophore receptor/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (5/8); new base A (2/2); total (7/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 = 4.05e-01 |
GGATTACCTGTCGCGGGTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCAAGAAGTCTCCGA > NC_002947/622961‑623227 | ggATTACCTGTCGCGGGTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGc > 1:467166/1‑150 (MQ=255) tGTCGCGGGTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGc < 2:442865/150‑1 (MQ=255) gTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCAt < 1:6173/150‑1 (MQ=255) cTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGaa > 2:352287/1‑150 (MQ=255) tCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGaaa > 1:435857/1‑150 (MQ=255) ccACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGaa < 1:350150/150‑1 (MQ=255) cagtaTGATGGTTGGGGCGGGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAACCATGAGAAATCTAGAAGGGTATCTGCACTATCGGGTATCAGTGc < 2:406098/150‑1 (MQ=255) ggTTGGAGCGCGGTTTGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACTGCG‑C‑CTTTGCTGAACAGGATAAAAGGTCGCGGCCTTGATTCTGAATCATGAAAAATCTAGAAATGTATCTGCACTATCGGGTATCAGTGCAATAAAGtt < 1:369409/150‑1 (MQ=255) aCGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGACTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATaaa < 1:121678/123‑1 (MQ=255) aCGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGACTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATaaa > 2:121667/1‑123 (MQ=255) aCGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGACTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATaaa > 2:121678/1‑123 (MQ=255) aCGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGACTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATaaa < 1:121667/123‑1 (MQ=255) tGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGttt > 1:465955/1‑150 (MQ=255) gCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCa > 1:378224/1‑150 (MQ=255) cTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCaa < 2:435857/150‑1 (MQ=255) gcgACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCAAGAAGTCTCCgg < 1:352287/150‑2 (MQ=255) gTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAg > 1:140876/1‑97 (MQ=255) gTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAg < 2:140876/97‑1 (MQ=255) | GGATTACCTGTCGCGGGTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCAAGAAGTCTCCGA > NC_002947/622961‑623227 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |