| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB122 | 2,780,706 | 0 | A | T | 42.3% | 16.0 / 11.5 | 26 | intergenic (+96/‑29) | ahpC_2/ahpF | Alkyl hydroperoxide reductase C/Alkyl hydroperoxide reductase subunit F |
| Reads supporting (aligned to +/- strand): ref base A (6/9); new base T (10/1); total (16/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.43e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGC > ACB122/2780586‑2780824 | gCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGCCCTAGCC‑TGGGCGTTTTTGTTTCCGAGTTTTGaaa > 6:111463/1‑135 (MQ=255)gCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAAAGCCCGGGCGACCTCGGC‑TGGGCGTTTTTGTTTGCCAGTTTTTaaa > 5:137094/1‑135 (MQ=255) gcaCCTCGTCGGAAATATCTATGGTCTTCCCAATACCGGGGGGGTTAACCACCCGCAAAAGCTTAATAACCTACCCCCCCGTAAAACGCCCGGGCGACCTCGCC‑TGGGGGTTTTTGTATCCGAGTTTTGAAAAAgg < 4:329380/134‑1 (MQ=255) gACCTCGTAGCCAAGAGCTAAGGTCATGCCAATAACGCGGGCGGTAAGCAGCCGCCAAAGCTTAAGTAACTCCCGCCCCGTAAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTTTATCCGAGTTTTGAAAAAgg < 6:206685/135‑1 (MQ=255) aCCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGaa > 6:82719/1‑136 (MQ=255) cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGCCCTCGCC‑TGGGCGGTTTTTTTTACGAGGTTTTAAAAAAGAATAGCCCGGAAGTTTGGCGCCAccggt > 7:29531/1‑131 (MQ=255) cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑GGGGCGGGTTTGGATCCGAGGTTTGAAAAAGGAATCGCCCGTATGTTGGTCGCCACGcgt > 4:242203/1‑133 (MQ=255) cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCACGCC‑TGGGCGTTTTTGTTTCCGACGTTTGAAAAACGAATAGCCCGTATGTTTGAAGCCACGcgt > 4:114401/1‑133 (MQ=255) cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCCGGCGACCACGCC‑TGGGGGGTTTTTTTTTCGAGTTTTTTAAAAAGTAAGGCCCGGATTTTGGGCGCCgcgcgg > 3:15052/1‑133 (MQ=255) ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTGTTTCCGAGTTTTTAAAAAGGAATCGCCCGGTAGTTTGACGCCACGcgtt > 7:319751/1‑132 (MQ=255) ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGGTTTTGTATCCGGGGTTTGGAAAAGGTATAGCCCGTATGTTGGTCGCCCCGCGTa > 3:105625/1‑135 (MQ=255) ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCCGGCGCCCACGCCGGGGGGGTTTTTTTTTCCTGGTTTGTAAAAATGAAAGGCCCG‑GTGTTGGGCGccccccgg > 1:122155/1‑128 (MQ=255) ccATAACGCGGGCGGTAAGAAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCGGGCGACCTCGCC‑AGGGCGTTTTTGTTTCCGAGTTTTGAAAAATGTATCGCGCGTATGTTGGAAGCCTCGCGTa > 8:339427/1‑135 (MQ=255) cggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGaa < 6:143621/136‑1 (MQ=255) gccgccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACtt < 5:82719/136‑1 (MQ=255) gccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGa < 6:322770/136‑1 (MQ=255) gccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTGGATCCGCGGTTTTAAAAACGCATCGCCCGGATGTTGGACGCCCCCCGGAAAACGCCAATGCAAACCGACTTg > 7:54459/1‑134 (MQ=255) gccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑GGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGa < 8:14509/136‑1 (MQ=255) gCCCGGTCAAAGGCGCGGGCGACCTCTCG‑TGGGCGTTTTTGAATCCGAGTTTTGAAAAAGGAATCTCCCGTATGTTGGACTCCAAGCTTTAATCGCTACTGAAAAGCTACTTGGAGCGGGTCACGCAGCCGATCga < 8:10255/136‑1 (MQ=255) gCCCCGTCAAACCCCCGGGCGACCTCGCC‑TGGGCGTTTTTGTTTACGAGTTTTGAAAAAAGAATCGCCCCTAAGTTGGACGCCCCGGGTAAATAGGAACGGAAAACCTAATTGGAGCGGGTCACCCCGCCGATCg > 2:52281/1‑135 (MQ=255) cAACCGCCCGGGCGCCCTCGCC‑TGGGCGTTTTTGTATCCGCGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAACCCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGtt < 4:211046/136‑1 (MQ=255) cGACCTCGCC‑TGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTTTGTTGGACGCCACGCTGAAATCGCCACTGACAACCTACTTGGAGCGGGTCACCCAGCCGCTCGAGATCGTTGCTTGCCTcgcc > 5:252108/1‑134 (MQ=255) cgcg‑TGGGGGTTTTTGTTTCCTATTTTTGAAAAAGTAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGgc < 3:114401/136‑1 (MQ=255) cGCC‑TGGGCGTTTTTGTTTGCGAGTTTTGAAAAAGGAATCGCCCGGTTGTTGGCCGCCCCGCGTAAATCGCAGCTGCAAACCTACTTGGCGCGGGTCGCCCCGGCGCTCGCGCTCGGTGCTGCCCGCGGCGCCgg > 2:317249/1‑135 (MQ=255) tGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCAcc < 7:219361/86‑1 (MQ=255) tGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCAcc > 8:219361/1‑86 (MQ=255) | GCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGC > ACB122/2780586‑2780824 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |