Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,843,140 | G→A | 100% | Q566* (CAG→TAG) | gacS ← | sensor protein GacS |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,843,140 | 0 | G | A | 96.2% | 78.0 / ‑4.2 | 26 | Q566* (CAG→TAG) | gacS | sensor protein GacS |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (13/12); total (13/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGG > NC_002947/1842993‑1843283 | tCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTaca < 1:820900/150‑1 (MQ=255) cacaGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAAc < 2:842174/150‑1 (MQ=255) cacaGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAAc < 1:103260/150‑1 (MQ=255) aaaGCACTACGCATTGGCAGTGGTAACGCTCGATTCCCCGATGATCGTGGCTCGGTTGCTCGGGCGATAGGTTACCCCGGTTGGCACCGCGAACCGCGCACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACTCCCTACAGCAACt < 2:240769/148‑1 (MQ=255) cGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTagg < 1:696711/150‑1 (MQ=255) aTTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTagga > 1:306064/1‑148 (MQ=255) aTTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTagga < 2:306064/148‑1 (MQ=255) gCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGACCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGc > 2:389529/1‑150 (MQ=255) atAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTc < 1:131137/150‑1 (MQ=255) gggCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGc > 1:31219/1‑122 (MQ=255) gggCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGc < 2:31219/122‑1 (MQ=255) gACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCAAACCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTgg > 1:62081/1‑150 (MQ=255) gTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGgcgcgc > 2:293130/1‑150 (MQ=255) ccAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGtt > 2:772461/1‑150 (MQ=255) ccAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGtt > 1:327979/1‑150 (MQ=255) gTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTg < 2:327979/150‑1 (MQ=255) ttGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTgg < 1:293130/150‑1 (MQ=255) cgcgAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATggc > 2:379037/1‑150 (MQ=255) tCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcgg < 2:557116/150‑1 (MQ=255) gAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTTCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcggcg > 1:401790/1‑150 (MQ=255) gCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTg > 1:521146/1‑91 (MQ=255) gCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTg < 2:521146/91‑1 (MQ=255) cTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGt > 1:474035/1‑142 (MQ=255) cTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGt < 2:474035/142‑1 (MQ=255) cccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTAACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCgg > 1:629911/1‑150 (MQ=255) ccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTgg > 1:753951/1‑138 (MQ=255) ccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTgg < 2:753951/138‑1 (MQ=255) | TCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGG > NC_002947/1842993‑1843283 |
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: ‑ |