Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,454,362 | 0 | T | G | 56.2% | 10.9 / 16.0 | 16 | V351G (GTG→GGG) | paaC | 3‑hydroxyacyl‑CoA dehydrogenase PaaC |
Reads supporting (aligned to +/- strand): ref base T (1/6); new base G (9/0); total (10/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.74e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAC > NZ_CP009273/1454278‑1454443 | ttgaGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGGAGTGATCGATAAAATGGCGGGGAAAggggtgg > 1:311402/1‑89 (MQ=255) caAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGGGACCATtgct > 2:335925/1‑90 (MQ=255) gCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGGGGGGACCATtgctgctg > 2:549197/1‑90 (MQ=255) gagaCGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGGGGTGACCATTGCTGCTGCa > 1:131328/1‑90 (MQ=255) gCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTg < 2:186476/90‑1 (MQ=255) gCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGGGGTGACCATTGCTGCTGCAGCGGTg > 2:307690/1‑90 (MQ=255) tGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGGGGGGACCATTGCTGCTGCAGCGGGGAACCCGGACt > 2:548883/1‑90 (MQ=255) cACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAgg < 1:270778/90‑1 (MQ=255) gtAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTAc < 1:380832/90‑1 (MQ=255) gtAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTAc < 1:429219/90‑1 (MQ=255) tCGATAAAATGGCGGGGAAGGGGGGGACCATTGCTGCTGCAGCGGGGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAAc > 2:461650/1‑90 (MQ=255) aTAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAAcagc < 1:548883/90‑1 (MQ=255) aaaaTGGCGGGCAAGGGGGGGACCATTGCTGCTGCAGCGGTGAAACCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCCAcagccg > 1:285812/1‑88 (MQ=255) gggCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAc < 2:73821/90‑1 (MQ=255) gggCAAGGGGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGGCCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAc > 1:347704/1‑90 (MQ=255) gggCAAGGGGGGGGCCATTGCTGCTGCAGCGGGGAACCCGGACTCAGGGACCCGCCAGGCCATTTATTACCTGCAACAGCAGGGCAAAAc > 1:441669/1‑90 (MQ=255) gggCAAGGGGGGGACCATTGCTGCTGCAGCGGGGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAc > 1:110673/1‑90 (MQ=255) | TTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAC > NZ_CP009273/1454278‑1454443 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |