| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,969,261 | G→T | 20.1% | L343I (CTA→ATA) | fliC ← | flagellin |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,969,261 | 0 | G | T | 20.1% | 39.5 / 3.6 | 20 | L343I (CTA→ATA) | fliC | flagellin |
| Reads supporting (aligned to +/- strand): ref base G (9/7); new base T (2/2); total (11/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01 | |||||||||||
GACCCAGGCCAGTCAAGGTACCCGCAGTACCACCTACGTCGAAGCCGGCAGCCGAGTTCAGTTTCAGTTGCCCAGTGGCACTGTCAATGGAAGCCGTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCGACAGCTGGGCACCACCGGCCTCGGCCAGACCCAAGGCA > NC_002947/4969114‑4969403 | gACCCAGGCCAGTCAAGGTACCCGCAGTACCACCTACGTCGAAGCCGGCAGCCGAGTTCAGTTTCAGTTGCCCAGTGGCACTGTCAATGGAAGCCGTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAggg < 2:434517/150‑1 (MQ=255) tCAAGGTACCCGCAGTACCACCTACGTCGAAGCCGGCAGCCGAGTTCAGTTTCAGTTGCCCAGTGGCACTGTCAATGGAAGCCGTTACACCTGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGTAGATAGGGGGTCAGTCGCac > 2:77136/1‑149 (MQ=255) aGGTACCCGCAGTACCACCTACGTCGAAGCCGGCAGCCGAGTTCAGTTTCAGTTGCCCAGTGGCACTGTCAATGGAAGCCGTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAAAca < 1:713333/150‑1 (MQ=255) tcagtttaatttGCCCAGTCACACGGTCGATGGAAGCCGTTACACCAGTCGTGCCGGTTTGAG‑ATTAGCCCTTCCCGACGACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGcaccac < 1:49232/150‑1 (MQ=255) tGCCCAGTGGCACTGTCAATGGAAGCCGTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAAt < 1:191054/150‑1 (MQ=255) ccTGTCAATGGAAGCCGTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTaga > 1:751874/2‑102 (MQ=255) aCTGTCAATGGAAGCCGTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTaga < 2:751874/102‑1 (MQ=255) gTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAg > 1:588606/1‑98 (MQ=255) gTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAg < 2:588606/98‑1 (MQ=255) cAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATGTGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCg > 1:15498/1‑150 (MQ=255) cAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATATGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCg > 1:414890/1‑150 (MQ=255) cAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATac > 1:588352/1‑118 (MQ=255) cAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATac < 2:588352/118‑1 (MQ=255) ttGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCCACACCGGCCCCCAATACACCATCACGCAGGTTGACAGTACATGTGTTGGACAGCTGGGCACCACCGGCCTCGGcc > 1:650597/1‑150 (MQ=255) aCTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCGACAGCTGGGCACCACCGGCCTCGGCCAGAccc > 2:592065/1‑150 (MQ=255) ccGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCaa < 1:555000/79‑1 (MQ=255) ccGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCaa > 2:555000/1‑79 (MQ=255) ccGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCGACAGCTGGGCACCACCGGCCTCGGCCAGACCCAAgg > 1:773171/1‑150 (MQ=255) ccGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCGACAGCTGAGCACCACCGGCCTCGGCCAGACCCAAgg > 2:133504/1‑150 (MQ=255) gAAGATGTGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCGACAGCTGGGCACCACCGGCCTCGGCCAGACCCAAGGCa < 2:15498/150‑1 (MQ=255) gAAGATATGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCGACAGCTGGGCACCACCGGCCTCGGCCAGACCCAAGGCa < 2:414890/150‑1 (MQ=255) | GACCCAGGCCAGTCAAGGTACCCGCAGTACCACCTACGTCGAAGCCGGCAGCCGAGTTCAGTTTCAGTTGCCCAGTGGCACTGTCAATGGAAGCCGTTACACCAGTCGTGCCGGTTTGAGCATTA‑ACCTTGCCGACTACCGAAGATAGGGTGTCAGTCGCAACCAATGCAATGTTGGTAGAGCCGAGCGTCAGGTTACCGGCACCACCGGCCCCCAATACACCATCACGCAGGTTGGCAGTACCGGTTTTCGACAGCTGGGCACCACCGGCCTCGGCCAGACCCAAGGCA > NC_002947/4969114‑4969403 |
| 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: ‑ |