Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,969,794 | 0 | A | C | 42.5% | 62.1 / 15.7 | 41 | G397G (GGT→GGG) | tap | methyl‑accepting protein IV |
Reads supporting (aligned to +/- strand): ref base A (22/1); new base C (0/17); total (22/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.03e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTG > NC_000913/1969651‑1969943 | aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGcc > 2:262686/1‑150 (MQ=255) aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGcc > 1:924419/1‑150 (MQ=255) aCCGAACTGACAATAGCGATCATGGTCGCGGCGGCGTTATTCTCCATTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCGGCCTGCGCGCTCCGGCTGGC‑AA‑GATTG‑CGCACTTCGCTTGCCa < 1:883783/151‑1 (MQ=255) cGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACt > 1:892406/1‑151 (MQ=255) tatCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAAcc > 2:292318/1‑151 (MQ=255) atCGATCTGGGTCGCGCCGGCGTTATCCACCATTTTCGACCCCTGCTGGACACGATTGCCTGACTCTTCGATCAGCCCTTTGTTCTCTTTTGCTGCCTGCGCGCTCCGGCTGGC‑AA‑GATTG‑CGCACTTCCCCTGCCACTACCGCAAAAcc < 1:406649/150‑1 (MQ=255) tGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTg > 1:855065/1‑151 (MQ=255) cgcgGGGGCTTTATTCCCCAGTTTCAACCCCTGCTGGACACGTTTGACTGACTCTTCGATCACCCTTTTGATCTTTTTTGCTGCCGGCGCGCTCCGGCTGGC‑AA‑GATGG‑CGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCt < 2:559933/146‑1 (MQ=255) gcggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCtctc > 2:457688/1‑151 (MQ=255) ggcggctTTATTCCCCAGTTCCAACCCCTGCTGGCCACGTTTGACTGTCTCTTCGATCGGCCCTTTGATCTTTTTTGCTGCCGGCGCGCTACGGCTGGC‑AA‑GATGG‑CGCCTTTCCCCTGCCACTACCGCAAAACCACGCCCCTGct < 2:857727/146‑1 (MQ=255) tATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑TGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGc > 1:33888/1‑151 (MQ=255) tCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGc > 1:851483/1‑151 (MQ=255) ccAGTTTCCACCCCGGCTGGCCACGATGTCCTGCCCCTTCGACCAGCCCTTTGTTTGTTTTTGCCCCCGGCCCGCTCCGG‑TGGCAAA‑GTTTG‑GGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc < 1:787615/146‑1 (MQ=255) cAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc > 1:289313/1‑151 (MQ=255) gAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg > 1:86716/1‑151 (MQ=255) tgctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc > 2:587194/1‑150 (MQ=255) gctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc > 1:950430/1‑149 (MQ=255) acccGATTGCCTGCCTCTTCATTCACCCCTTTTATCTCTTTGCCTCCCGGCGCGCCCCGGCTGGA‑AA‑GTTGG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCa < 1:426805/142‑1 (MQ=255) ttGACTGACTTTTCGCTCACCCTTTTGATCTTTTTTGCCCCCGGCGCGCTACGGCTGGA‑AA‑GTTTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAata < 1:10709/150‑1 (MQ=25) ctgcctCTCCGATCGGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCGGGA‑AA‑GATTG‑CGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTgg < 2:81674/150‑1 (MQ=37) actCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTg > 2:953661/1‑151 (MQ=255) tctTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTgg > 1:713447/1‑150 (MQ=255) tCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa > 1:890713/1‑149 (MQ=255) aTCAGCCTTTTTCTCTTTTTTTTCCCCTGCGCGCTACGGCTGGC‑AA‑‑ATTGCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCaa < 1:514269/151‑1 (MQ=18) ccttttGTTCTTTTTTGCTCCCGGCGCGCTACGGCTGGA‑AA‑GATTG‑CCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCg < 2:6232/148‑1 (MQ=25) ccttttGTTCTTTTTTCCTGCGTCCGCGCACCGGC‑GGC‑AAGGTTTG‑CCCCCCTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg < 1:778105/140‑1 (MQ=18) ctctTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACg > 1:204068/1‑151 (MQ=255) ttttGCTGCCTGCGCGCTCCGGCTGGA‑AA‑GTTTG‑CGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt < 2:289553/150‑1 (MQ=255) tttGCCGCCGGCGCGCTCCGGCGGGA‑AA‑GATGC‑CCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg < 1:150842/150‑1 (MQ=25) tgcctgcGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt > 1:843654/1‑150 (MQ=255) cctgcGCGCTACGGCTGGC‑AA‑GTTTG‑CGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGAt < 2:374954/145‑1 (MQ=255) cccgcgcccTACGGCTGGC‑AA‑TTTCG‑GGCCCTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt < 2:579948/145‑1 (MQ=25) gcgcTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCa > 2:475104/1‑150 (MQ=255) cGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt > 2:90984/1‑151 (MQ=255) aa‑GTTTG‑CGCTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtg < 2:203085/151‑1 (MQ=255) aa‑GTTTG‑CGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGc < 2:892538/145‑1 (MQ=255) aa‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtg > 1:934394/1‑151 (MQ=255) ggttgg‑cgcCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt < 2:643226/137‑1 (MQ=255) cTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCg < 2:543359/146‑1 (MQ=255) ttCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTcc > 1:237778/1‑151 (MQ=255) cACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTg > 2:696876/1‑151 (MQ=255) | ACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGC‑AA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTG > NC_000913/1969651‑1969943 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |