| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_007779 | 3,018,219 | 0 | T | G | 65.0% | ‑1.7 / 34.7 | 40 | S135A (TCG→GCG) | ssnA | chlorohydrolase/aminohydrolase |
| Reads supporting (aligned to +/- strand): ref base T (2/12); new base G (26/0); total (28/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.63e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-05 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGA > NC_007779/3018124‑3018319 | tctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGctc > 1:1134068/1‑101 (MQ=255) aTTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTcc < 3:1469415/101‑1 (MQ=255) aTTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTcc > 3:1128725/1‑101 (MQ=255) ttACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGGTccccc > 1:473149/1‑100 (MQ=255) tACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTcccc > 1:30281/1‑99 (MQ=255) cAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACAt < 3:1452642/101‑1 (MQ=255) cGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCcccacttc > 3:1555770/1‑97 (MQ=255) ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACATTgcg > 3:126904/1‑101 (MQ=255) ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCGCTcccccattcg > 3:1379708/1‑93 (MQ=255) gACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTgcgc > 3:1444322/1‑101 (MQ=255) aCTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGcgcg < 1:351796/101‑1 (MQ=255) tGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATGGCGCG‑Ac < 4:160991/101‑1 (MQ=255) ttGTTCCTTGGAAGCGATTAAGCGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCGCTCCACATTGCGCG‑ACGcct > 4:134540/1‑99 (MQ=255) tGTTCCCTGGAAGCGATTCAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACCTTGCGCG‑GCGCAtt > 4:83738/1‑101 (MQ=255) tGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACATTGCGCG‑ACGccct > 1:702448/1‑98 (MQ=255) gTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTTCGCG‑ACGCAttt > 3:367370/1‑101 (MQ=255) tCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCCTTGCGCG‑GCGCGctttt > 1:578042/1‑101 (MQ=255) ccTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTaa < 3:465740/101‑1 (MQ=255) tGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACCTTGCGCG‑CCGCATTTTTtaaa > 3:957876/1‑101 (MQ=255) ttACGAGCGGATGTACCTCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑g > 4:154006/1‑101 (MQ=255) ggATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGGCCTCCCCCTTTCGCGCCCCCCTTTTTAAAAATTGGCCTCgcgcgcg > 3:1085397/1‑101 (MQ=255) aTGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGa < 3:907585/101‑1 (MQ=255) aCATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTGGTCCT‑GCGNGCTAGGACCTg < 4:20072/101‑1 (MQ=255) aTCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGtt < 3:1074932/101‑1 (MQ=255) aTCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTGCGCG‑GCGCCTTTTTTTAAAATGGGCC‑GCGCGCGCTGACCTGtt > 3:1428362/1‑101 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTGCGCG‑ACGCATTTTTTAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTg > 1:958708/1‑101 (MQ=255) tATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGaaa < 1:244620/101‑1 (MQ=255) aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAac < 1:1430425/101‑1 (MQ=255) gATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACCTTGCGCG‑ACGCATTTTTTAAAGTTGGGCT‑GCGCGCGATGACCTGTTTTTTAactac > 3:823526/1‑101 (MQ=255) aTCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACATTGCGCG‑ACGCATTTTTTAAAATTGGCCT‑GCGCGCGATGACCTGTTTTGAAactact > 3:1523054/1‑101 (MQ=255) cACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCCTTGCGCG‑ACGCATTTTTTAAAATTTGGCT‑GCGCGCGATGACCCGTTTTTTAactactca > 3:344265/1‑99 (MQ=255) cACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCCCCCCCTTGCGCG‑GCGCACTTTTTAAAATTGGCCT‑GCGCGGGGTGGCCCGTTTTTTAACTACTGa > 3:1367198/1‑101 (MQ=255) cATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCAATTCGCG‑ACGCATTTTTTAAAATTGGGCC‑GCGCGCGATGACCCGTTTTGAAACTACTGACCg > 1:280350/1‑101 (MQ=255) aTGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTGCGCG‑ACGCATTTTTTAAAATTGGCCT‑GCGCGCGGTGACCTGTTTTTAAAATACTGACCGt > 1:1074745/1‑101 (MQ=255) gTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCa > 1:473469/1‑101 (MQ=255) gTATATCGGCGGGTCCCTCTCCACATTGCGCG‑ACTCTTTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCa < 3:709479/101‑1 (MQ=255) tataTCGGCGGGGCGCTCTCCACATTGCGCG‑ACGCATTTTTTAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTTAAACTACTGACCGTAACAACGGCATCaa > 1:875878/1‑101 (MQ=255) tatCGGCGGGGCGCTCCCCCCAATGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAag > 4:40050/1‑101 (MQ=255) cggcggGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGt < 1:1118114/101‑1 (MQ=255) gTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAACTACTGAGCGTAACAACGGCATCAAAGAGTTGCAgg > 3:10053/1‑101 (MQ=255) tCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa > 3:1019732/1‑101 (MQ=255) | TCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCG‑ACGCATTTTTAAAAGTTGGCCT‑GCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGA > NC_007779/3018124‑3018319 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |