Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 2,758,314 | A→C | 100% | K374Q (AAA→CAA) | ygbK → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,758,314 | 0 | A | C | 93.6% | 153.4 / 1.5 | 47 | K374Q (AAA→CAA) | ygbK | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (1/2); new base C (23/21); total (24/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.09e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CAGCAATACGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATAAAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAG > REL606/2758170‑2758458 | cAGCACTACGGTGCACAAACAGCCAGGCAGGCAGTAGAAACACTGTCTTCTCAACTAGCGGCGCGGTTAGCAGCGGACGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGGCCTCCGGCGGAGGCACACAGAGCCTGGGAATACAAggg < 1:332503/150‑1 (MQ=255) cAGCAATACGGTGCACAAACAGCCAGTCAGGTAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGATAGCAGCGGAAGGCGTGATACGCTTTATTGTAGCAGGCGGTGAGACCTCAGGCGTAGTCACACAGAGCCTGGGAATACAAggg > 1:636619/1‑150 (MQ=255) cAGCAATACGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAggg > 2:226450/1‑150 (MQ=255) cAGCAATACGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAAAACAAggg > 2:313582/1‑150 (MQ=255) gCAATACGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGtt > 2:41606/1‑150 (MQ=255) gCAATACGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGtt > 1:277187/1‑150 (MQ=255) aCGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCata < 2:138047/150‑1 (MQ=255) aCGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCata < 2:603982/150‑1 (MQ=255) aCGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATAAAGGGGTTTCata < 2:465654/150‑1 (MQ=255) caAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCa > 2:420751/1‑150 (MQ=255) gCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGc > 2:175064/1‑150 (MQ=255) gaaGAACCACTGTTTTCTCAACTATCCCCGCGGTTAGCAGCGGAAGGCGGGACACGCCTTAATGTCGCAGGCGGAGCGACCTCCGGCGTAGTCAAACACAGTCGGGGAATACCAGGGATTCATATTTGCCCAACCGTTTCCCCCGGCGTg > 1:96494/3‑150 (MQ=255) ttCGCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGc < 1:136016/150‑1 (MQ=255) gcggcgCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAg < 2:374218/150‑1 (MQ=255) cggcgCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTGCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑c < 1:247336/150‑1 (MQ=255) ggcgCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGCGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑cc < 1:422335/150‑1 (MQ=255) ggcgCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑cc < 2:406201/150‑1 (MQ=255) ggcgCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑cc < 2:447326/150‑1 (MQ=255) ggcgCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑cc < 1:121450/150‑1 (MQ=255) ggcgCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑cc < 1:419555/150‑1 (MQ=255) ggcgCGGTTAGCAGCGGAAGGCGTAACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑cc < 1:321020/150‑1 (MQ=255) ggcgCAGTCAGCTGCGGACAGCGTGACACCCCTTATTGTCGCAGGCGGTGCGACCTCCGGCGTAGTCACACAGGGCCGGGGAAGACAAGGGTTTCATGTTGGCCCAACCATGTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑cc < 1:511841/150‑1 (MQ=255) gcgcGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCt > 2:418495/1‑150 (MQ=255) cgGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGtc < 2:638487/150‑1 (MQ=255) ggTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGGGCCTGGGGATACAAGGGTTTCATATTGTCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGAACGGGATAAG‑CCTGtct > 1:673924/1‑150 (MQ=255) ggTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATAAAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGtct < 2:536981/150‑1 (MQ=255) gacaggcagcGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGtct < 2:402033/145‑1 (MQ=255) gTTAGCAGCGGAAGGCGGGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCAAACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGtctc < 1:229566/150‑1 (MQ=255) tAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCAc > 1:104830/1‑150 (MQ=255) gcagcGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGtctctgtc > 1:400875/1‑146 (MQ=255) gCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAGCCCTGtctcttatacacatctc > 2:9777/1‑137 (MQ=255) cGTGAAACGCGTTATTGTAGCGGGCGGTGAGTCCGCCGGCGTAGTCCCACAGACCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATc < 1:204844/150‑1 (MQ=255) gtagCAGGCGGTGAGACATCCGACGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCACGACCCAGTCCCCCGGCGTGCAGTGGGAAAACGCACTGGATAAT‑CCTGTCTAACTCGCCCTTAAATCTGGCAAATTCGGtgtt > 1:178183/4‑148 (MQ=255) ggCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGcc > 2:61538/1‑150 (MQ=255) ggagagaCCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTATCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAGCGCACTGGATACC‑CGTGTGTCACTCGCCCTTAACTCTGGCCACGTCGgtggtggatcattt > 1:438978/4‑142 (MQ=255) ggTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCttt > 2:100240/1‑150 (MQ=255) cTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACg > 2:595740/1‑150 (MQ=255) aGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAag < 2:105135/150‑1 (MQ=255) gTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAaga > 2:470752/1‑150 (MQ=255) tCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAagag < 1:99944/150‑1 (MQ=255) cacaGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGCGt > 1:526018/1‑150 (MQ=255) cacaGAGCCTGGGAATAAAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAagagaag > 1:41358/1‑148 (MQ=255) gagCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACCAGCCCAGCGAGAGttttt > 1:630926/1‑150 (MQ=255) gagCCTGGGAATACAAGGGTTTCATATTCGCCCAACCACTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCCCCCTTCAATCTGGCCACTTCGTGCATGACGCCTTTTTTTCCCGCGCCCAAAGAGAGtgctt > 1:216088/1‑146 (MQ=255) cagggAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCCCTGGATAAG‑CCAGTCTCACTCGCCCTTAAATCTGGACACTGCGGTGATGACGCCTTTTTTTCACGAGACCAAAGAGTGTTTTTATCa < 1:208529/148‑1 (MQ=255) cTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCa < 2:527533/150‑1 (MQ=255) cTGGGAATAAGAGGGTCTCATACTGGCCCAACCATTTCCCCCGGCGTGCCGATGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAAGCTGGCAACTTCGGTGATGTCGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCa < 1:127979/150‑1 (MQ=255) ggAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGa > 1:293566/1‑150 (MQ=255) gAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGCGCCCAAAGAGGGTTTTTATCATGac > 1:213620/1‑149 (MQ=255) gAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAg > 2:65828/1‑150 (MQ=255) | CAGCAATACGGTGCACAAAAAGCCAGTCAGGCAGTAGAAACACTGTTTTCTCAACTAGCGGCGCGGTTAGCAGCGGAAGGCGTGACACGCTTTATTGTCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATAAAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTAAACGCACTGGATAAG‑CCTGTCTCACTCGCCCTTAAATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAG > REL606/2758170‑2758458 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |