Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,969,794 | 0 | A | C | 40.8% | 76.7 / 10.3 | 49 | G397G (GGT→GGG) | tap | methyl‑accepting protein IV |
Reads supporting (aligned to +/- strand): ref base A (24/5); new base C (1/19); total (25/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.98e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCC > NC_000913/1969647‑1969941 | ggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCt < 2:287246/151‑1 (MQ=255) aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCa < 1:1098446/151‑1 (MQ=255) aCAATATCGATCATGGTCGGGGCGGTGTTATTCACCAGTTCCGACCCCTGCGGCACACGATTGACTGCCTCTCCGATCACCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGAAA‑GGTTG‑CGCACTCCACCTGCCACTACCGCaa < 1:456199/151‑1 (MQ=255) cAATATCGATCATGGTCGCGGCGGCTTTATTCACCAGTTCCGACCCCTGCTGGACACGATTGACTGCCTCTTCGATCAGCCCTTTGCTCTCTTTTGCTGCCTGCGCGCTGCGGCGGGCAA‑GATTG‑CGCCCTTCCCCTGCCACTa < 1:181600/144‑1 (MQ=255) aTCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGcc > 2:963057/1‑151 (MQ=255) cATGGTCGCGGCGTCTTTATTCCCCAGTTTCGACCCCGGCTGGACACGTTTGCCTGCCTCTTCGATCCGCCTTTGGCTCTCTTTCGCTCCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCCCTTCCCCTGCCACTACCGCAAAACCAc < 1:128849/146‑1 (MQ=255) aTGGTCGCGGCGGCGCTATTCACCAGTTTCGAACCCTGCTGGACACGCTTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGcccc > 1:895305/1‑150 (MQ=255) gTCGGGGCGGCTTTATTCCCCAGTTTCGACCCCTGCTGGACACGTTTGACTGACTCTCCGATCAGCCCTTTGCTCTTTTTGCCTCCCGGCGCTCTACGGCGGGCAA‑GTTGG‑CGCTCTTCACCTGCCACTACCGCAAAACCACGCCCCTg < 1:417791/149‑1 (MQ=255) cgcgGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACCGACTCTTCGCTCAGCCCTTTGATCTCTTTTGCTCCCTCCGCGCTACGGCCGGCAC‑GATTG‑CGCCCTTCACCTGCCACTACCGCAAAACCACGCCCCTGccc > 1:441312/1‑148 (MQ=255) gcggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCtctc > 2:820823/1‑151 (MQ=255) cggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTAGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc > 1:264237/1‑151 (MQ=255) ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg > 1:157822/1‑149 (MQ=255) ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcg > 1:611574/1‑147 (MQ=255) ttCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GCTTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGc > 1:690634/1‑149 (MQ=255) aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc > 2:132259/1‑150 (MQ=255) tttCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCAc > 1:871258/1‑150 (MQ=255) aCCCTCTTGGCCACCATTGTCTGCCTTTACTACCCGCCCTTTGGTGTCTTTTCCTTCTTGCGC‑CCACCGCCGGCAAGGTTTG‑CGCTTTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGTCGCTTCCACTGCCg < 2:446870/149‑1 (MQ=255) gctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GCTTG‑CGCACTTCACCTGCCACTACAGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCAACTGCCGCATTc > 1:962107/1‑149 (MQ=255) ctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAgg > 2:107255/1‑151 (MQ=255) gaccacttttGCCTGCCTCTCCGACCACCCCTTTGATCTCTTTTGCTCCCGGCCCGCAACGGC‑GGCAA‑TGTTGCCCCACTTCCCCTGCCCCTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc < 1:747209/138‑1 (MQ=255) gaacacttttGACTGCCTCTTCGATCACCCTTTTGATCTTTTTTGCGCCCTGCGCGCTCCGCCGGGAAA‑GATTG‑CGCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAggg < 2:802888/142‑1 (MQ=255) aTTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCGGGCAA‑GATTG‑CGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa < 2:1083750/148‑1 (MQ=255) ttGACTGACTTTTCGATCACCCTTTTGATCTTTTTTGCTCCCGGCGCGCTACGGGGGGCAA‑GATTG‑CACACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGaa < 2:1000643/148‑1 (MQ=255) ttGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAatat > 2:59464/1‑151 (MQ=255) tgtctctTCGACCCGCCCTTTGCTCCTTTTTGCCCCCCGGGCTCTACGGCGGGAAA‑GATGG‑CGCACTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTg < 2:96350/145‑1 (MQ=21) tgcctctTCGACCACCCTTTTTATCTTTTTTGCTCCCTGCGCGCTACGGCTGGCAA‑GTTTG‑CCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt < 2:883278/144‑1 (MQ=25) ctctTCGATCCCCCCTTTGATCTCTTTTTCTGCCGCCGCGCTACGGCGGGAAA‑GATTG‑CGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCtt < 1:791046/150‑2 (MQ=25) ttCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa > 1:921393/1‑150 (MQ=255) tccGATCACCCTTTTATTCTTTTTTGCTGCCTGCGCCCTCCGGCTGCAAA‑GTTGG‑CCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTc < 2:706688/143‑1 (MQ=21) gATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc > 1:858138/1‑150 (MQ=255) aTCACCCCTTTGATCTTTTGTGCTCCCGGCGCGCTACGGCGGGAAA‑GTTTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa < 2:848152/150‑1 (MQ=25) ccgCCTTTTGATCGCTTTGCCGGCCTGCGCGCTGGGGGGGGAAG‑GTTGG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt < 1:1024637/149‑1 (MQ=18) tttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg > 2:434953/1‑151 (MQ=255) tttGATCTCTTTTGCTCCCGGCGCGCTCCGGCGGGCAA‑GATTG‑CCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc < 2:634804/150‑1 (MQ=37) tcttTTTTGCTCCCTGCGCGCTACGGCTGGAAA‑GATTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa < 2:817989/150‑1 (MQ=37) ttttttGCTGCCTGCGCGCTCGGGCTGGCAA‑GTTTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa < 2:292604/143‑1 (MQ=37) tctTTTCCTCCCGGCGCGCTCCGGCGGGAAA‑GATTG‑CCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa < 2:762698/148‑1 (MQ=25) ctTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt > 1:1050037/1‑151 (MQ=255) ttttGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg > 2:970191/1‑151 (MQ=255) tttGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg > 1:470351/1‑150 (MQ=255) ttcctgccggcGCGCTCCGGCGGGCAA‑GATTG‑CGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGa < 2:390227/147‑1 (MQ=37) ttGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGAt > 1:50719/1‑151 (MQ=255) cgcgcTCCGGCTGCCAA‑GATTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg < 2:567191/148‑1 (MQ=255) tACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt > 1:69045/1‑151 (MQ=255) cGGCTGGAAA‑GTTTG‑CCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt < 2:819186/149‑1 (MQ=25) cAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt > 1:860424/1‑151 (MQ=255) cgcACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCg > 2:124251/1‑150 (MQ=255) cACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtc > 2:362177/1‑151 (MQ=255) ttCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTcc > 1:634698/1‑151 (MQ=255) | GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCC > NC_000913/1969647‑1969941 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |