Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,312,206 | 0 | A | . | 72.4% | 40.6 / 38.7 | 29 | coding (41/399 nt) | yciA | acyl‑CoA thioesterase |
Reads supporting (aligned to +/- strand): ref base A (4/4); new base . (9/12); total (13/16) | |||||||||||
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.00e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGA > NC_000913/1312133‑1312282 | cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTaaacaag < 2:235231/80‑5 (MQ=255) aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAg < 2:235232/81‑1 (MQ=255) aaccaaccACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAggg < 1:235233/80‑1 (MQ=255) accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGa > 2:235255/1‑81 (MQ=255) accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGAc < 2:235234/81‑1 (MQ=255) ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGAc < 2:235235/81‑1 (MQ=255) ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACg < 1:235236/81‑1 (MQ=255) ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACg > 1:235269/1‑81 (MQ=255) aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGtt < 2:235237/81‑1 (MQ=255) ccaccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTa > 1:235270/1‑81 (MQ=255) aTGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa < 2:235238/81‑1 (MQ=255) tCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACa > 1:235260/1‑81 (MQ=255) tCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACa < 1:235239/81‑1 (MQ=255) aCCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACAt > 1:235268/1‑81 (MQ=255) cattggcattggTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATgg > 2:235263/1‑80 (MQ=255) attggcattggTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTa < 2:235240/81‑1 (MQ=255) ggcattggTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaa < 2:235242/80‑1 (MQ=255) ggcattggTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaa < 1:235241/80‑1 (MQ=255) cattggTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAcc > 1:235266/1‑81 (MQ=255) ttggTATCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGa < 1:235243/81‑1 (MQ=255) aTCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGGCGTTATGTGTTTTAGACATGGTAAAACCGACTTa > 2:235265/1‑80 (MQ=255) aTCGGCGGGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTa > 2:235264/1‑80 (MQ=255) ggcggGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAg > 1:235276/1‑81 (MQ=255) ggcggGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAg > 1:235257/1‑81 (MQ=255) cggGCATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGa < 1:235244/81‑1 (MQ=255) ggCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGaa < 2:235245/81‑1 (MQ=255) cATGGCTAAAGTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAAt > 2:235271/1‑81 (MQ=255) gCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATca > 2:235273/1‑80 (MQ=255) gTACGT‑AAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGcgcg < 2:235246/81‑1 (MQ=255) cGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGa < 1:235247/80‑1 (MQ=255) | CATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGA > NC_000913/1312133‑1312282 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |