Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 2,775,879 | Δ1 bp | intergenic (‑49/+200) | insJ‑3 ← / ← cysH | IS150 hypothetical protein/phosphoadenosine phosphosulfate reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,775,878 | 0 | G | . | 98.5% | 322.2 / ‑0.5 | 66 | intergenic (‑48/+201) | insJ‑3/cysH | IS150 hypothetical protein/phosphoadenosine phosphosulfate reductase |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base . (35/30); total (36/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01 |
CGACTTCAAGGCGCTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCT > REL606/2775788‑2775973 | cGACTTCAAGGCGCTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGaaa > 2:1402339/1‑101 (MQ=18) aCTTCAAGGCGCTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAata > 1:1487266/1‑101 (MQ=255) aaGGCGCTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGcc < 2:239215/101‑1 (MQ=255) aaGGCGCTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGcc > 2:916019/1‑101 (MQ=255) aGGCGCTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCt > 2:49278/1‑101 (MQ=255) ggCGCTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTg < 2:1510289/101‑1 (MQ=255) gcgcTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GCTGATAGAAATATCGCCTgg < 2:1264364/101‑1 (MQ=255) tttCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCtt < 1:1187480/101‑1 (MQ=255) cAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTg > 1:427895/1‑101 (MQ=255) cAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTg > 2:1323014/1‑101 (MQ=255) cAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTg > 1:1614000/1‑101 (MQ=255) cAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTg > 1:769184/1‑101 (MQ=255) cAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTg > 1:802940/1‑101 (MQ=255) aaGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTcc > 1:70283/1‑101 (MQ=255) gggTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt > 2:1014242/1‑101 (MQ=255) gggTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt < 2:1296015/101‑1 (MQ=255) gTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt < 1:1024198/101‑1 (MQ=255) gTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt < 2:326464/101‑1 (MQ=255) gTATTTTGGCTTTGACATATTAGGGGCTATCCCATTTCAGCGTCCAACAAAATAGGTGCAGAACA‑GTTGACAGAAATAGCGCCGGGGCATTTTGTACATCt > 2:1269747/1‑101 (MQ=14) tATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg > 1:1254681/1‑101 (MQ=255) aTTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTgg < 1:171002/101‑1 (MQ=255) tttGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc < 1:1402345/101‑1 (MQ=255) tGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGa > 1:544213/1‑101 (MQ=255) tGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGa > 2:1052986/1‑101 (MQ=255) aCATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt < 1:968341/101‑1 (MQ=255) aCATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt < 2:265811/101‑1 (MQ=255) aCATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt < 1:965172/101‑1 (MQ=255) aCATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt < 1:649360/101‑1 (MQ=255) atatTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc > 2:775459/1‑101 (MQ=255) atatTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc > 1:102590/1‑101 (MQ=255) atatTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc > 2:1622524/1‑101 (MQ=255) tatTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt > 1:1026502/1‑101 (MQ=255) atTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCtt > 2:1035627/1‑101 (MQ=255) ttAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTa < 1:1269915/101‑1 (MQ=255) ttAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTa < 1:370781/101‑1 (MQ=255) ggggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc > 1:811464/1‑101 (MQ=255) ggggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACg < 2:975607/101‑1 (MQ=255) ggggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACg > 2:856663/1‑101 (MQ=255) ggggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACg < 1:360589/101‑1 (MQ=255) ggCTATTCCATTTCATCGTCGAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc < 1:1374933/101‑1 (MQ=255) ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc > 2:1410103/1‑101 (MQ=255) cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAg > 1:160024/1‑101 (MQ=255) cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAg > 2:1247453/1‑101 (MQ=255) cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAg < 1:1164423/101‑1 (MQ=255) tATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGa > 2:1118977/1‑101 (MQ=255) aTTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc > 1:1106619/1‑101 (MQ=255) tCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAg > 2:1155600/1‑101 (MQ=255) ccATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc < 1:1097893/101‑1 (MQ=255) aTTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt < 2:823356/101‑1 (MQ=255) aTTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt < 2:66131/101‑1 (MQ=255) aTTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt > 1:1537459/1‑101 (MQ=255) aTTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCAAGACAGCCt > 2:1644577/1‑101 (MQ=255) tttCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc < 1:1438390/101‑1 (MQ=255) tttCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc > 1:232337/1‑101 (MQ=255) tCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa < 1:1035632/101‑1 (MQ=255) cGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCAc > 2:115090/1‑101 (MQ=255) tCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCAccc < 1:1293162/101‑1 (MQ=255) caacaaAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGa > 2:220974/1‑101 (MQ=255) acaaAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACg < 2:317719/101‑1 (MQ=255) caaAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGc < 1:1622530/101‑1 (MQ=255) aaaaTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGcc > 2:133734/1‑101 (MQ=255) aaTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCAt < 2:540182/101‑1 (MQ=255) aTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCAtt > 2:1005682/1‑101 (MQ=255) ggTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCta > 1:1141745/1‑101 (MQ=255) ggTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCta < 2:1266378/101‑1 (MQ=255) gTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCtat < 1:1076662/101‑1 (MQ=255) cAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATAcc > 2:4481/1‑101 (MQ=255) aGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCt < 1:1277006/101‑1 (MQ=255) | CGACTTCAAGGCGCTTTTCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCT > REL606/2775788‑2775973 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |