Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 782,311 | 0 | A | C | 51.3% | 47.1 / 60.0 | 37 | G247G (GGT→GGG) | gpmA | 2,3‑diphosphoglycerate‑dependent phosphoglycerate mutase |
Reads supporting (aligned to +/- strand): ref base A (15/3); new base C (0/19); total (15/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.71e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAA > NZ_CP009273/782222‑782400 | aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTa > 2:418412/1‑90 (MQ=255) aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTa > 1:1587685/1‑90 (MQ=255) aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTa > 2:1896717/1‑90 (MQ=255) aaTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCt < 1:19371/90‑1 (MQ=255) aTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg > 1:1142611/1‑90 (MQ=255) aTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg > 2:617523/1‑90 (MQ=255) tCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGc > 2:923041/1‑90 (MQ=255) tCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGc > 2:533082/1‑90 (MQ=255) gTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgc > 2:1350369/1‑90 (MQ=255) gTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgc > 2:1596502/1‑90 (MQ=255) gTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgc > 1:705550/1‑90 (MQ=255) tAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgccg > 1:1520909/1‑90 (MQ=255) aaaaaGCCGCCTCCCTTGCAGTCGTTTTTCTCATTTTAAACGAATGCCTTTTCCTTCCCTTTCCCTTGGTTTGCAACCGCCGCTGCTTTc < 1:617523/90‑1 (MQ=255) aaaGCCGCCTCACTTGCAGTCGGTTTTCTCTTTTAAAACGAATGACGTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGc < 1:1619877/90‑1 (MQ=255) aaaGCCGACTCCCTTGCAGTCGGCTTTCTCATTTTAAACAAATGCCGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGc < 2:2008372/90‑1 (MQ=255) ctggCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc < 2:28868/87‑1 (MQ=255) cTTGCGGTCGGTTTTCTCTTTTAAATCGAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc < 1:411335/90‑1 (MQ=255) cTTGCAGTCGGTTTTCTCTTTTAAAACGAATGCCTTTTACTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc < 1:923041/90‑1 (MQ=255) cTTCCAGTCGGCTTTCTCTTTTTAAACGAATGCCTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc < 1:1896717/90‑1 (MQ=255) ggTCGGTTTTCTCTTTTTAAACAAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt < 1:533082/89‑1 (MQ=255) aGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTCCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt < 2:1946082/90‑1 (MQ=255) aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt < 2:734734/90‑1 (MQ=255) tCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCCCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTa < 2:1084672/90‑1 (MQ=255) gtcTTTCTCATTTAAAACAAATGACTTTTCCTTCCTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc < 2:1222350/88‑1 (MQ=255) ggttttCTCATTTTAAACGAATGCCTTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc < 2:1587685/87‑1 (MQ=255) ggttttCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc < 2:208209/87‑1 (MQ=255) ggCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc < 2:1942165/90‑1 (MQ=255) gttttCTCTTTTTAAACGAATGACTTTTTCTTCCCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc < 1:1308949/88‑1 (MQ=255) gttttCTCTTTTCAAACAATTGTTTTTTCCTTCGCTTTCCCTTGTTTTGAACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc < 2:817558/88‑1 (MQ=255) gCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc > 2:923641/1‑90 (MQ=255) ttCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGa > 1:300798/1‑90 (MQ=255) aaTGACTTTTCCTTCGCTTTCCCCTGGTTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc < 2:1835825/90‑1 (MQ=255) aTGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg > 1:1347727/1‑90 (MQ=255) acttttACTTCCTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt < 1:841046/87‑1 (MQ=255) ttACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTg < 2:1489958/44‑1 (MQ=255) ttACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTg > 1:1489958/1‑44 (MQ=255) ttCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtc > 2:617832/1‑90 (MQ=255) tCGCTTTACCCTGGTTTGCAACCGCCGCTGCCTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct > 2:710670/1‑90 (MQ=255) tttCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCCGCGGTTTGAAATTCtcgtc < 2:891151/90‑1 (MQ=255) ttACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg > 1:1070514/1‑90 (MQ=255) aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCTGCGGTTTGAAATTCTCGTCGaa > 2:1103121/1‑90 (MQ=255) | AGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAA > NZ_CP009273/782222‑782400 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |