| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 6,086,626 | C→A | 19.1% | E25* (GAG→TAG) | PP_5337 ← | LysR family transcriptional regulator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 6,086,626 | 0 | C | A | 19.1% | 50.0 / 6.4 | 21 | E25* (GAG→TAG) | PP_5337 | LysR family transcriptional regulator |
| Reads supporting (aligned to +/- strand): ref base C (9/8); new base A (2/2); total (11/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 = 3.34e-01 | |||||||||||
GACGAGAAATAGCTGACCCGCCTCGGTCAGCTCGATGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTGCTTTCAAGGTTCACGAAGCACATCTCCCGGCGCGCACCAAAATGGAACACGCTCGAATATCAATTGGCGTCAC > NC_002947/6086495‑6086757 | gACGAGAAATAGCTGACCCGCCTCGGTCAGCTCGATGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCg > 2:556708/1‑136 (MQ=255) tAGCTGACCCGCCTCGGTCAGCTCGATGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGaa > 2:448052/1‑136 (MQ=255) ccGCCTCGGTCAGCTCGATGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCg > 1:451722/1‑136 (MQ=255) cGCCTGGGTCAGCTCGATGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCACCGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCgg < 1:373293/136‑1 (MQ=255) cGCCTCGGTCAGCTCGATGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCgg < 1:384510/136‑1 (MQ=255) aGCTCGATGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCa < 1:375800/136‑1 (MQ=255) tGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACt < 2:194401/136‑1 (MQ=255) cGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCtt > 1:234349/1‑136 (MQ=255) cGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCtt > 1:392731/1‑136 (MQ=255) ttCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTATGCGGCCt < 1:498552/88‑1 (MQ=255) ttCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTATGCGGCCt > 2:498552/1‑88 (MQ=255) ttCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTATGCGGCCt > 2:495732/1‑88 (MQ=255) ttCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTATGCGGCCt < 1:495732/88‑1 (MQ=255) cACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTg > 2:504435/1‑136 (MQ=255) aCCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTGc > 2:537719/1‑136 (MQ=255) gCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTGCTTTCAAGGTTCACGAAGCACATCTCCCGGCGCGCACCAAAATgg > 2:303364/1‑136 (MQ=255) gggTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTGCTTTCAAGGTTCACGAAGCACATCTCCCg > 2:222556/1‑106 (MQ=255) gggTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTGCTTTCAAGGTTCACGAAGCACATCTCCCg < 1:222556/106‑1 (MQ=255) aaCGCCGCTCTGCTGGCTGTGAAAAACTGCGGGTAGCGGCCATAGCACTGAAGTATTCCAGGCATTTGCTTTATAGGTTCAAGAATCACATCTGCCGGCGCGCTGCAAATTGCAACTCGCTCGAATATCAAttttc < 2:97098/136‑4 (MQ=255) aaCGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTGCTTTCAAGGTTCACGAAGCACATCTCCCGGCGCGCACCAAAATGGAACACGCTCGAATATCAATTGGc < 1:448052/136‑1 (MQ=255) cgcTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTGCTTTCAAGGTTCACGAAGCACATCTCCCGGCGCGCACCAAAATGGAACACGCTCGAATATCAATTGGCGTcac < 1:527241/136‑1 (MQ=255) | GACGAGAAATAGCTGACCCGCCTCGGTCAGCTCGATGGGTGTGCGAGAGCGGTTCACCAGTTGCAGGCCCAACGCCGCCTCAAGGCTGCGAATTCGCCGGCTGAAGGCCGGTTGGGTCACGAAACGCCGCTCTGCGGCCTGGGAAAAACTGCGGGTAGCGGCCAGAGCACTGAAGTCTTCCAGCCATTTGCTTTCAAGGTTCACGAAGCACATCTCCCGGCGCGCACCAAAATGGAACACGCTCGAATATCAATTGGCGTCAC > NC_002947/6086495‑6086757 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |