Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,806,632 | 0 | C | T | 44.8% | 16.5 / 13.1 | 29 | S181L (TCG→TTG) | coaBC | bifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase CoaBC |
Reads supporting (aligned to +/- strand): ref base C (7/9); new base T (13/0); total (20/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.21e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCG > NZ_CP009273/3806504‑3806760 | agagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCaa < 1:57513/138‑1 (MQ=255) cttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTTCCCCTCCACGACCtt > 2:609208/1‑137 (MQ=255) tGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAAAATCTGAACATTATGATTAccg > 1:537086/1‑139 (MQ=255) ggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCATTTTTTGCCCGCCAACAACCTGAAACATATGAAAATTATGATTACCGCCgg > 2:163481/1‑139 (MQ=255) cAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCGGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGccc < 1:163481/139‑1 (MQ=255) cAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGaa < 1:609208/139‑1 (MQ=255) aGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGTGCCTTTTTTGCCCCGCAAACACCTGAAACATCTGAACATTATGATTTACGCCCCCCCGACGCGGCGACCGCt > 2:412802/1‑139 (MQ=255) gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTCCCCGTCAAAGACCTGAAACATCTGGAAATTATGATTACCGCCGGCCCGCCGCGTGAAACGccc > 2:60853/1‑136 (MQ=255) gtggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATc < 2:537086/138‑1 (MQ=255) ggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCCTCAACGAACTGAAACATCTGAACATTATGGTTACCCCCCGCCCGACGCGTGAACCGCCCGATCCgg > 1:515715/1‑139 (MQ=255) atatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTTCCGCCGGCCCGACGCGGGAAACGCCCGCTCCCGGGCg > 2:302789/1‑139 (MQ=255) atatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCCCAACGACCTTGAACAACTGAACATTATGATTACCGCCGGCCCGACGCGTGGAACGCTCGATCCGGTGCg > 2:14838/1‑139 (MQ=255) tCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTta < 1:302789/139‑1 (MQ=255) tGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCGTTTTTTGCCCCTCCACGACCCTAAACAACTGAACATTATGATTTCCGCCCGCCCGACGCGTGAACCCCCCCATCCCGGGCGGTTTAtctc > 2:233161/1‑138 (MQ=255) tGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGATGCGTGAACCGCTCGATCCGGTGCGTta > 2:156452/1‑132 (MQ=255) tGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGATGCGTGAACCGCTCGATCCGGTGCGTta < 1:156452/132‑1 (MQ=255) ggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACAt < 2:85614/74‑1 (MQ=255) ggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAAAAt > 1:85614/1‑74 (MQ=255) aTGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCCACGACCTGAAAAATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCCCGATCCGGTGCGTTTTATCTCTAATca > 1:423648/1‑137 (MQ=255) aTGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCCTCAACGACCTGAAAAATCTGAACATTATGATTACCGCCGGCCCGACGCGGGAACCGCCCGATCCGGTGCGTTATATCTCTAATca > 1:349353/1‑137 (MQ=255) ttGTGGATATGGCGGTAGCGCCTTTTTTGCCCCTCAACGACCTGGAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCCATCCGGTGCGTTATATCTCTAATCCCAGCTCCGGCAAGATGGGgtt > 2:99016/1‑139 (MQ=255) tgtgGATATGGCGGTAACGCATTTTTTGCCCGTCAACAACCTGAAAAAACTGAACATTATGATTACCGCCGGGCCGACGCGTGAAACCCCCCATCCGGTGGGTTATATCTCTAATCACAGCTCCGGGAAAAAGGGtttt > 1:47835/1‑139 (MQ=255) tgGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCCCCCGATCCCGTGCGTTTTATCTCTAATCACAGCTCCGGCCAGATGGGgttttc > 1:545934/1‑137 (MQ=255) atGGCGGTAGCGCATTTTTTGCCCGGCCACGACCTGAAACATCTGGACATTTTGATTTCCGCCGGCCCGACGCGTGAAACGCCCGATCCGGTGCGGTTTATCTCTAATCACAGCTCCGGCAAGATGGGGTTTGCTATcg > 2:481272/1‑139 (MQ=255) ggCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCCCGATCCGGTGCGTTATATCTCTAATCACAGCCCCGGCAAGATGGGTTTTGCTATcgcc > 2:106935/1‑139 (MQ=255) gCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGaa < 2:247782/75‑1 (MQ=255) gCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGaa > 1:247782/1‑75 (MQ=255) gagagcATTTTTCGGCCGGCAACGAACCGTAAAATCTGAACATTTTGATTAACGCCGGCCCGACACGAGGACCGCACGATCCGGTGCGTTATATCTTTAAACATAAATCCGGCCAGAAGGGGTTTTCTATcccccccgc > 2:597617/5‑139 (MQ=255) gcgcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgccgccg < 1:470707/139‑1 (MQ=255) | AGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCG > NZ_CP009273/3806504‑3806760 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |