Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 1,329,281 | C→G | intergenic (‑241/‑139) | yciN ← / → topA | hypothetical protein/DNA topoisomerase I |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,329,281 | 0 | C | G | 95.0% | 245.5 / 3.2 | 60 | intergenic (‑241/‑139) | yciN/topA | hypothetical protein/DNA topoisomerase I |
Reads supporting (aligned to +/- strand): ref base C (0/3); new base G (32/25); total (32/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.57e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01 |
TAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACT > REL606/1329209‑1329352 | tAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGAcc < 1:1576641/76‑1 (MQ=255) gATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTg > 1:1252601/1‑76 (MQ=255) tGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAAtt > 1:1491919/1‑76 (MQ=255) tGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAAtt > 1:795987/1‑76 (MQ=255) tGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTCAAtt > 2:825339/1‑76 (MQ=255) gCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAAttt < 2:155425/76‑1 (MQ=255) ttCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAAc < 1:481948/76‑1 (MQ=255) tCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACa > 2:172688/1‑76 (MQ=255) tCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACa < 2:125149/76‑1 (MQ=255) tCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACa < 1:1426101/76‑1 (MQ=255) tCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACa < 1:188863/76‑1 (MQ=255) cGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACag < 1:86380/76‑1 (MQ=255) gACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACaga < 2:1040582/76‑1 (MQ=255) aaaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCa < 1:1027510/76‑1 (MQ=255) aaaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCa > 2:1060514/1‑76 (MQ=255) aaaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCa < 1:867888/76‑1 (MQ=255) aaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCAc < 1:1060186/76‑1 (MQ=255) aTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGa < 1:737247/76‑1 (MQ=255) ttCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGAc < 1:909528/76‑1 (MQ=255) tCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACa < 1:1559979/76‑1 (MQ=255) gtGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGt > 2:223336/1‑76 (MQ=255) tGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGtt > 2:1311058/1‑76 (MQ=255) tGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGtt > 2:917017/1‑76 (MQ=255) gCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTg > 1:900106/1‑76 (MQ=255) ataGCGCCTGTCGGCAAAGACCTGTTAACTCAGAGACCTGAATTTTCGTGAACAGAGGCACGACAAGGGGTTGata > 2:1391138/1‑76 (MQ=255) ataGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGAGAACAGAGTCACGACAAGGGGTTGata > 2:610244/1‑76 (MQ=255) gcgcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATcc < 2:227640/76‑1 (MQ=255) cgcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCg < 2:1093982/76‑1 (MQ=255) cgcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCg > 2:1165621/1‑76 (MQ=255) gcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGc < 2:1435587/76‑1 (MQ=255) gcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGc > 1:471925/1‑76 (MQ=255) tGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCaga > 1:299307/1‑76 (MQ=255) tGTAGGCCAAGACCGGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCaga > 1:101968/1‑76 (MQ=255) gTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCagag > 2:905197/1‑76 (MQ=255) gTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCagag > 2:1670707/1‑76 (MQ=255) tAGGCCAAGACCTGTTAACTCAGTGACCTGTATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCagaga < 1:162138/76‑1 (MQ=255) gCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCg < 2:379588/76‑1 (MQ=255) cAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAg > 2:1362654/1‑76 (MQ=255) aaGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGt > 2:904969/1‑76 (MQ=255) aaGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGt > 1:705033/1‑76 (MQ=255) aaGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGt > 2:1545966/1‑76 (MQ=255) aaGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGt > 1:1541154/1‑76 (MQ=255) gACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTcc < 2:1283427/76‑1 (MQ=255) gACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTcc < 2:218363/76‑1 (MQ=255) aCCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCa < 2:652793/76‑1 (MQ=255) aCCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCa > 2:1022277/1‑76 (MQ=255) tGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCatat > 1:189867/1‑76 (MQ=255) tGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCatat < 1:1626277/76‑1 (MQ=255) ngTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCatat > 1:725883/2‑76 (MQ=255) gTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATc > 1:1130467/1‑76 (MQ=255) tAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCgg > 2:1042733/1‑76 (MQ=255) tAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCgg < 1:1085855/76‑1 (MQ=255) tAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCgg < 2:753146/76‑1 (MQ=255) tAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCgg > 1:578281/1‑76 (MQ=255) aaCTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGt > 1:467384/1‑76 (MQ=255) aCTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTa < 2:1531726/76‑1 (MQ=255) cTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTaa < 2:169481/76‑1 (MQ=255) cTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTaa < 2:541168/76‑1 (MQ=255) cAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACt > 1:654268/1‑76 (MQ=255) cAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACt > 2:982935/1‑76 (MQ=255) | TAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACT > REL606/1329209‑1329352 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |