Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 1,343,489 | G→A | 100% | L542L (CTG→TTG) | yciR ← | modulator of Rnase II stability |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,343,489 | 0 | G | A | 91.7% | 63.4 / ‑1.3 | 24 | L542L (CTG→TTG) | yciR | modulator of Rnase II stability |
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (6/16); total (6/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGAAATATTTAC > REL606/1343342‑1343624 | gAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAAtt < 2:464073/150‑1 (MQ=255) aaTATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAGTtc < 1:418696/150‑1 (MQ=255) aaTATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATtc < 2:461241/150‑1 (MQ=255) tCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAActc < 1:441841/150‑1 (MQ=255) tttgtGGCATCCATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGGACCTAGTTGGCTAAATTGTTGAATCACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCGGTc < 1:210978/146‑1 (MQ=255) ttGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTc < 2:496311/150‑1 (MQ=255) cgatcgGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACAt < 1:299002/150‑1 (MQ=255) gttcgGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGAACCTAGTTGGCTAAATTGTTGAATAACGGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATc < 1:17542/148‑1 (MQ=255) gatcgGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAACAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATc < 2:638959/150‑1 (MQ=255) atcgGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCt > 1:554001/1‑150 (MQ=255) aaGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAggg < 2:31774/150‑1 (MQ=255) taagagacAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGc < 2:398067/143‑1 (MQ=255) agacactcAGATGCGAAAATTAAGAGAAACCGGTACCAAAATCAGACAGTTGTACTTGCGCACCTAGTTTGCTAAAATGTTGAATAACAGACGGTGGCAATTCATCATTCTCAATCAGACAAATCTCTGTCAGTTCAACATCTATAGGGc < 1:498471/143‑1 (MQ=255) aGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGc < 2:25447/150‑1 (MQ=255) aGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGc < 2:199524/150‑1 (MQ=255) aGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGc < 1:277119/150‑1 (MQ=255) ttGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGcctgtctctta > 1:396917/1‑140 (MQ=255) aGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGaaa < 1:61255/150‑1 (MQ=255) gTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATgg < 2:555608/150‑1 (MQ=255) ggcACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATc < 1:580628/149‑1 (MQ=255) gCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACGACTCCCTGTCAGTTCCACATCTATAGGGCAGTCTTCAACATTGAGCTCCTGGAGAACCTTTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACg > 1:112404/1‑150 (MQ=255) gttgAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGaaa > 2:346674/1‑150 (MQ=255) aaCAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGAAATATTTAc > 2:37780/1‑150 (MQ=255) aaCAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGAAATATTTAc > 2:3376/1‑150 (MQ=255) | GAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGAAATATTTAC > REL606/1343342‑1343624 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |