Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,520,988 | 0 | T | G | 57.1% | 5.1 / 10.3 | 14 | T161P (ACT→CCT) | sgcX | M42 family metallopeptidase |
Reads supporting (aligned to +/- strand): ref base T (3/3); new base G (8/0); total (11/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.49e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAG > NZ_CP009273/4520902‑4521074 | gAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGGGAc > 2:283845/1‑90 (MQ=255) gTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATAAAGAGGGACAAGCGGAccc > 1:75368/1‑89 (MQ=255) gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGGGACAGGCGGACCCACCTg > 2:228379/1‑90 (MQ=255) gCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGGGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtcttt > 2:202302/1‑90 (MQ=255) gCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGGGACAGGCGTACCCACCCGAATGCCCCTCCGGATAGGGGCgtcttt > 1:110724/1‑90 (MQ=255) tCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGGGACAGGGGGACCGACCCGAATGCCCATCCGGATAGCGTCGTCTTTGTCTTTGGcg > 1:195609/1‑89 (MQ=255) gTTTGCCAGGAGTTGCGGCGGGGTATACAGAGGGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTTTCTTTTGCAcc > 1:21945/1‑90 (MQ=255) tGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCCTTGTCTTTGGCACCAat > 2:249832/1‑90 (MQ=255) cAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATc < 1:136260/90‑1 (MQ=255) gTTGCGGCGGGGTATAAAGAGGGACAGGCGGACCGACCCGAATGCCCATCCGTATAGCGGCGTCTTTGTCTTTTGCACCAATATCAATcc > 1:247317/1‑90 (MQ=255) agagTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGaa < 1:283845/90‑1 (MQ=255) gagTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCa < 1:195461/62‑1 (MQ=255) gagTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCa > 2:195461/1‑62 (MQ=255) gagTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAg > 1:231081/1‑90 (MQ=255) | GAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAG > NZ_CP009273/4520902‑4521074 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |