Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,291,076 | 0 | A | C | 63.0% | 8.2 / 28.1 | 27 | V141G (GTG→GGG) | yraQ | permease |
Reads supporting (aligned to +/- strand): ref base A (10/0); new base C (0/17); total (10/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-06 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTCCCAGCAGCGTG > NZ_CP009273/3290988‑3291165 | ccAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGacac > 2:747845/1‑90 (MQ=255) aaGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGc > 1:363656/1‑90 (MQ=255) aaGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGc > 2:1000384/1‑90 (MQ=255) cATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTcc > 2:442982/1‑90 (MQ=255) acCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCcgcg > 2:527150/1‑90 (MQ=255) aGCGCCGCCGGGTTTAACCCCGGTTTCCCCACCCAGAACCCCGGCCCCCCCCCCACCGCCCCCTGTGGGCGCCGCTTCCCCGCCGCGAcc < 1:983545/90‑1 (MQ=255) gCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGcc > 1:353270/1‑90 (MQ=255) ttttACCCCCGGATTGCCCACCAAGAATGCCAGCCCCCCCCCCACCGACCCCTTTTGGCGACGCATTCCCGCCGCGACCGGAGCCGcaca < 1:424551/89‑1 (MQ=255) gTTTAACCCCGGATTGCCCACCCAAAATGCCAGCGCCCCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGcaca < 2:682135/90‑1 (MQ=255) ttAACCCCGGATTCCCCACCCAAAATGCCAGCGCCCCCCCCACCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagc < 1:339780/90‑1 (MQ=255) taaaCCCCGGATTGCCCACCCAGAATGCCAGCGCCCCCCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagc < 2:829678/88‑1 (MQ=255) taccacGGGTTTCCCCATCCCGAAGCCCCCCCCCCCCCCCACCGACCCCTTTTGGCGCCGCCTTCCCGCCGCGACCGGAGCCGCAcagca < 1:897323/87‑1 (MQ=255) tAACCCCGGTTTGCCCACCCAGAACGCCGGCGCCCCCCCCCCCGCCCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagca < 1:440444/90‑1 (MQ=255) tAACCCCGGATTGCCCACCCAGAATGCCAGCGCCCCCCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagca < 2:21064/90‑1 (MQ=255) tAACCCCGGATGCCCCACCCAGAATGCCAGCGCACCGCCCATCGACCCCTGTGGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagca < 2:587807/90‑1 (MQ=255) tAACACCGGATTGCCCATCCAGAATGCCAGCGCCCCCCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagca < 2:67729/90‑1 (MQ=255) tAACACCGGATTGCCCACCCAGAATCCCAGCCCACCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagca < 2:714889/90‑1 (MQ=255) ccGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTac > 1:67729/1‑90 (MQ=255) ggATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTacac > 2:644526/1‑90 (MQ=255) ggATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTacac > 1:1065499/1‑90 (MQ=255) ccACCCAGAATCCCAGCCCCCCCCCCATCGACCCCTGTTGGCGCCGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGc < 2:912868/90‑1 (MQ=255) aGAACGCCAGCCCCCCCCCCACCCACCCCTGTGGGCGCCGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCa < 1:1042428/90‑1 (MQ=255) cgccccgcccATCGCCCCCTGTGGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAg < 2:930258/90‑1 (MQ=255) cgccccgcccACCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAg < 2:511503/90‑1 (MQ=255) ccccccGCCCACCGACCCCTGTGGGCGCCCCATCCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAg < 1:660704/86‑1 (MQ=255) gcACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGc > 1:974671/1‑90 (MQ=255) ccGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGtt > 2:384549/1‑90 (MQ=255) cccATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTccc < 2:418227/90‑1 (MQ=255) acacCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTCCCAGCAGCg > 1:800424/1‑90 (MQ=255) acacCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTCCCAGCAGCg > 1:694290/1‑90 (MQ=255) acCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTCCCAGCAGCGTg > 2:1006522/1‑90 (MQ=255) | CCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTCCCAGCAGCGTG > NZ_CP009273/3290988‑3291165 |
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: ‑ |