Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,291,063 | 0 | A | C | 55.6% | 16.1 / 20.9 | 18 | G145G (GGT→GGG) | yraQ | permease |
Reads supporting (aligned to +/- strand): ref base A (7/1); new base C (0/10); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCAAAACCCCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAG > NZ_CP009273/3290979‑3291149 | cgcAAAACCCCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGcc < 1:145184/90‑1 (MQ=255) ccAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCACATCGacac > 1:210562/1‑90 (MQ=255) aGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCt > 2:544756/1‑90 (MQ=255) gCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTg > 1:533134/1‑90 (MQ=255) aGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTgg > 1:447142/1‑90 (MQ=255) aGCGTCGCCGGGTTTAACCCCGGATTGCCCATCCAGAATGCCAGCGCCCCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGAcc < 2:202330/90‑1 (MQ=255) aGCGTCGCCGGGTTTAACCCCGGATTGCCCACCCAGAATGCCGGCGCCCCCCCCACCGACCCCTGTTGGCGACGCATTCCCGCCGCGAcc < 1:544756/90‑1 (MQ=255) aGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCCCCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGAcc < 2:72770/90‑1 (MQ=255) aGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCCCCGCCCACCGACCCCTGTTGGCGACGCATTCCCGCCGCGAcc < 2:28713/90‑1 (MQ=255) gCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCg > 1:164879/1‑90 (MQ=255) gCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGcc > 2:318335/1‑90 (MQ=255) gCCGGGTTAAACCCCGGATTGCCCACCCAGAATGCCAGCGCCCCCCCCACCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGcc < 1:159567/90‑1 (MQ=255) ttttAACCCCGGATCCCCCATCCAGAATGCCAGCGCCCCCCCCATCGCCCCCTGTGGGCGACGCATTCCCGCCGCGACCGGAGCCGcaca < 2:538192/89‑1 (MQ=255) gTTTACCCCCGGATTGCCCATCCAGAATGCCGGCGCCCCGCCCATCGCCCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGcaca < 2:301182/90‑1 (MQ=255) ccGGATTCCCCATCCAGAATGCCAGCCCCCCCCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTac < 1:125326/90‑1 (MQ=255) cGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTaca > 1:343279/1‑90 (MQ=255) acccAGAATGCCAGCGCCCCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGccc < 2:365289/88‑1 (MQ=255) cAGCGCCCCGCCCATCGACCCCTTTTGGCGCCGCATCCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGaaaa < 2:373724/90‑1 (MQ=255) cgccccgcccATCGACCCCTTTGGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAg < 2:371981/90‑1 (MQ=255) | CGCAAAACCCCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAG > NZ_CP009273/3290979‑3291149 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |