Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,344,212 | 0 | A | C | 54.5% | 17.2 / 29.8 | 22 | S725A (TCA→GCA) | arcB | aerobic respiration two‑component sensor histidine kinase ArcB |
Reads supporting (aligned to +/- strand): ref base A (10/0); new base C (0/12); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACA > NZ_CP009273/3344129‑3344297 | ttCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgc > 1:129814/1‑90 (MQ=255) ttCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgc > 2:358545/1‑90 (MQ=255) cccGACTTTTTCTCCCCGGGCCGGAAGGTAAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCcgccgc < 1:342420/89‑1 (MQ=255) ttATCTCCCCGGGCGGGAGGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTa < 1:280503/90‑1 (MQ=255) ttATCTCCCCGGGCGGGAAGGTCAGGAGACTAATTTTGCTGCCCCAGTTGTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTa < 1:383323/90‑1 (MQ=255) tccccAGGCCGAAGGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAAtttt < 1:358545/88‑1 (MQ=255) cccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt > 2:93312/1‑90 (MQ=255) cAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGt > 2:189762/1‑90 (MQ=255) aGGCGGAAAGGTCGGGAAACTAATTTTCCTCCCCCAGTTTTTGCAGAGGGCGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGtc < 2:19406/90‑1 (MQ=255) gAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcct > 2:313148/1‑90 (MQ=255) gAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcct > 2:79167/1‑90 (MQ=255) cgggAGACTGATTTTGCTGCCCCAGTTGTTGCAGGGGGCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa < 1:499730/88‑1 (MQ=255) gagaCTAATTTTCCTGCCCCAGTTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTCTAATTTTATGTCCTTCCTCAACAATGc < 2:129814/90‑1 (MQ=255) agaTTGATTTTCCTGCCCCAGTTTTTGCAGTGGGCCAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc < 2:475761/90‑1 (MQ=255) agaCTGATTTTGCTGCCCCAGTTGTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc < 1:115904/90‑1 (MQ=255) cTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt > 2:459937/1‑90 (MQ=255) tGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATcc > 2:249369/1‑90 (MQ=255) aCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcgc < 1:313148/90‑1 (MQ=255) tGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc > 2:281697/1‑90 (MQ=255) gaaGATGGCGAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg < 2:219210/88‑1 (MQ=255) cGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc < 1:284930/90‑1 (MQ=255) aCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACa > 1:492408/1‑90 (MQ=255) | TTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACA > NZ_CP009273/3344129‑3344297 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |