Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,738,361 | 0 | T | G | 57.9% | 0.4 / 10.2 | 19 | G66G (GGT→GGG) | yiaN | 2,3‑diketo‑L‑gulonate transporter large permease YiaN |
Reads supporting (aligned to +/- strand): ref base T (0/8); new base G (10/1); total (10/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC > NZ_CP009273/3738283‑3738440 | ccAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATAATGaa > 2:37910/1‑90 (MQ=255) ccAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGaa > 2:4544/1‑90 (MQ=255) tCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATgcgg > 2:99486/1‑90 (MQ=255) tGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATgcggggg > 2:213544/1‑88 (MQ=255) aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGCGGGGGGCCTGtt > 2:2831/1‑89 (MQ=255) aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGCGGGGGGCCTGTc > 2:108595/1‑90 (MQ=255) tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGCGGCCTGTGAAAAc > 1:200890/1‑90 (MQ=255) gTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGc < 1:99486/90‑1 (MQ=255) gTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGATTGTCAAAACGc < 1:213544/90‑1 (MQ=255) gAGCCGATAGCTTCTCCCTGCTTGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACATGAATGCGGGGGGCTTTTCAAAAAGCAttgttt > 1:260505/1‑89 (MQ=255) gCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACATGAAAGCGGGGGGCCTGTCAAAAAGCATTTTTGAc > 2:140525/1‑90 (MQ=255) ccGATAGCTTCCCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATAATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGccc < 2:272848/90‑1 (MQ=255) ctgGCGATTCCTTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACtggtg < 1:214801/90‑1 (MQ=255) ctgGCGATTCCGTTCTTTGTGCGGGGGGGGAAAATCAGAAAGGGGGGGGGCTTGACAAAACGCTTTGTTGCCCTGCCGATGAAACtggtg < 2:202466/90‑1 (MQ=255) ttgCGATTCCGTTTTTTGTGCGGGGGGGTAAAATCATGAAGGCGGGCGGTTTGTAAAAACGCTTTGTTGCCCTGCCGATGAAACtggtgg < 2:260505/88‑1 (MQ=255) ccGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATaaa < 1:108595/90‑1 (MQ=255) ttttGTGCGGGCGGGTAAAACAATGAAGGGGGGGGGTTTGTAAAACCGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAgg < 2:200890/89‑1 (MQ=255) ttGTGCTGGCGGGGGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGGCCTGCCGATGAAACTGGGGGGGGATAAAACAGGCg > 1:75989/1‑90 (MQ=255) tgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc < 1:37910/90‑1 (MQ=255) | CCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC > NZ_CP009273/3738283‑3738440 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |