Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 228,181 | 0 | C | T | 35.9% | 52.9 / 10.1 | 40 | F191F (TTC→TTT) | yafD | endonuclease/exonuclease/phosphatase family protein |
Reads supporting (aligned to +/- strand): ref base C (9/16); new base T (14/0); total (24/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.56e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGT > NZ_CP009273/228044‑228316 | gACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGCTGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGCGGTTTc > 2:576306/1‑138 (MQ=255) acacCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCa < 1:853311/138‑1 (MQ=255) ccGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAGTTCAATggc > 1:342305/1‑137 (MQ=255) tggtggTTAATATACACGCCGTCAACTGCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTACAATGCCTGGAGCCGTA‑g < 1:28215/139‑1 (MQ=255) ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGTCTgg > 2:13364/1‑130 (MQ=255) ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTgg < 1:13364/130‑1 (MQ=255) tggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAAAGCCCGGGGCCGGG‑GAAg > 1:185724/1‑139 (MQ=255) ttAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGAt > 1:337027/1‑139 (MQ=255) aaTATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATTGCTGGAg > 1:467811/1‑125 (MQ=255) aaTATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAg < 2:467811/125‑1 (MQ=255) atatACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGaa < 2:312478/139‑1 (MQ=255) atatACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGa < 1:468670/138‑1 (MQ=255) acGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATGACCGGGGCCGTAGGAAGGAGGaagcggt > 2:753266/1‑134 (MQ=255) acGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTCAATTCCTTGAGCCGTA‑GAAGGATGAAcgcggt > 1:27781/1‑137 (MQ=255) acGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGTTTTAATTCCTGGAGCCGTA‑GAAGGGTGAACGCGtt > 2:442248/1‑139 (MQ=255) gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGTTTTAATTCCTGGAGGCGGA‑GGAGGATGAACGCGGTATATc > 2:135244/1‑139 (MQ=255) gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGTTTTAATGCCTGGAGCCGGA‑GGAGGATGGACGCGGTttatt > 2:16042/1‑138 (MQ=255) tCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCCATGCCTGGGGCCGCA‑GAAGGGTGAAAGCGGTATATCg > 1:330336/1‑139 (MQ=255) aaCTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCt < 1:16042/139‑1 (MQ=255) tCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGGAATGTCGCTgcg > 1:849063/1‑139 (MQ=255) taGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAg < 2:217409/139‑1 (MQ=255) taGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAg < 1:442248/139‑1 (MQ=255) aGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGAGCCGGA‑GAAGGGTGAACGCGTTATATCGCTTTGTGCGGGGAATGTCGCGGCGCCAGGTGCg > 2:564774/1‑139 (MQ=255) aGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTTCGCGGGGAATGTCGCTGCGCCAGGTGCg > 1:192395/1‑139 (MQ=255) aGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATAATGCTTTGCGCGGGGAAGGTCGCTGCGCCCGGGGCg > 2:400450/1‑139 (MQ=255) aGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTATATGCCTGGGGACGCA‑GAAGGAGGGACGAGGTATATTGCTTTTCGCGGGGAAACTGGCCGCGGCCGGTGCg > 2:617666/1‑139 (MQ=255) aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGTTTTAATGCCTGGAGCCGTT‑GAAGGAAGAACGCGTTATATTCCTTTTCGCGGGGAAAGTCGGTGCGCCAGGGGGGGTTTAcc > 1:503644/1‑139 (MQ=255) cTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCg < 1:509879/139‑1 (MQ=255) gCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAgcg < 1:412369/139‑1 (MQ=255) aTCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAAAGCCCGGGGCCGGA‑GAAAGATGACCGCGCTTTATTGCTTTTCGCGGGGAATGTCGCTGCGCCCGGGGCGGTTTTTCGATTATACGCGCCg > 2:714127/1‑139 (MQ=255) aTCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTTAAAGGCCGGAGGCGGA‑GAAGGATGAAAAGGGTATATTGCCTTGTGCGGGGAATGTCGCGGCGCCAGGGGCGTTTTAACGATGATTAGCGCCg > 1:226140/1‑139 (MQ=255) cAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCCGGAGCCGTA‑GAAAGATGGACGCGTTATATCGCTTTGCGCGGGGAAAGTCGCTGCGCCAGGGGCGTTTTATCCATTATAAGCGCCggc > 2:803874/1‑137 (MQ=255) tAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGc < 1:803874/139‑1 (MQ=255) acaGCGGCCCGGTCATTATGGCGGGAGGTTTTAAAGCCTGGAGCCGTA‑GAAGGAGGAAAGCGTTATATTGCTTTGCGCGGGGAAAGTCGCTGCGCCAGGTGGGGTTTTCCGGTGATCAGCGCCGTCGGGCGTTTTGTCg > 1:471490/1‑139 (MQ=255) acaGCGGCCCGGTCATTATGGCGGGAGATTTCAAAGCCTGGAGCCGTA‑GAAGGATGAACGCGGTATATCGCTTTGCGCGGGGAATTTCGCCGCGCCAGGGGCGTTTTACCGGTGATCAGCGCCGTCGGGCGTTTGGTCg > 1:943111/1‑139 (MQ=255) gCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCt < 2:185724/139‑1 (MQ=255) gCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCt < 1:853639/139‑1 (MQ=255) cGGTCATTATGGCGGGAGACTTCAATGCCTTGGGCCGGG‑GGAGGATGAACGCCGTGTATAGCTTCTCGCCGGAAAATGCGCCGCGCCCGGGCGGTTTTTCCGGTGATCAGCGCCGTTAGGCGTTTGGTGGCCCGCTCGa > 2:858460/1‑139 (MQ=255) gTCATTATGGCGGGAGATTTTAATGCCTGGAGGCGTA‑GAAGGATGGAAGCGGTATATTGCTTTGCGCGGGGAATGTCGCTGCGCCAGGGGCGGTTTACCGGTGATCAGCGCCGTCGGGGGGTTGGGGGCCCGCCCGGtt > 1:163345/1‑139 (MQ=255) gagaTTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCg < 1:135244/139‑1 (MQ=255) tttCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGt < 2:909603/139‑1 (MQ=255) | GACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTA‑GAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGT > NZ_CP009273/228044‑228316 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |