Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,344,212 | 0 | A | C | 51.5% | 35.5 / 51.4 | 33 | S725A (TCA→GCA) | arcB | aerobic respiration two‑component sensor histidine kinase ArcB |
Reads supporting (aligned to +/- strand): ref base A (15/1); new base C (0/17); total (15/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.54e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAT > NZ_CP009273/3344124‑3344298 | aTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAc > 2:255533/1‑90 (MQ=255) cccGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTTTTGCAGAGGGCGTACCCCCCCTGCCCCcgccgc < 1:945310/89‑1 (MQ=255) ccGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCa > 2:543678/1‑90 (MQ=255) aCGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCt > 2:619153/1‑90 (MQ=255) tttttctcCCCGGGCGGGAGGGTCGGGAGACTGATTTTGCTGCCCCATCTTTTGCGGAGGGCGAACCCCCCCTGCCCCCGCCGCACCttt < 2:880803/86‑1 (MQ=255) tttATCTCCCCGGGCCGGAAGGTCGGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCttt < 1:44786/89‑1 (MQ=255) cTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAttt > 2:1127789/1‑90 (MQ=255) tCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTa > 1:397653/1‑90 (MQ=255) cccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt > 1:384211/1‑90 (MQ=255) gggCGGGAGGGCCAGGAGACTAATTTTCCGCCCCCATTTTTTGCGGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGtc < 2:859822/89‑1 (MQ=255) gggCCGGAGGGTCAGGAGACTGATTTTGCCGCCCCAGCTGTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGtc < 1:948785/89‑1 (MQ=255) aGGCGGGAGGGTAAGGAGACTGAATTTGCTGCCCCAGTTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtc < 1:823186/90‑1 (MQ=255) gCCGGAGGGTAGGGAGACTAATTTTGCTCCCCCAGTTGTTGAAGAGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcct < 1:770867/90‑1 (MQ=255) tCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca > 2:116409/1‑90 (MQ=255) agaCTGATTTTGCTGCCCCAGTTGTTGCAGGTGGCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc < 1:153620/90‑1 (MQ=255) agaCTGATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc < 1:110122/90‑1 (MQ=255) agaCTAATTTTGCTGACCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc < 1:144199/90‑1 (MQ=255) agaCTAATTTTGCTCCCCCAGTTTTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc < 1:996970/90‑1 (MQ=255) cTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt > 2:358126/1‑90 (MQ=255) cTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt > 1:168825/1‑90 (MQ=255) cTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt > 1:768335/1‑90 (MQ=255) cTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt > 1:226899/1‑90 (MQ=255) gattttGCTGCCCCGGCTTTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt < 1:718561/87‑1 (MQ=255) ccAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg > 1:899543/1‑90 (MQ=255) aGCTGTTGAAGGTGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc < 1:116409/90‑1 (MQ=255) gtGTTGAAGGTGGCGTACCCCCCCCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAg < 1:543678/89‑1 (MQ=255) gaagagGGCGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg < 1:1022239/88‑1 (MQ=255) aTGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc > 1:621657/1‑90 (MQ=255) ccACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc < 1:753619/83‑1 (MQ=255) ccACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc > 2:753619/1‑83 (MQ=255) aCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACa > 2:821083/1‑90 (MQ=255) cTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt < 1:446331/90‑1 (MQ=255) cTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt < 2:621657/90‑1 (MQ=255) | ATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAT > NZ_CP009273/3344124‑3344298 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |