Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 2,352,815 | G→A | 100% | A92V (GCA→GTA) | rpsI ← | 30S ribosomal protein S9 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 2,352,815 | 0 | G | A | 83.3% | 23.1 / 0.6 | 12 | A92V (GCA→GTA) | rpsI | 30S ribosomal protein S9 |
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (4/6); total (6/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAAAATTAATGATT > USA300TCH1516_ALE/2352670‑2352959 | tttGTATATATCGTCCGGCAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCtt > 1:266499/1‑149 (MQ=255) tttGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCtt > 1:176415/1‑149 (MQ=255) aCGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGt > 1:81548/1‑149 (MQ=255) gaATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATag < 1:283306/149‑1 (MQ=255) aGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGa < 1:99465/149‑1 (MQ=255) ttttttACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGt < 1:111450/149‑1 (MQ=255) aTACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATa > 1:90235/1‑149 (MQ=255) gTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTa > 1:226377/1‑149 (MQ=255) gTGGGTCACGAGTAAGTAATCCAGCGCGTTATAAAGAACCTCTGGATGCAGGATCTACTGCTAATAATGGACGAGCGAGTCCGTGACGGATAGCTTAAGCTAGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTa > 1:44493/1‑149 (MQ=255) ggTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCt < 1:269249/149‑1 (MQ=255) aaaGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATggttggtt < 1:282661/149‑1 (MQ=255) aTCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAAAATTAATGAtt < 1:144053/149‑1 (MQ=255) | TTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAAAATTAATGATT > USA300TCH1516_ALE/2352670‑2352959 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |