| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,357,894 | 0 | C | A | 40.0% | 20.8 / 16.7 | 25 | pseudogene (79/186 nt) | ypaA | protein YpaA |
| Reads supporting (aligned to +/- strand): ref base C (12/3); new base A (0/10); total (12/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.08e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGG > NC_000913/2357830‑2357957 | aCGGCTGCGGCAGGGAGGACAACAAATTGGCTGGCCGGGAGGTAAATTTGAAGGTTTGCATGAAc > 1:4569536/1‑65 (MQ=25) gaaggaCATCAAATTGGCTGGCAGGAAGGGAAATTAGAAGGTTTGCATGAACAAGCCATTAAAAt > 1:3196034/1‑65 (MQ=255) aggaCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTg > 1:3689479/1‑65 (MQ=255) ggaCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGc > 1:3854022/1‑65 (MQ=255) ggaCATCAAATTGGCTGGCAGGAAGGGAAATTAGAAGGTTTGCATGAACAAGCCCTTAAAATTGc > 1:2352853/1‑65 (MQ=38) ggaCATAAAGTTGGGTGGGAGGAAGAAAAATTAGAGGTTTTGCATAAAAAAGCCATTAAAATTGc < 1:3700181/65‑1 (MQ=16) gaaCATAAAATTGGCGGGAAGGAAGGTAAATAAGAAGGTTTGCATGAACAAGCCATTAAAATTGc < 1:1865787/63‑1 (MQ=21) aCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCtt > 1:1261214/1‑65 (MQ=255) taaaaTTGGCTGGCGGGAAGGTAAATAAGAGGGTTTGCATAAAAAAGCCATTAAAATTGCTTTgc < 1:4063325/63‑1 (MQ=25) tCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgc > 1:2464737/1‑65 (MQ=255) tCAAATGGGCGGGAAGAAGGGTAAATTAGAAGTTTTGCAAAAAAAAGCCATTAAAATTGCTTTgc < 1:4279733/65‑1 (MQ=16) ttGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATg > 1:1413397/1‑65 (MQ=255) ttGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAAAAAGCCATTAAAATTGCTTTGCGCATg < 1:1948157/65‑1 (MQ=255) ttGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATg < 1:1519899/65‑1 (MQ=255) ttGGCTGGCAGAAAGAAAAATTAGAAGGTTGGCAAAAAAAAGCCATTAAAATTGCTTTGCGCATg < 1:1172760/65‑1 (MQ=18) ttGGCTGGCAGAAAAGTAATTTAGAAGTTTTCCATAAACAAGCCATTAAAATTGCTTTGCGCATg < 1:2190488/65‑1 (MQ=18) ttGGCTGGAAGGAAGGAAAATTAGAAGGTTTGCATAAACAAGCCATTAAAATTGCTTTGCGCATg < 1:108793/65‑1 (MQ=21) tgggctggcAGGAAGGTAAATTAGAAGTTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATg < 1:2819623/63‑1 (MQ=25) tgggcgggaAGGAAGTTAAATTAGAAGGTTTGAAAAAAAAAGCCATTAAAATTGCTTTGCGCATg < 1:1548098/63‑1 (MQ=18) ggctggcAGGAAGGTAAATTAGAAGGTTTGAATGAAAAAGCCATTAAAATTGCTTTGCGCATGCt < 1:72915/65‑1 (MQ=25) ggcgggcgGGAAGATAATTTAGAAGTTTTGAATAAAAAAGCCATTAAAATTGCTTTGCGCATGCt < 1:1166995/65‑1 (MQ=12) ggcgggcAGGAAGGAAAATAAGAGGGTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATGCt < 1:951261/65‑1 (MQ=18) ggcgggaAGGAAGTTAAATTGAAAGGTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATGCt < 1:4141956/65‑1 (MQ=17) ggcgggaAGGAAGGAAAATAAGAAGGTTTGAAAAAAAAAGCCATTAAAATTGCTTTGCGCATGCt < 1:2925745/65‑1 (MQ=12) ggcaggcaGGAATATAAATTAGAGGTTTTAAATGAACAAGCCATTAAAATTGCTTTGCGCATGCt < 1:4149138/65‑1 (MQ=17) tgggAGGAAGGAAAATTAGAAGGTTAAAATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGa < 1:473480/65‑1 (MQ=17) tggcAGGAAGGTAAATTAGAAGGTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGa < 1:3580910/65‑1 (MQ=255) gcAGGAAGGTAAATTAGAAGTTTTGCAAAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc < 1:4262085/65‑1 (MQ=21) gcAGGAAGGTAAATTAGAAGGTTTGAATGAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc < 1:67678/65‑1 (MQ=255) ggaaggaaAATAAGAAGGTTTGAATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAgg < 1:1826199/65‑1 (MQ=17) ggAAGGTAAATTAAAGGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAgg < 1:3651968/65‑1 (MQ=255) gacggTAAATTGAAGGGTTTGCAGGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg < 1:9158/62‑1 (MQ=14) gAATGTAAATTAGAAGGTTCGCATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg < 1:3630763/65‑1 (MQ=18) gAAGGTAAATTAGAAGGTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg < 1:117187/65‑1 (MQ=255) aaaGTTAAATTAGAGGTTTTGAATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg < 1:387362/64‑1 (MQ=17) aaaGATAAATAAGAAGTTTTGAATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg < 1:4529459/64‑1 (MQ=17) ggTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCtt > 1:143418/1‑65 (MQ=255) ggTAAATTAGAAGGATTGCATGAACAAGAAATTCAAATTAATTTGTGCATGCAGGAGCAGGGCtt > 1:4567157/1‑65 (MQ=9) ggTAAATTAAAGGTTTTGCAAAAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCtt < 1:3021995/65‑1 (MQ=17) aaaTTAGAAGGTTTGAATAAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGa < 1:2484262/65‑1 (MQ=255) aaTTAGAGGTTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt < 1:838074/65‑1 (MQ=17) aaTTAGAAGGTTTGAAAAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt < 1:3365152/65‑1 (MQ=17) aaTAAGAAGTTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt < 1:3992722/65‑1 (MQ=14) agtAGAAGTTTTGAATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATc < 1:1768314/63‑1 (MQ=17) ttAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATcc > 1:3381482/1‑64 (MQ=255) tAGAAGGTTTGCATAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGt < 1:1538586/65‑1 (MQ=255) aaGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGAc > 1:147225/1‑65 (MQ=255) ggTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa > 1:4141125/1‑65 (MQ=255) ttttGCAAAAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAg < 1:3094578/64‑1 (MQ=17) ttGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGt < 1:1558709/65‑1 (MQ=255) cATGAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCt < 1:4242684/65‑1 (MQ=255) aTGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTc > 1:3127239/1‑65 (MQ=255) aCAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCgg > 1:2592506/1‑65 (MQ=255) | ACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGG > NC_000913/2357830‑2357957 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |