Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP002094 | 500,260 | 0 | G | C | 38.9% | 79.7 / 27.3 | 54 | Q121E (CAA→GAA) | LLNZ_02620 | peptide binding protein |
Reads supporting (aligned to +/- strand): ref base G (25/8); new base C (0/21); total (25/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.12e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCAACTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCATCTTT > CP002094/500162‑500353 | aTCAACTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTg > 1:475296/1‑99 (MQ=255) tCAATTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgc < 1:1140867/99‑1 (MQ=255) tCAACTTTTTTTTTCCCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgc < 1:759513/99‑1 (MQ=255) tCAACTTTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgc < 2:598302/99‑1 (MQ=255) tCAACTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgc > 1:1376705/1‑99 (MQ=255) tCAACTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgc < 2:1395829/99‑1 (MQ=255) tCAACTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgc > 1:3964/1‑99 (MQ=255) aCTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgccag > 2:560170/1‑99 (MQ=255) aCTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgccag < 1:1069928/99‑1 (MQ=255) ctttttttttCCCGGCCATTAATT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgccagc < 1:1616963/96‑1 (MQ=255) cTGTTTTTTTCCCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgccagc < 1:1500250/99‑1 (MQ=255) cTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGTGGT‑CAACGGTTCTTTgccagc > 2:300125/1‑99 (MQ=255) cTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTgccagc < 1:294789/99‑1 (MQ=255) ggCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATc > 2:1493406/1‑99 (MQ=255) ccATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTCCCAGCCATAAACAAAATCAg < 2:789691/99‑1 (MQ=255) ccATTATTT‑CGCTGGAATTAGCAACTGGGGTAAAGGAATAAGCATATTGTAACCCTGTTTGGGGGT‑CAACGTTTCTTTCCCAGCCATAAACAAAATCAg < 1:596834/99‑1 (MQ=255) ccATTATTT‑CGCTGGAATTAGCAACGGGTGTCAGGGCATAAGCATATTGGGAACCAGTTTTGGGGT‑CAACGGTTCTTGGCCAGCCATAAACAAAATCAg < 2:972217/99‑1 (MQ=255) ccttattTCCGCGGGAATTAGCAACTGGTGAGAAGCCATACGCATTTGGTGACCCAGTTTGGGGGC‑CACCGGTTTTTTCCCAGCCATAAACAAAATCAg < 1:560170/97‑1 (MQ=255) ttattT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAATGGTTCTTTGCCAGCCATAAACAAAATCAGCgg > 1:33322/1‑99 (MQ=255) tattT‑CGCTGGCATTAGCAACTGGTGTAAAGGCATAAGCATTTTGTGAACCAGTTTTGGGGT‑CACCGGTTTTTTCCCAGCCATAAACAAAATCAGCGGc < 2:1212570/99‑1 (MQ=255) ttt‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCg > 2:142144/1‑99 (MQ=255) cGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTa < 2:1328829/99‑1 (MQ=255) cggacATAAGCACGTGGTGTCAAGACATAAGCATATTGTGACCCATTTTGGGGGT‑CAACGGTTTTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAAc < 1:1400519/95‑1 (MQ=255) cTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAAc > 1:776604/1‑99 (MQ=255) ggCATAAGAAACGGGTGTAAAGAAATAAGCATTTTGTGAACCATTTTGGGGGC‑CACCGTTTCTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACtt < 1:872977/99‑1 (MQ=255) gCATTAGCAACTGGTGTCAAGGCATAAGCATTTTGTGAACCATTTTTGGGGT‑AAAGGTTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTa < 2:812126/99‑1 (MQ=255) cATTAGAAACGGGTGTAAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CACCGGTTCTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAg < 2:669049/99‑1 (MQ=255) aTTAGCAACTGGTATAAAGGAATAAGCATATGGTAAACCTGTTTTGGGGC‑CACCTGTTCTTGGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGc < 1:1096823/99‑1 (MQ=255) gCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCg > 1:1067986/1‑99 (MQ=255) gCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCg > 1:1022132/1‑99 (MQ=255) cAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGt > 2:1642097/1‑99 (MQ=255) cAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGt > 1:989144/1‑99 (MQ=255) aCTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGttt > 2:756830/1‑99 (MQ=255) tGGTGTAAAGAAATAATCATTTTGTAAACCATTTTTGGGGT‑CAGCGTTTCTTCCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGa < 1:554515/99‑1 (MQ=255) gttgtCAAGGAATAAGAATTTTGTGAACCGGTTTGGGGGT‑CAGCTGTTTTTGGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGAc < 1:1406939/97‑1 (MQ=255) gtaaaGGTATAAGCATATTGTGAACCTGTTTTGGGGT‑CACCTGTTCTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCAt < 1:1546698/96‑1 (MQ=255) gtCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCAt < 2:162195/99‑1 (MQ=255) ggcAAGGAATAAGCATATTGTGACCCGTTTTTGGGGT‑CACCGTTTCTTCCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCAt < 1:153011/97‑1 (MQ=255) tCAAGGGATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCACCCATAAACAAAATCAGCGGCCGTTAACTTGGCTCCGTTTGACCAtt > 2:723460/1‑99 (MQ=255) tCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCAtt > 1:638297/1‑99 (MQ=255) aaGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCAtttt > 1:1070577/1‑99 (MQ=255) ggCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTaa > 2:372369/1‑99 (MQ=255) aaaaaGCTTTTTGTGACCCATTTTTGGGGT‑AAAGGTTTTTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAAcc < 2:1633217/96‑1 (MQ=255) aaaaGCATATTGT‑AACCCGTTTTTGGGTCCAACGTTTCTTCCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAAccc < 1:800586/97‑1 (MQ=255) aGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGtc > 2:1297996/1‑99 (MQ=255) aGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGtc > 1:11952/1‑99 (MQ=255) aGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGtc > 2:1599289/1‑99 (MQ=255) aGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGtc > 1:901533/1‑99 (MQ=255) gCATTTTGTGAACCAGTTTTGGGGT‑CAACGTTTCTTCCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGtct < 1:718666/99‑1 (MQ=255) cATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGtctc > 2:693751/1‑99 (MQ=255) ttGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCaa > 2:946153/1‑99 (MQ=255) ttGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCaa > 1:1695/1‑99 (MQ=255) ggacccATTTTTGGGGT‑CAACTTTTCTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGt < 2:1068842/95‑1 (MQ=38) cAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACAt > 1:820985/1‑99 (MQ=255) gTTTGGGGGT‑CAACGTTTCTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATta < 2:1155442/99‑1 (MQ=255) ttGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATaa > 2:796699/1‑99 (MQ=255) gT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGtttt > 1:1256608/1‑99 (MQ=255) cAAGGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTAc < 2:1297744/99‑1 (MQ=255) aCGGTTTTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCa < 2:119442/99‑1 (MQ=255) aCGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCa > 2:1580491/1‑99 (MQ=255) cGGTTTTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCAt < 2:911756/99‑1 (MQ=255) cGGTTCTTCCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCAt < 1:799117/99‑1 (MQ=255) tttttttCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCATCt < 1:207266/95‑1 (MQ=255) tttCTTTCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCATCt < 2:1590531/98‑1 (MQ=255) ttcttttCCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCATCt < 1:1473338/95‑1 (MQ=255) gTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCATCt > 1:283376/1‑99 (MQ=255) tCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCATCttt < 2:114447/99‑1 (MQ=255) | ATCAACTGTTTTTTTACCGGCCATTATTT‑CGCTGGCATTAGCAACTGGTGTCAAGGCATAAGCATATTGTGAACCAGTTTTGGGGT‑CAACGGTTCTTTGCCAGCCATAAACAAAATCAGCGGCCGTTAACTTAGCTCCGTTTGACCATTTTAACCCGTCTCGCAAAGTTACATTATAAGTTTTACCATCTTT > CP002094/500162‑500353 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |