Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,737,017 | 0 | A | C | 52.6% | 3.1 / 20.2 | 19 | T309P (ACC→CCC) | yiaK | 3‑dehydro‑L‑gulonate 2‑dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/10); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAA‑CCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGAT > NZ_CP009273/3736942‑3737097 | atcatGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGaa > 2:337021/1‑90 (MQ=255) ttACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAA‑ccgccg > 2:48431/1‑90 (MQ=255) ttACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAA‑ccgccg > 1:509194/1‑90 (MQ=255) aCTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTTGCCGAAAA‑CCGCCGTAACGGc > 1:577684/1‑90 (MQ=255) aCTAGTGCCGAGCGTGCTGACAAAAATCAGCCCATTCGCTTCCCCGGCCATAATTTTACTCCCCTGCTGGCCAAAAA‑CCGCCGTAACGGc < 1:183188/90‑1 (MQ=255) agggCCGGGCGTGCGAAAAAAAATCAGGCCTTCCCCTCCCCCGCCCATAATTTTCCTCCCCGGCGGGCAAAAAA‑CCGCCGTAACGGCATc < 2:509194/87‑1 (MQ=255) tGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAA‑CCGCCGTAACGGCATCAc > 1:387756/1‑90 (MQ=255) tGCCGAGCGTGCGAACAAAAACCAGCCCTTCCGCTCCCCCGGCCATGATTTTACTCCCCTGCGGGCCAAAAA‑CCGCCGTAACGGCATCAc < 2:220712/90‑1 (MQ=255) ggCCGAGCGTGCTGAAAAAAATCAGCCCATTCGTTTCCCCGCCAATGAATTTACTCCCCTGCTGGCCAAAAA‑CCGCCGTAACGGCATCAc < 1:238447/89‑1 (MQ=255) cggggCGTGCTAAAAAAAATCACCCCTTTCCCTCCCCCGCCCATAATTTTACTCCCCTGC‑GGCCAAAAACCCGCCGTAACGGCATCACTg < 2:507465/86‑1 (MQ=255) aGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAA‑CCGCCGTAACGGCATCACTGTtg > 1:542693/1‑90 (MQ=255) gCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAA‑CCGCCGTAACGGCATCACTGTTGATGAc > 1:214936/1‑90 (MQ=255) gCTGAAAAAAACCAGGCCATTCCCTTCCCCGGCCATGATTTTACTCCCCTGCTGGCCAAAAA‑CCGCCGTAACGGCATCACTGTTGATGAc < 1:7133/90‑1 (MQ=255) aaaCCAGCCCTTCCGCTCCCCCGCCCATAATTTTCCTCCCCTGC‑GGCCAAAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGtgtg < 2:511130/90‑1 (MQ=255) aGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAA‑CCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCa > 1:576752/1‑90 (MQ=255) ccTTCCCCGCCCATAATTTTCCTCCCCGGCGGGCAAAAAA‑CCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAgg < 1:48431/89‑1 (MQ=255) gCCATGAATTTACTACCCTGCTGGCCGAAAA‑CCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATgag > 1:134942/1‑90 (MQ=255) gCCATAATTTTCCTCCCCGGCGGGCAAAAAC‑CCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATgag < 2:576752/90‑1 (MQ=255) atttttCCTCCCCTGGGGGCAAAAAC‑CCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGAt < 2:454081/87‑1 (MQ=255) | ATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAA‑CCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGAT > NZ_CP009273/3736942‑3737097 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |