Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 86,384 | C→A | intergenic (‑51/+42) | leuA ← / ← leuL | 2‑isopropylmalate synthase/leu operon leader peptide |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 86,384 | 0 | C | A | 91.3% | 68.9 / NA | 23 | intergenic (‑51/+42) | leuA/leuL | 2‑isopropylmalate synthase/leu operon leader peptide |
Reads supporting (aligned to +/- strand): ref base C (0/2); new base A (12/9); total (12/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACC > REL606/86311‑86455 | aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAacc > 2:84959/1‑76 (MQ=255) aaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACcg > 2:725475/1‑76 (MQ=255) gANTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACTGCGCANTg < 2:1095901/76‑1 (MQ=255) gACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATg > 1:805419/1‑76 (MQ=255) tGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGgcgc < 2:654531/76‑1 (MQ=255) gCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGtt < 2:1048835/76‑1 (MQ=255) cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGttt < 2:1255923/76‑1 (MQ=255) cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGttt > 2:413898/1‑76 (MQ=255) cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGttt < 2:261139/76‑1 (MQ=255) tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGgttt > 1:893685/1‑76 (MQ=255) gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGGTTtttgtt > 1:1055589/1‑76 (MQ=255) tGGGTCCTTCTCTTTTTTAGAGCCCCTCGCTGCGGGCAAAAAAAAACCCGCGCAAGGGCGCGGGTTTTTTGTTTGa < 1:867197/76‑1 (MQ=37) ggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACt > 2:563609/1‑76 (MQ=255) gTCCTTGTCTCTTTAAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAAGGGCGCGGGTTTTTTGTTTGACTg < 2:748751/76‑1 (MQ=255) tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGc > 2:1158807/1‑76 (MQ=255) ccTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAACCCGCCAAGGGGCGGGGTTTTTTTTTTTCTGCg > 1:1215499/1‑76 (MQ=255) tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAACACCGCGCCCTTGCGCGGGTGGTTTGTTTGTCTGCGTGc > 1:753191/1‑76 (MQ=255) gTCTCTTTTATAGCGCCCCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCt > 1:113511/1‑76 (MQ=39) tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGc > 2:557548/1‑76 (MQ=255) gcgcCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTGTTTTGTTTGACTGCGTGCTGGCTTAATGCTg < 1:1270169/76‑1 (MQ=255) ccTCGCTTCGTGCATAAAAAAAACCGCGCAATGGCGCGGGGTTTTTGTTTCACTACGTGGTGGCTTTATGCTGGAt > 2:598694/1‑76 (MQ=255) gCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCg > 2:648321/1‑76 (MQ=255) gCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTGTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCg < 2:258325/76‑1 (MQ=255) cTTCGGGCATAAGAACAACAGCGCAATGGTGCGTGTTTTTTGTATGAGTGCGTGCTGGCTTAATGCTGGATGCCGc > 1:1296211/1‑76 (MQ=255) cTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGc < 1:245700/76‑1 (MQ=255) ttCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTGTTTTGTTTGACTGTGTGCTGGCTTAATGCTGGATGCCGCt < 2:1263564/76‑1 (MQ=255) tCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCNCTc < 1:1201941/76‑1 (MQ=255) cATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCg > 2:451918/1‑76 (MQ=255) aTAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGt > 2:1209941/1‑76 (MQ=255) aaaaaaaaCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCt > 2:143977/1‑76 (MQ=255) aaaaagACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTa < 2:49454/76‑1 (MQ=255) aaagaCCCGCGCAGTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAc < 1:437915/76‑1 (MQ=255) aaaaACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGTTGGCTTAATGCTGGATGCCGCTCACTCGTCTAcc > 2:1279381/1‑76 (MQ=255) | AAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACC > REL606/86311‑86455 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |