Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 69,900 | 0 | A | C | 56.2% | 21.0 / 33.4 | 32 | V370G (GTG→GGG) | thiP | thiamine/thiamine pyrophosphate ABC transporter permease ThiP |
Reads supporting (aligned to +/- strand): ref base A (8/6); new base C (0/18); total (8/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.86e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCAC > NZ_CP009273/69821‑69977 | gttgttGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGc > 1:441899/1‑90 (MQ=255) tgAGCAGAAAAAAGAAGCCGGTAGCCAGCACAACCCCCGGCAGGGCGAGAATCAACAGGCCGCTCATCTCCAGCACCGGACCCGCCAGCa < 2:513148/90‑1 (MQ=255) aaaaaaAGCCGGTACCCAGCAAAACCCCCGGCATGGCGAGGATCAACATGCCGCCCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGcc < 2:806166/90‑1 (MQ=255) aaaaGAAGCCGGTAGCCGGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGcc < 2:227054/90‑1 (MQ=255) aaaaGAACCGGGTACCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGcc < 1:368186/90‑1 (MQ=255) aaaaGAACCCGGTAGCCACCACAACCCCCGGCAGGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGcc < 1:265680/90‑1 (MQ=255) gCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCg > 1:285760/1‑90 (MQ=255) gCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCg > 1:664788/1‑90 (MQ=255) ccGGTGCCCCGCCCAACCCCCGGCGGGGGGGGGATCAACAGGCCGCTCATCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGc < 1:5327/90‑1 (MQ=255) ccGGTAGCCACCACAATCCCCGGCATGGCGAGGATCAACATGCCCCTCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGc < 2:587902/90‑1 (MQ=255) gTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGt > 1:95638/1‑90 (MQ=255) gCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGAc > 1:153409/1‑90 (MQ=255) cccaACCCCCGCCAGGGCGAGGACCAACAGGCCCCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 1:362579/88‑1 (MQ=255) cacaACCCCCGGCAGGGCGAGGATCAACATGCCCCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 2:664788/90‑1 (MQ=255) cacaACCCCCGGCAGGGCGAGGACCAACATGCCGCTCATCTCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 1:279919/90‑1 (MQ=255) caccccccGCCAGGGGGAGCATCAACAGGCCCCCCCTCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGc < 1:714780/86‑1 (MQ=255) caATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGc < 2:225349/90‑1 (MQ=255) aTCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTTCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTc > 2:279919/1‑90 (MQ=255) cGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATag > 1:455199/1‑90 (MQ=255) ggCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 2:153409/90‑1 (MQ=255) ggCATGGCGAGGATCAACATGCCCCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 2:441899/90‑1 (MQ=255) ctgggCGAGAAACAACATCCCCCTCACCTCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagcag < 2:631464/87‑1 (MQ=255) cagggCGAGGACCAACATGCCCCTCACCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagcag < 1:320581/87‑1 (MQ=255) ggagaATCAACAGCCCCCTCATCCCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 2:686822/89‑1 (MQ=255) cGAGGATCAACACGCCGCTCATCTCCACCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 2:399552/90‑1 (MQ=255) ggATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATgg > 2:287199/1‑90 (MQ=255) aTCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc < 2:53890/90‑1 (MQ=255) cATGCCGCTCATCCCCACCACCTAACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCAc < 1:110922/90‑1 (MQ=255) ccctcctcCCCCGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 1:400721/88‑1 (MQ=255) cGCTCATCTCCACCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:377030/90‑1 (MQ=255) cGCTCATCCCCAGCACCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:837457/90‑1 (MQ=255) cTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGcac < 2:95638/90‑1 (MQ=255) | GTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCAC > NZ_CP009273/69821‑69977 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |