Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,775,868 | 0 | A | C | 53.3% | 22.6 / 39.9 | 30 | V73G (GTA→GGA) | ypjA | adhesin‑like autotransporter YpjA/EhaD |
Reads supporting (aligned to +/- strand): ref base A (10/4); new base C (0/16); total (10/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCT > NZ_CP009273/2775779‑2775954 | tAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATa > 2:1234689/1‑90 (MQ=255) tattaCCCTGCCCCAAAAAATTTTTATGCCTCCCTATATCTTAATACCCCCCAGATTAAAAAAGACTTCCATTAGATCCCCCGCCATGaa < 1:1177154/90‑1 (MQ=255) taGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACAt > 2:443507/1‑90 (MQ=255) tcccccAAAAAATTTTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGTTTTCAATTAGATCCCCCGCCATGAACATTctgc < 1:15002/88‑1 (MQ=255) gCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgct > 2:42245/1‑90 (MQ=255) gCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgct > 2:1177154/1‑90 (MQ=255) acaAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCa > 2:777677/1‑90 (MQ=255) aaTTGTTTTGCCTTCCTATATCTTAATACCCACCAAATTAAATAAGACTTCAATTAGACCCCCCGCCATGAACATTCTGCTGGCCATGGt < 1:39138/90‑1 (MQ=255) gTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGat < 2:820228/90‑1 (MQ=255) gccTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat < 2:864149/88‑1 (MQ=255) aCTTCCTATATTTTGATACCCCCCAGATTAAATAAGTTTTCAATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATg < 1:1107614/90‑1 (MQ=255) cTTCCTATATCTTGATATCCACCAGATTAAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGa < 2:714414/90‑1 (MQ=255) ctCTTGATACCCCCCAGATTAAATAAGTTTTCAATTATATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAtt < 1:927276/89‑1 (MQ=255) atCTTGATATCCACCAGATTCAATAAGACTTCCAATAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAtt < 2:177738/90‑1 (MQ=255) gATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATg > 1:275721/1‑64 (MQ=255) gATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATg < 2:275721/64‑1 (MQ=255) gATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATttgttg > 1:944689/1‑90 (MQ=255) gATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATttgttg > 2:305456/1‑90 (MQ=255) tCCACCAAATTAAATAATTCTTCAATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTAcc < 1:682276/90‑1 (MQ=255) gATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCa > 2:1024814/1‑90 (MQ=255) ttCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgc < 2:1066118/43‑1 (MQ=255) ttCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgc > 1:1066118/1‑43 (MQ=255) tCAATAAGACTTCTATTATATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATcc < 1:777677/90‑1 (MQ=255) ttCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGt < 2:1087551/90‑1 (MQ=255) ttCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGt < 1:1024814/90‑1 (MQ=255) tCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGtt > 2:718479/1‑90 (MQ=255) ccATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTg < 2:621871/90‑1 (MQ=255) caaTTATATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTg < 2:889608/88‑1 (MQ=255) aTTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGcc < 1:396426/90‑1 (MQ=255) ttAGACACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCa < 1:533915/90‑1 (MQ=255) gATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCt < 2:534269/90‑1 (MQ=255) | TAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCT > NZ_CP009273/2775779‑2775954 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |