Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,147,295 | 0 | G | C | 20.0% | 119.7 / 18.8 | 45 | R600P (CGC→CCC) | metG | methionyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base G (30/6); new base C (0/9); total (30/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.65e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTCAGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGGTCGTCACACCATTATG > REL606/2147195‑2147395 | aTTCAGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTgcg > 1:214730/1‑101 (MQ=255) cAGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCt < 2:1438832/101‑1 (MQ=255) aGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTg < 1:410910/101‑1 (MQ=255) ggAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGa < 2:1592372/101‑1 (MQ=255) cgatctaTCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGac < 2:1337369/95‑1 (MQ=255) aaaCCATCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGaca > 1:1331843/1‑101 (MQ=255) ccATCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTg > 1:1594527/1‑101 (MQ=255) aTCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGacagatcg > 1:1337375/1‑96 (MQ=255) aTCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGa > 2:732933/1‑101 (MQ=255) aTCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGa > 1:1198567/1‑101 (MQ=255) ccTTTGACTACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAGGTTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGAtctc < 2:1450999/101‑1 (MQ=255) ccTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGAtctc > 2:281757/1‑101 (MQ=255) cGCAAAAGTTGCCCTGCGCGTAGCGCTGATCAAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAC‑CTGCTGCCCCTGACACTGGATCTCGGCGGTGAAAAAc < 2:80814/101‑1 (MQ=255) caaaaGTTGCCCTGCGCGTAGCGCTGATCGAAACCGCAGAGTTTGTTGAAGGTTCTG‑ACAAC‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGc < 1:807744/99‑1 (MQ=255) cAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCa > 1:579602/1‑101 (MQ=255) aaaGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCaa > 2:980369/1‑101 (MQ=255) atgTTGCCCGGCGCGGAGCGCTGAACAAAAACGCAGATTTTGTTAAAGTTTCTG‑ACAAA‑CTGCTGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAAt < 1:1616858/99‑1 (MQ=255) aatttGACCGGCGCGTGGCGCTAATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑AAAAA‑CTCCGGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAAt < 2:775443/98‑1 (MQ=255) aaGTGGCCCGGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTA‑ACAAA‑CTGCTGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAAt < 1:1242403/101‑1 (MQ=255) aGTTGCCCTGCGCGTAGCGCTGATCAAAAACGCAGAGTTTGTGGAAGGTTCTG‑ACAAC‑CTGCTGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATg < 2:1012562/101‑1 (MQ=255) gTTGCCCTGCGGGTAGCCATGATCGAAAACGCAGAGTTTGTTGAAGTTTCTG‑AAAAA‑CTGCCGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGt < 2:494182/101‑1 (MQ=255) gTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGt > 1:1543881/1‑101 (MQ=255) tgcccTGCGGGGAGCGGTGATAAAAAACGCAGAGTTTGTGGAGGGTTCTG‑ACAAT‑CTGCCGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGtct < 2:538849/98‑1 (MQ=255) tgcccTGCGCGTAGCGCTGATCAAAAACGCAGATTTTGTTGAAGGTTCTA‑ACAAA‑CTGCTGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGtct < 2:473375/98‑1 (MQ=255) tGACCTGCGCGTAGCGCTGATCGAAAACGAAGAGGTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGCGGAAAACGCAATGtct > 2:883392/1‑101 (MQ=255) aCCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGtcttc > 1:41701/1‑101 (MQ=255) ccggcgcGAACCGCTAATCAAAACCGCAGAGTTTTTTGAAGGTTCTG‑ACAAC‑CTCCTGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGtcttct < 1:441676/98‑1 (MQ=255) cggcgcGTAGCGCTGATCGAAACCGAAGAGTTTGTTGAAGTTTCTG‑ACAAA‑CTCCGGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTtctc < 1:1261102/99‑1 (MQ=255) cTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTtctc > 2:904522/1‑101 (MQ=255) cTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTtctc > 1:221312/1‑101 (MQ=255) cTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTtctc > 1:1324145/1‑101 (MQ=255) cgTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTa > 1:1704308/1‑101 (MQ=255) gTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTAt > 1:609477/1‑101 (MQ=255) gcgcgaaacGAAAACGCAGATTTTGTTGAAGGTTCTA‑ACAAA‑CTGCTCCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCg < 2:1672796/93‑1 (MQ=255) gcTGATAAAAAACGCAGATTTTGTTGAAGTTTCTG‑ACAAC‑CTGCTCCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGtt < 2:677938/101‑1 (MQ=255) ccgaagagAAACCGAAGATTTTGTTAAAGGTTCTA‑ACAAA‑CTGCTGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGtt < 2:1688201/94‑1 (MQ=255) ttgataaaaacCGAAGATTTTGTTGAAGTTTTTG‑ACAAG‑CTGCTGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTc < 1:765409/94‑1 (MQ=255) cTGATGAAAAACGTAGATTTTGTTGAGGCTTGTG‑ACAAG‑CTGTCGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTc < 1:1644381/101‑1 (MQ=255) tgagagAAAACGCAGATTTTGTGGAAGTTTCTA‑ACAAA‑CTGCTGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCt < 2:1615985/96‑1 (MQ=255) tGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCt > 2:431541/1‑101 (MQ=255) gATCGAAACCGGAGATTTTGTTGAAGGTTCT‑‑ACAAACCTGCGGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTg < 1:1493332/101‑1 (MQ=255) gATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGCGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTg > 2:1120290/1‑101 (MQ=255) cGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCt > 1:1084526/1‑101 (MQ=255) gCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGa > 2:398294/1‑101 (MQ=255) cAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGAt > 1:28919/1‑101 (MQ=255) ttttttGAAGGTTCTA‑ACAAC‑CTGCCGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCa < 2:1025722/101‑1 (MQ=255) tgAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCAc > 2:1550909/1‑101 (MQ=255) tgAAGGTTCTG‑AAAAC‑CTGCTCCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCAc < 2:221312/101‑1 (MQ=255) aaGGTTCAGAACAAC‑CTGCGGCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACt < 2:213244/101‑1 (MQ=255) gTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGAtt > 2:1114113/1‑101 (MQ=255) gTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGAtt > 1:382837/1‑101 (MQ=255) cTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGgt > 1:1084474/1‑101 (MQ=255) cTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGgt > 1:109464/1‑101 (MQ=255) tG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGgtc > 1:40779/1‑101 (MQ=255) aCAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGgtcgt > 2:1554213/1‑101 (MQ=255) cAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGgtcgtc > 2:1226805/1‑101 (MQ=255) cAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGgtcgtc > 2:1436243/1‑101 (MQ=255) ctgctcCCCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGGTCGTCacac < 2:213103/101‑1 (MQ=255) gcCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGGTCGTCACACCATTAtg > 1:202721/1‑101 (MQ=255) | ATTCAGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTGCGCGTAGCGCTGATCGAAAACGCAGAGTTTGTTGAAGGTTCTG‑ACAAA‑CTGCTGCGCCTGACACTGGATCTCGGCGGTGAAAAACGCAATGTCTTCTCCGGTATTCGTTCTGCCTACCCGGATCCGCAGGCACTGATTGGTCGTCACACCATTATG > REL606/2147195‑2147395 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |