Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,393,061 | 0 | T | G | 33.3% | 106.1 / 35.9 | 67 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (6/38); new base G (22/0); total (29/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.07e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392986‑4393136 | aaaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggt < 1:1729352/76‑1 (MQ=255) aaaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggt < 2:407496/76‑1 (MQ=255) aaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtg < 1:1036431/76‑1 (MQ=255) aaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtgg < 2:106429/76‑1 (MQ=255) gCAGGGCGGTTCCGATGTTATGAGCAACGCCCCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTgggg < 2:1438702/76‑1 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGgggggggg > 2:1502121/1‑74 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGgggggggg > 2:1194185/1‑74 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGgggggggg > 2:1117206/1‑74 (MQ=255) gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑t < 1:565409/76‑1 (MQ=255) gCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑Taa < 1:1117206/76‑1 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATg < 1:1784085/76‑1 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGG‑TAAATg > 2:1815980/1‑76 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGG‑TAAATgg > 1:38509/1‑76 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGG‑TAAATgg > 2:664173/1‑76 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGA‑TTAAtag > 2:1801919/1‑74 (MQ=255) tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGt < 1:566189/76‑1 (MQ=255) tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGC‑TAAATggg > 2:439803/1‑75 (MQ=255) tCCGATGTTATGAGCAACACCACACGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGt < 1:1868408/76‑1 (MQ=255) cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGCTGGCTCTTATCGGCTGGGGGGGGG‑TAAATGGttt > 1:372213/1‑76 (MQ=255) cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGG‑TAAAAGGgtt > 1:965262/1‑76 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGtttt < 2:1473379/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGtttt < 2:226829/76‑1 (MQ=255) aTGTTATGAGCAACCCCCCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGttttt < 2:235377/76‑1 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCA‑TAAATGGttttt > 1:1059853/1‑76 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTtc < 1:249364/76‑1 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTtct < 1:819023/76‑1 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGA‑TAAATGGGTTTtct > 2:974630/1‑76 (MQ=255) ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGG‑TGAAGGGGTTTtctt > 2:1330584/1‑75 (MQ=255) ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGC‑TAAATGGGTTTtttc > 2:1330596/1‑76 (MQ=255) ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGA‑TAAATGGGTTTtttc > 2:1438469/1‑76 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCg < 1:664173/76‑1 (MQ=255) tGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGt < 2:1125556/76‑1 (MQ=255) tGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGA‑TAAATGGTTTTTTTCggg > 2:1862451/1‑75 (MQ=255) gAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATTAAAGGGTTTTTTTCGGt > 2:620719/1‑76 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCg < 1:736872/76‑1 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCg < 2:1740225/76‑1 (MQ=255) aaCACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGG‑TAAATGGGTTTTTTCGGTTTCGc > 2:147701/1‑76 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGA‑TAAATGGGTTTTTTCGGTTCCGcc > 2:1593066/1‑75 (MQ=255) caccacCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCaa < 1:1438469/76‑1 (MQ=255) caccacCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGA‑TAAATGGTTTTTCTCGGGTCCGcca > 2:1054217/1‑74 (MQ=255) accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGG‑TAAATGGTTTTTTTCGGTTCCGCaaa > 1:1184507/1‑76 (MQ=255) accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGA‑TAAATGGGTTTTTTCGGTTTCGCaaa > 2:1164658/1‑76 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCg < 2:981260/76‑1 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATg < 2:206731/76‑1 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCA‑TAAATGGTTTTTTTCGGTTCCGCCAAGCGATg > 1:1425015/1‑76 (MQ=255) tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCG‑TAAATGGTTTTTCTCGGTTCCGCGAAGCGCTgcg > 1:374736/1‑76 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCAt < 2:192973/76‑1 (MQ=255) gtgCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATc < 1:1410487/74‑1 (MQ=255) agCTTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCt < 2:106211/76‑1 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCA‑TAAATGGTTTTTCTCCGTTCCGCAAAGCGATGCGCATCCGGt > 1:860127/1‑76 (MQ=255) cTGGAAGATGGCTCTTATCGGCTGGGGGGGCA‑TAAATGGGTTTTTTCGGTTCCGCCAAGCGATGCGCATCTGGTGc > 1:1418277/1‑76 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGtaat < 2:1609223/76‑4 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGGGGGGCA‑AAAATGGTTTTTCTCGGATCCGCAAAACAATGCGCATCTGGGGCt > 2:4399/1‑76 (MQ=255) gAAGATGGCTCTTATCGGCTGGGGGGGCA‑TAAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTgg > 2:1930143/1‑76 (MQ=255) aaGATGGCTCTTATCGGCTGGTGGGGGA‑TAAATGGTTTTTTTCGGTTCCGCAAAGCGATTCGGCTCTGGTGCTTgc > 2:467542/1‑76 (MQ=255) aaGATGGCTCTTATCGGCTGGTGGGGCA‑CAAATGGTTTTTCTCCGGGCCGCCAAGCCCCGCGCAACTGGCGCTTgc > 1:209459/1‑76 (MQ=255) gATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcgc < 2:36668/76‑1 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 2:1640620/76‑1 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 2:300491/76‑1 (MQ=255) ctctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc > 2:605870/1‑76 (MQ=255) ctctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc < 2:1423927/76‑1 (MQ=255) ctctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc < 2:616240/76‑1 (MQ=255) ctctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc < 2:1138262/76‑1 (MQ=255) ctctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGGGCTGGCGCAGAcc > 2:1795737/1‑76 (MQ=255) tctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcg < 1:1362101/76‑1 (MQ=255) tctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcg < 1:656423/76‑1 (MQ=255) ctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgc < 2:143231/76‑1 (MQ=255) ctTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgc < 2:1441934/76‑1 (MQ=255) tATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCgg < 2:949932/76‑1 (MQ=255) tATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCgg < 1:1916051/76‑1 (MQ=255) aTCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCggg < 1:1545404/76‑1 (MQ=255) aTCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCggg < 2:109204/76‑1 (MQ=255) cGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg < 1:1831033/76‑1 (MQ=255) ggCTGGCGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtgg > 1:1082253/1‑76 (MQ=255) gCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtggt > 1:206232/1‑76 (MQ=255) gCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtggt > 1:260452/1‑76 (MQ=255) cTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc < 2:982914/76‑1 (MQ=255) cTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc > 2:1210471/1‑76 (MQ=255) tggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 1:1331521/76‑1 (MQ=255) tggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:21637/76‑1 (MQ=255) tggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt > 2:1568417/1‑76 (MQ=255) tggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:152128/76‑1 (MQ=255) tggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 1:498339/76‑1 (MQ=255) ggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg < 2:190946/76‑1 (MQ=255) ggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg < 2:19087/76‑1 (MQ=255) ggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg > 2:1867030/1‑76 (MQ=255) ggtggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg < 2:1176845/76‑1 (MQ=255) tggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc < 1:1099313/76‑1 (MQ=255) tggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc > 1:1508860/1‑76 (MQ=255) tggGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc < 1:577133/76‑1 (MQ=255) | AAAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCA‑TAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392986‑4393136 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |