Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,245,234 | 0 | G | T | 21.7% | 247.4 / 24.9 | 83 | L25F (TTG→TTT) | ycgY | uncharacterized protein YcgY |
Reads supporting (aligned to +/- strand): ref base G (15/50); new base T (18/0); total (33/50) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.39e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAGATTA > NC_000913/1245186‑1245282 | cAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTgg < 1:1532723/50‑1 (MQ=255) cAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTgg < 1:4870937/50‑1 (MQ=255) aGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGc < 1:6652049/50‑1 (MQ=255) aGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGc < 1:885178/50‑1 (MQ=255) aGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGc < 1:4985441/50‑1 (MQ=255) aGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGAGAGGCGG‑TTTTGGc > 1:2786889/1‑50 (MQ=255) gTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGCt > 1:6118347/1‑50 (MQ=255) gTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGCt > 1:1055676/1‑50 (MQ=255) gTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGcc > 1:1469850/1‑49 (MQ=255) gTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCt < 1:6432840/50‑1 (MQ=255) gTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑GTTTTGcc > 1:5790248/1‑49 (MQ=255) tAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGGTg > 1:4647845/1‑50 (MQ=255) tAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTg < 1:4368651/50‑1 (MQ=255) aGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGCTg > 1:361481/1‑50 (MQ=255) aGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGATTTTGGcgc > 1:7071352/1‑48 (MQ=255) aGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTtgctgc > 1:4342065/1‑50 (MQ=255) aGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑GTTTTGGTgc > 1:710741/1‑50 (MQ=255) gCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGCTtcg > 1:5383831/1‑50 (MQ=255) gCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTgcg < 1:3934388/50‑1 (MQ=255) ccAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGGTgcgc > 1:2786458/1‑50 (MQ=255) ccAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTgcgc < 1:5671158/50‑1 (MQ=255) cAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCt < 1:5636404/50‑1 (MQ=255) aGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGCTGCGCtt > 1:1729801/1‑50 (MQ=255) aGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGCTGCGCtt > 1:717765/1‑50 (MQ=255) aGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCtt < 1:217527/50‑1 (MQ=255) tCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTAt < 1:1638964/50‑1 (MQ=255) tCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTAt < 1:6387915/50‑1 (MQ=255) aGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGGTGCGCTtttt > 1:1539500/1‑49 (MQ=255) aGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTAtt < 1:3606805/49‑1 (MQ=255) aGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTAtt < 1:6174490/49‑1 (MQ=255) aGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTAtt < 1:2848798/49‑1 (MQ=255) aGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTg < 1:3634074/50‑1 (MQ=255) aGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTg < 1:5980805/50‑1 (MQ=255) ggCGTTAACAATGGGGCGGATAGTGAGGCGG‑ATTTTGGTGCGGTTAttt > 1:6774626/1‑48 (MQ=255) gCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGCTGCGCTTTTTTaa > 1:1663939/1‑50 (MQ=255) gCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGGTGCGCTCAtttt > 1:2465788/1‑47 (MQ=38) gCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGaa < 1:3181144/50‑1 (MQ=255) gCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGaa < 1:3638015/50‑1 (MQ=255) gTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTTGCTGCGCTTTTTGaaaa > 1:3300658/1‑48 (MQ=255) gTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑a < 1:3582804/50‑1 (MQ=255) ttAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑At < 1:4286494/50‑1 (MQ=255) aaaaaTGGGGCGGATAGTGAGGCGG‑TTTTTGCTGCGCTTATTTAAA‑Attt > 1:5355267/1‑50 (MQ=255) aaaaTGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTa < 1:3978075/50‑1 (MQ=255) aaaaTGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTa < 1:4438422/50‑1 (MQ=255) aaaaTGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTa < 1:970079/50‑1 (MQ=255) aaaTGGGGCGGATAGTGAGGCGG‑TTTTTGCTGCGCCTTTTTAAC‑Atttt > 1:5418734/1‑48 (MQ=37) aaaTGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAg < 1:1762720/50‑1 (MQ=255) aaTGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGa < 1:646222/50‑1 (MQ=255) aTGGGGCGGATAGTGGGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGaa < 1:5191606/50‑1 (MQ=255) aTGGGGCGGATAGTGAGGCGG‑TTTTTGGTTGGCGTTTTGAAC‑AATTTGaa > 1:218259/1‑50 (MQ=255) aTGGGGCGGATAGTGAGGCGG‑TTTTTGCTGCCCTTTTTTAAC‑ATTTAGaa > 1:5195756/1‑50 (MQ=255) aTGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGaa < 1:2984509/50‑1 (MQ=255) ggggCGGATAGTGAGGCGG‑TTTTTGGCGCGCTTATTTAAC‑ATTTAAAAtt > 1:2593550/1‑50 (MQ=255) ggggCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTTTTGAACAATTTAAAAt > 1:4967839/1‑50 (MQ=37) gggCGGATAGTGAGGCGG‑TTTTTGCTGCGCTTATTGAAC‑ATTTAGAAtt > 1:3904887/1‑49 (MQ=255) gggCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAAtt < 1:3860606/49‑1 (MQ=255) gggCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTg < 1:4083159/50‑1 (MQ=255) ggCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGa < 1:3601961/50‑1 (MQ=255) gCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGAt < 1:2292070/50‑1 (MQ=255) aTAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAAt < 1:678317/50‑1 (MQ=255) tAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATc < 1:6054590/50‑1 (MQ=255) tAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATc < 1:3848833/50‑1 (MQ=255) gTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATcc < 1:3015332/49‑1 (MQ=255) gTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCt < 1:1857687/50‑1 (MQ=255) gTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCt < 1:172366/50‑1 (MQ=255) gTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCt < 1:5761751/50‑1 (MQ=255) tGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCtt < 1:5299192/50‑1 (MQ=255) aGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCt < 1:6222273/50‑1 (MQ=255) aGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCt < 1:4296465/50‑1 (MQ=255) aGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCt < 1:5684380/50‑1 (MQ=255) aGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCt > 1:180788/1‑50 (MQ=255) ggCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTg < 1:7244597/50‑1 (MQ=255) gCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTgg < 1:6991394/50‑1 (MQ=255) cGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTgg > 1:6778210/1‑49 (MQ=255) gg‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTgg > 1:652895/1‑48 (MQ=255) gg‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTgg > 1:1755992/1‑48 (MQ=255) gg‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTgg > 1:2819442/1‑48 (MQ=255) gg‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCa < 1:1862337/50‑1 (MQ=255) g‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTgg > 1:2476969/1‑47 (MQ=255) g‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAg < 1:6152070/50‑1 (MQ=255) g‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAg < 1:1595981/50‑1 (MQ=255) g‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAg < 1:4843859/50‑1 (MQ=255) ttttGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAgc > 1:6211059/1‑49 (MQ=255) ttttGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAGa < 1:919595/50‑1 (MQ=255) ttGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAGAt > 1:3343377/1‑49 (MQ=255) tGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAGAtt < 1:4959608/49‑1 (MQ=255) ggCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAGAtt > 1:7054836/1‑48 (MQ=255) ggCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAGATTa < 1:342894/49‑1 (MQ=255) ggCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAGATTa < 1:1086286/49‑1 (MQ=255) | CAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGG‑TTTTGGCTGCGCTTATTGAAC‑ATTTAGAATTGATCAATCCTTCTGGCAGATTA > NC_000913/1245186‑1245282 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |