Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,061,105 | 0 | A | G | 52.2% | 10.0 / 27.4 | 23 | L511P (CTG→CCG) | yihQ | sulfoquinovosidase |
Reads supporting (aligned to +/- strand): ref base A (7/4); new base G (0/12); total (7/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.35e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCC > NZ_CP009273/4060987‑4061230 | gCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATg > 2:241340/1‑139 (MQ=255) ccGCCAAACGGCAGGTTGTCGCCAGGACGTTTCCCTTGGTGGGTGCGCATCACCGGCGTGAAGGCGCTGAAACCGCACCGGCGCAGCAGCAGCTCTTTGCTGCGTTTCATCTCAAACGGGGTGGTGTAACCGCCAATGt < 2:263792/139‑1 (MQ=255) aatcTGCCATTTGTCGCAAGGACTTTTCCCTTCGGGGGTGCGCATCATCGGCGTAAAGGCGCTGAAACCGCCCCAGCGCAGCAGCGCCTTTTTCCTGCTCTTCATCTAAAACGGGGTGGTGAAACCGCCAATGTCGCtg < 1:614598/136‑1 (MQ=255) aaccTGCCAGTGGTCGCCAGGGCGTTTCCCTTCGGGGGGGCGCATCATCGGCGTGAAGGCGCTGAAACCGCCCCGGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACGGGGTGGTTTATCCGCCAATGTCGCtg < 1:501216/136‑1 (MQ=255) aaaCTGCCGGTGGTCCCCAGGACGTTTCCCTTCGGGGGTGCGCATCATCGGCGTGAGGGCGTTAAAACCGCCCCAGCGCAGCAGCAGCTTTTTGCTGCGTTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCtg < 1:493220/139‑1 (MQ=255) ccAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCCCCAGCGCAGCAGCAGCTTTTTGCTGCGTTTCATCTAAAACAGGGGGGTGTAACCGCCAATGTCGCTgtggtg < 1:530782/138‑1 (MQ=255) ttGTCGCCGGGACGGTTCCCTTGGGGGGGGCGCATCATCGGCGTGAAGGCGCTAAAACCGCCCCAGCGCAGCAGCAGCTTTTTGCTGCGTTTCATCTCAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGc < 1:381481/139‑1 (MQ=255) ccAGGACGGTTCCCTTGGGGGGTGCGCATCATCGGGGTGAAGGCGCAAAAACCGCCCGGGGGCGGCAGGAGCTTTTGGCTGCGCTTCATCTCAAGCGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc < 2:285281/139‑1 (MQ=255) ccAGGACGGTTCCCCTCGTGGGTGAGCACCATCGGGGTGAGGGCGGTGAAATCGCCCCAGCGCAGCACCAGCTTTTTGCTGCGCTTCATCCCAAACAGGGTGGTATAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc < 2:504206/139‑1 (MQ=255) gagggTTACCTGCGTGGGGGCGCTTCCGGGGGGGGAAGGCGCTAAAACCGCCCCGGCGCGGCGGCAGTTTTTTGCGGCGCTTCATCAAAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGt < 2:625729/136‑1 (MQ=255) cGGTTCCCGTGGGGGGTGCGCATCATGGGGGGGAGGGCGCTAAAACCGCCCCGGCGCAGCAGCGGCTTTTTGCGGCGTTTCATCTCAACCGGGGGGTTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCa < 2:446091/139‑1 (MQ=255) cGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGa > 1:574489/1‑139 (MQ=255) cgcATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGcgccg > 1:434472/1‑139 (MQ=255) ggtaagggCGCTGAAATCCCCCCGGCGAAGCGGCACTTTTTTGCTGCTTTCCATCTCAAACGGGGTGTTTTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAAcc < 1:636452/133‑1 (MQ=255) cGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAAcc > 1:400182/1‑139 (MQ=255) aaGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACg > 1:504648/1‑139 (MQ=255) ccAGGGCGGCGGAAGCTTTTTGCTGCGCTTAATCAAAAAGGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCgaga < 1:478328/139‑1 (MQ=255) aGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGAGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGAcc > 1:229948/1‑138 (MQ=255) cgcAGCGGCAGCTCTTTGCGGCGCTTCATCACAAACGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTcc < 2:31013/139‑1 (MQ=255) cgcAGCAGCAGCTTTTTGCGGCCCTTCATCTAAAACGGGGGGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTcc < 1:408620/139‑1 (MQ=255) cagcagcTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTc > 2:338498/1‑138 (MQ=255) cctacATCTAAAACGGGGGGGTGAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccgc < 2:400182/135‑1 (MQ=255) cTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccgcc < 1:338498/139‑1 (MQ=255) | GCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCC > NZ_CP009273/4060987‑4061230 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |