| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,648,148 | 0 | A | C | 27.6% | 52.1 / 13.0 | 29 | G1380G (GGT→GGG) | yfhM | alpha2‑macroglobulin |
| Reads supporting (aligned to +/- strand): ref base A (21/0); new base C (0/8); total (21/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.33e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACC > NC_000913/2648081‑2648213 | gTAGCAGTTTATTTTCTTCCAGCAAGGAGGGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCa < 1:1287727/68‑1 (MQ=255) gCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCAccc > 1:201148/1‑68 (MQ=255) aGTTTGTTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAg > 1:2437410/1‑68 (MQ=255) gTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGc > 1:2035472/1‑68 (MQ=255) ttcttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATcc > 1:2357098/1‑68 (MQ=255) ccAGCAAGGGGAGCATTAGCGCGTTGCCGCGCAGTGACCTACCGAAATCCCCCAGCCATATCCGCtca < 1:1742790/68‑1 (MQ=255) ccAGCAAGGAGAGCTTTAACGCGTTTTCGGGCAGTAAACTACCGTAATCCCCCAGCCATATCCGCtca < 1:582893/68‑1 (MQ=255) ccAGCAAGGAGAGCTTGAACGCGTTGGCGCGCAGTAAACTCCCGTAATCCCCCAGCCATATCCGCtca < 1:3142780/68‑1 (MQ=255) ccAGCAAGGAGAGCATTACCGCTTTGTCGCGCAGTAAACTACCGTAATCCCCCAGCCATATCCGCtca < 1:2343565/68‑1 (MQ=255) ccAGCAAGGAGAGCATTAACGCTTTGTCGCGCATTGAACTCCCGTAACCCCCCAGCCATATCCGCtca < 1:1711486/68‑1 (MQ=255) ccAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTCCCGTAATCCCCCAGCCATATCCGCtca < 1:2273336/68‑1 (MQ=255) ccAGCAAGGAAGGCATTAACGAGGTGGCGCGCAGTGAACTACCTTAACCCCCCAGCCATATCCGCtca < 1:4156331/68‑1 (MQ=255) agaGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTAAACTCCCGTAATCCCCCAGCCATATCCGCtca < 1:869308/66‑1 (MQ=255) aGCAAGGAGAGCATTAACGCTTTGCCGCGCAGTGACCTACCGTAATCCCCCAGCCATATCCGCtcatc < 1:2162303/68‑1 (MQ=255) gCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtcatca > 1:2220083/1‑68 (MQ=255) agggAGAGCATTACCGCGTTGCCGCGCAGTAACCTCCCGTAATCCCCCAGCCATATCCGCTCATCACt < 1:440655/66‑1 (MQ=255) aaGGATGGCATTACCGCGTTGCCGGGCAGTAAACTACCGTAATCCCCCAGCCATATCCGCTCATCACt < 1:1143228/68‑1 (MQ=255) aaGGAGCGATTTAACGGGTGGGCGGGAAGGGAACTACGGTAATCCTCCAGCCATATCCGCTCATCACt < 1:2810076/68‑1 (MQ=255) aaGGAGAGGTGTAACGGGTTGTCGCGCAGTGAACTACCGAACTCCCCCAGCCATATCCGCTCATCACt < 1:3290751/68‑1 (MQ=255) aaGGAGAGCATTACGGCTTTTTCGCGCAGGAATCTCCCGTAATCCCCCAGCCATATCCGCTCATCACt < 1:3182325/68‑1 (MQ=255) aaGGAGAGCATTAACGCTTTGTCGCGCAGTAAACTCCCGTAACCCCCCAGCCATATCCGCTCATCACt < 1:635742/68‑1 (MQ=255) aaGGAGAGCATTAACGCGTTGCCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACt < 1:908958/68‑1 (MQ=255) aaGGAGAGCATAACCGTGTGGCCGCGCAGTGAACTCCCGTAACCCCCCAGCCATATCCGCTCATCACt < 1:782240/68‑1 (MQ=255) ggAGAGATTTAGCGGTTTGTGGCGAAGTAACCTCCCTAAATCCCCCAGCCATATCCGCTCATCACTAt < 1:2696951/68‑1 (MQ=255) agagCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATAACTATTa > 1:564267/1‑68 (MQ=255) aTTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCAATATCCGCTCATCACTATTACgcgg > 1:4058370/1‑68 (MQ=255) cgcgTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGttt > 1:148968/1‑68 (MQ=255) cgcgTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGttt > 1:2074641/1‑68 (MQ=255) cgcgTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGttt > 1:3964268/1‑68 (MQ=255) gTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTca > 1:4028307/1‑68 (MQ=255) gTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTca > 1:3071362/1‑68 (MQ=255) gcAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagc > 1:1072848/1‑68 (MQ=255) ccAGTAAACACCCGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagc < 1:3023012/67‑1 (MQ=255) cTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGc > 1:3191942/1‑68 (MQ=255) aCCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGctt > 1:2851022/1‑68 (MQ=255) cGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGcttct > 1:543033/1‑68 (MQ=255) cGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGcttct > 1:3589350/1‑68 (MQ=255) aaTCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTca > 1:2001403/1‑68 (MQ=255) aaTCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTAca > 1:3417879/1‑68 (MQ=255) tCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcacc > 1:402663/1‑68 (MQ=255) tCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcacc > 1:2621614/1‑68 (MQ=255) tCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcacc > 1:1935863/1‑68 (MQ=255) tCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcacc > 1:1073985/1‑68 (MQ=255) | GTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACC > NC_000913/2648081‑2648213 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |