Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,606,322 | 1 | . | G | 63.6% | 121.6 / 85.8 | 55 | noncoding (80/87 nt) | leuQ | tRNA‑Leu |
Reads supporting (aligned to +/- strand): ref base . (10/10); new base G (27/8); total (37/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.13e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGT > NC_000913/4606189‑4606450 | atCGCTTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGCACACTACCACATGAAGCTAGCGCGTCTACCATTTCCGCCACCTTCGCCCAGTGCTAGCAATCTCATCGTGGTTTTTGGTGCG‑A‑‑ggggg‑gg < 2:1007207/141‑1 (MQ=18) cGCCACGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑ggggg‑gggg < 2:1369195/141‑2 (MQ=255) ggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTa > 1:1706175/1‑141 (MQ=255) gACTGGAACCCCCACGTCCGTAAGAACACTACCACCTGCATCTCGCGTGTATACAATTTCCCCCACCGTCGCCCAGTGCGAGCAATTACAACGTGTTATTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGa < 2:285805/141‑1 (MQ=25) ccccccCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATACCGCCACCTTCGCCCAGTGCGAGCATTATCAACGTGGTTTTTGGTGCGTG‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAac < 2:1303586/141‑1 (MQ=37) cccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAac < 2:1147438/141‑1 (MQ=255) ccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAaca > 2:1323302/1‑141 (MQ=255) cccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAacac > 1:618217/1‑141 (MQ=255) aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCCATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCAGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa > 1:57624/1‑141 (MQ=255) acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa < 1:1387876/141‑1 (MQ=255) acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa < 1:1387902/141‑1 (MQ=255) acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa < 1:1388126/141‑1 (MQ=255) acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa < 1:1494391/141‑1 (MQ=255) acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt > 1:1085050/1‑141 (MQ=255) cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc > 1:740550/1‑141 (MQ=255) cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc > 1:195810/1‑141 (MQ=255) tAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc > 1:1108668/1‑141 (MQ=255) tAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc > 1:1396362/1‑141 (MQ=255) gcgcGTCTACCCATTCCTCCACCTTCGCCCATTGCGAGAAATTTCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc < 1:792792/141‑1 (MQ=18) gcgcGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt > 1:656423/1‑141 (MQ=255) gcgcGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc > 1:1367985/1‑141 (MQ=255) gcgcGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc > 1:1367993/1‑141 (MQ=255) gcgcGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc > 1:512457/1‑141 (MQ=255) gcttCTACCAAGTCCTCAACCTTCGCCCAGTGCCAGCAATATCACCGGGGTTTGTGGTGCG‑AGGGGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc < 2:1958812/138‑1 (MQ=17) gtCCGCCACCTTCGCCCCGTGCTAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAt < 1:1765522/140‑1 (MQ=255) gtCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGCTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAt < 1:1302120/140‑1 (MQ=255) ccGCCACCTTCGCCCAGTGCTAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAtt < 2:1592847/141‑1 (MQ=255) ccGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAtt < 2:741545/141‑1 (MQ=255) gCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCt > 1:203658/1‑141 (MQ=255) ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaaa > 1:178642/1‑141 (MQ=255) ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaaa > 1:178678/1‑141 (MQ=255) cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc > 1:1597085/1‑141 (MQ=255) cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc > 1:291648/1‑141 (MQ=255) cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc > 1:291667/1‑141 (MQ=255) aGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCt > 1:693707/1‑141 (MQ=255) gTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCt > 1:875422/1‑141 (MQ=255) gTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCt > 2:450903/1‑141 (MQ=255) gTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCt > 2:447267/1‑141 (MQ=255) gTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCt > 1:1176438/1‑141 (MQ=255) gTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACGTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCt > 2:1777167/1‑141 (MQ=255) gTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTa > 1:658883/1‑141 (MQ=255) gAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTGGCATACCATCAAGTCTTAAAAAGAATTGCTACCa > 2:1664826/1‑141 (MQ=255) gAGCAATATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCa > 2:1664809/1‑141 (MQ=255) aaTATCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGa > 1:1136/1‑141 (MQ=255) tatCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg > 1:1731785/1‑141 (MQ=255) tatCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg > 2:1531966/1‑141 (MQ=255) tatCAACGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg > 1:1731840/1‑141 (MQ=255) aCGTGGTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAt > 1:1414052/1‑141 (MQ=255) aCGTGGTTTTTGGTGCG‑A‑‑GGGGGTGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAt < 1:1627555/141‑1 (MQ=255) aCGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAtt > 1:1975831/1‑141 (MQ=255) ggTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa > 2:859683/1‑141 (MQ=255) ggTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa > 2:859701/1‑141 (MQ=255) ggTTTTTGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa > 2:931611/1‑141 (MQ=255) tttttGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt < 2:656423/141‑1 (MQ=255) tttttGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt < 1:967790/141‑1 (MQ=255) tttttGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt < 1:567570/141‑1 (MQ=255) attttGGTGCG‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt < 1:966085/140‑1 (MQ=255) | ATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGT > NC_000913/4606189‑4606450 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |