Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 2,393,232 | A→G | 19.7% | S89P (TCA→CCA) | moeA ← | Molybdopterin molybdenumtransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 2,393,232 | 0 | A | G | 19.7% | 53.2 / 6.7 | 33 | S89P (TCA→CCA) | moeA | Molybdopterin molybdenumtransferase |
Reads supporting (aligned to +/- strand): ref base A (12/13); new base G (3/3); total (15/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 = 7.14e-01 |
GGTTTACGAATTGTAAATGTATCTTCTAGTTCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTACAAT > USA300TCH1516_ALE/2393096‑2393367 | ggTTTACGAATTGTAAATGTATCTTCTAGTTCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCt < 1:23997/141‑1 (MQ=255) tGTAAATGTATCTTCTAGTTCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTCCACCAATATg > 1:190881/1‑141 (MQ=255) ttCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGGACCTGCACCCATACGATCAATCACTTTAAACTc > 2:154794/1‑141 (MQ=255) ttCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTc > 2:158085/1‑141 (MQ=255) ttCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCACCAACTAATTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTc > 2:216739/1‑141 (MQ=255) tGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACGAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTc < 1:87081/141‑1 (MQ=255) cAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTTGTGATCCCCAACTAATTTATCAGAAACTGAACCCCCCCCACTATGATAAATCACTTTAAACTAAAAGCGATTCtcaa > 2:96252/1‑138 (MQ=255) aacaGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGcccc < 2:238005/141‑1 (MQ=255) tCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACAAATTAATCCTGAACTGAACAAGCAACAATATGATCAATCACTTTAAACTCAATGCGATACTAACAACTTGCCCCTTGTaaa > 2:141017/1‑141 (MQ=255) tGCGCCATTCGGTATTTGTGCTACAATGGAAAAACGCACCGCTTCGTGATCCACAACTAATTTATCGGAAACTGAACCTGCACCAATATTATCAATCACTTTAATCTCAATGCGATTCTGACCACTTGCCCCTTGTGAATc < 1:56465/141‑1 (MQ=255) ttAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCACCTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTACACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTg > 2:85696/1‑141 (MQ=255) tttGTGCTCTAGTCATAATACGCACCGCTTCGTGATCCCCTACTAACTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTg < 2:213916/141‑1 (MQ=255) tttGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTg < 2:45353/141‑1 (MQ=255) tgtgCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCa > 1:268770/1‑141 (MQ=255) tgtgCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACACCACCCCCAATTCGATTCCGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCa > 2:267867/1‑141 (MQ=255) gtgCTCCAGTCATAAGACGCACCGCTTCGCGATCCCCAACTAATTTATCAGAAACTGTACCTGCACCAATATGATCAATCACTTTAAACTCCATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCaa < 2:260785/141‑1 (MQ=255) tgataCGGACAGCTTCGTGATTACCAACTAGGTTATCAGAAACTGCACTTGCACCAATATGATCAATCACTTGAAGCTCAGTGCGATTGTGGCGACTTGGCCCTTGTGTTTCAACACTGGGTATTGCAGAGCCataatatg < 2:255494/139‑3 (MQ=255) aataCGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAgg > 1:209004/1‑141 (MQ=255) ccGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTAt > 1:221697/1‑141 (MQ=255) tgCTTCTTGATACTCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATc < 2:135470/140‑1 (MQ=255) ttCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCaaa > 1:63765/1‑141 (MQ=255) ttCGTGATCCCCAACTAATTTATCAGAAACTGAACCCTCACCAATATGATCAATCACTTTAAACTCAATTCTTTTCTGACCACTTTCCNCTTGTGAATCAACACTCCGAATTGCAAAACCACCATAAGGCCCTTTCTCaaa > 2:75429/1‑141 (MQ=255) tCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGt < 1:269099/141‑1 (MQ=255) tcccAACTAATTTATCAGAAACTTAACCTACATCGATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACGTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGta < 2:203328/140‑1 (MQ=255) cAACTAATTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTatat < 1:158085/141‑1 (MQ=255) cAACTAATTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTatat < 1:216739/141‑1 (MQ=255) cAACTAAGTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTatat < 1:154794/141‑1 (MQ=255) aCTAATTCATCAGTCACTGAACCTGCACCAGTATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCa < 1:131470/141‑1 (MQ=255) tttATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGt > 2:198593/1‑141 (MQ=255) aTCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTACCAAACCTGGGTGTATCATAAGTAGc > 1:107652/1‑141 (MQ=255) aGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTAc < 1:157659/141‑1 (MQ=255) aGAAACTGAACCGGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTAc < 2:47865/141‑1 (MQ=255) aaCTGAACCTGCACCAATATGATCAATCCCTTTAAACTCAATGCGATTCTGACCACTTGCCCCCTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTACAat < 2:209822/141‑1 (MQ=255) | GGTTTACGAATTGTAAATGTATCTTCTAGTTCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTACAAT > USA300TCH1516_ALE/2393096‑2393367 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |