Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,100,865 | 0 | A | T | 25.0% | 92.4 / 11.3 | 44 | Y402F (TAT→TTT) | ycdB | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (10/23); new base T (11/0); total (21/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.60e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGTTTTT > REL606/1100790‑1100939 | aaGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTa < 2:1963865/76‑1 (MQ=255) aaGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTa < 2:1685844/76‑1 (MQ=255) cTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTAtttt > 1:96400/1‑76 (MQ=255) ttCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTAttttt < 2:1589185/76‑1 (MQ=255) ttCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTAttttt > 2:1085284/1‑76 (MQ=255) tCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGTATACGTTAAACCTATCGGCGGA‑GGTTCtttttt > 2:1286593/1‑76 (MQ=255) ccTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTAttttttt > 1:289636/1‑76 (MQ=255) cTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTttttttttt > 1:1209879/1‑75 (MQ=255) tGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTgc < 2:1659588/76‑1 (MQ=255) cAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTTTTTTTTTGCGCt > 2:1515516/1‑76 (MQ=255) aaaaaaGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTTTTTTTTTGTGCTGCCggg > 1:262587/1‑76 (MQ=255) aaaGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGt < 1:97727/76‑1 (MQ=255) aaaGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGt < 2:157908/76‑1 (MQ=255) aaaGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTTTTTTTTTTGGGTGTCgggggg > 2:233094/1‑75 (MQ=255) aaGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTg < 2:569092/76‑1 (MQ=255) aaGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTTTTTTTTTGCGCTGCCggggggg > 2:27373/1‑74 (MQ=255) aaGGCACAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTg < 1:1076342/76‑1 (MQ=255) gCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAg < 1:1402560/76‑1 (MQ=255) tCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTTTTTTTTTGTGCCGCCGGGGGGGGAggg > 2:738566/1‑75 (MQ=255) cAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTTTTTTTTTTTGGTGGCGGGGGGGGAgggg > 1:1383549/1‑74 (MQ=255) tGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcgc‑g < 1:1201264/76‑1 (MQ=255) tGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTATTTTTTTGCGCTGCCGGGGGGGAAGGAggc‑g > 2:490849/1‑76 (MQ=255) gggCGAAGCGCTGGAGGAATACGTTAACCCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcgc‑g < 1:1285248/75‑1 (MQ=255) gAAGCGCTGGAGGAATACGTTAAAACTATCGGCGGC‑GGTTTTTTTTTTTTCCTGCCGGGGGTGGAGGGCGC‑GAACg > 2:1671622/1‑76 (MQ=255) aaGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGGGGAGGGCGC‑GAACGa > 2:1851226/1‑76 (MQ=255) aaGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTTTTTTTTTGCGCTGCCCGGGGGGGAGGGCGC‑GGACGa > 2:1449381/1‑76 (MQ=255) aaGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGGTATTTTTTTGCGCTGCCGGGGGGGGAGGACGC‑GAACgg > 2:212015/1‑75 (MQ=255) cgcTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGatta < 1:254630/76‑1 (MQ=255) gcTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGattat < 1:149821/76‑1 (MQ=255) tGGAGGAATACGTTAAACCTATCGGCGGCGGGTTTTTTTTTTTGGCGCCCGGGGGGGAGGGGGGC‑G‑ACGGTTAttt > 1:1173899/1‑76 (MQ=255) tGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAATGATTAttt < 2:709765/76‑1 (MQ=255) ggaATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGa < 2:1062904/76‑1 (MQ=255) cGTTAAACCTATCGGCGGC‑GGTTTTTTTTTTTCCCCGCCGGGGGGGAAGGAAGCGGAACGTTATTTTCGCAAggcg > 2:1492831/1‑73 (MQ=255) cGTTAAACCTATCGGCGGC‑GGGTTTTTTTTTGCGCTGCCGGGGGGGGAGGGCGC‑GAACGATTATTTTGGGAGcgcg > 2:838333/1‑76 (MQ=255) ttAACCCTATCGGCGGG‑GGTTA‑TTTTTTGCGCGGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTa < 1:860415/76‑1 (MQ=255) tAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACTATTATTTCGGAAGCGCGTTa < 1:1323260/76‑1 (MQ=255) tAAACCTATCGGCGGC‑GGGTTTTTTTTTTGGCTGCCGGGGGGGGAGGGGGC‑GGAAGATTTTTTTGGAAAcgcggtt > 1:1433132/1‑73 (MQ=255) tAAACCTAGCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTa < 1:733648/76‑1 (MQ=255) aaaCCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGGGAAGGGCGC‑GAACGATTATTTCGGAAGCGCGGTAt > 2:48302/1‑76 (MQ=255) agccTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTAtt < 1:1268987/74‑1 (MQ=255) aaCCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGGGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTAtt > 2:1795650/1‑76 (MQ=255) aCCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTg < 2:321651/76‑1 (MQ=255) ccTATCGGCGGC‑GGGTTTTTTTTTTCGCTGCCGGGGGGGGAGGGCGC‑GGACGATTATTTTGGGAGCGCGGTTTTGc > 2:1770037/1‑76 (MQ=255) tATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCgg < 2:518230/76‑1 (MQ=255) cgggggc‑ggTTATTTTTTTGCGCGGCGGGGGGTGAAGGACGC‑GAACGA‑TATTTCGGAAGCGCGTTATTGCGGGttt < 1:1792578/76‑1 (MQ=255) ggcggc‑ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGttt < 1:1887003/76‑1 (MQ=255) gcggc‑ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTa < 1:1580849/76‑1 (MQ=255) gcggc‑ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTa < 2:1739743/76‑1 (MQ=255) gggc‑ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTaa < 1:1934303/75‑1 (MQ=255) cggc‑ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTaa < 2:1991371/76‑1 (MQ=255) ggc‑ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAAt < 1:744047/76‑1 (MQ=255) ggc‑ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAAt < 1:749309/76‑1 (MQ=255) ggc‑ggTTAGTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAAt < 2:160758/76‑1 (MQ=255) gc‑ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATg > 2:999479/1‑76 (MQ=255) gtggTTATTTTTTTGCGCTGCGGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATg < 1:1509644/74‑1 (MQ=255) ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGtt < 2:1345110/76‑1 (MQ=255) ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGtt > 1:203997/1‑76 (MQ=255) ggTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGtt > 1:1225372/1‑76 (MQ=255) gTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGttt < 2:1214426/76‑1 (MQ=255) gTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGAGTTAATGttt > 2:771603/1‑76 (MQ=255) ttATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGtttt > 1:1771066/1‑76 (MQ=255) ttATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGtttt > 1:1284987/1‑76 (MQ=255) tATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGttttt < 2:1596335/76‑1 (MQ=255) | AAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGC‑GGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGC‑GAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGTTTTT > REL606/1100790‑1100939 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |