Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,161,106 | A→G | E84G (GAA→GGA) | rpoB → | DNA‑directed RNA polymerase subunit beta |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,161,106 | 0 | A | G | 94.5% | 200.7 / 0.1 | 55 | E84G (GAA→GGA) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (29/23); total (30/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.85e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 |
GCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCG > REL606/4161037‑4161177 | gCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCa < 1:379887/76‑1 (MQ=255) aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAAt < 2:1353664/76‑1 (MQ=255) ggTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTg > 2:1257684/1‑76 (MQ=255) ggTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGGCAAATCCggg > 1:302385/1‑74 (MQ=255) aaTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCg > 1:960977/1‑76 (MQ=255) aaTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCg > 1:2523/1‑76 (MQ=255) aTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTGCAGGGATGTCAAATCCGTGGCGt > 2:808053/1‑76 (MQ=255) ccGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGAc > 1:675480/1‑76 (MQ=255) aGCTGCAATACGTCAGCTTCCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTa < 2:1332394/76‑1 (MQ=255) cTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTAtt > 1:237089/1‑76 (MQ=255) tGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTc > 2:199146/1‑76 (MQ=255) tGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTc > 2:553050/1‑76 (MQ=255) aTACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTTCGCa > 2:952124/1‑76 (MQ=255) tACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCAc < 2:1563589/76‑1 (MQ=255) tACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCAc > 2:329197/1‑76 (MQ=255) tACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCAc > 2:460004/1‑76 (MQ=255) tACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCAc > 2:988939/1‑76 (MQ=255) tACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCAc > 1:478432/1‑76 (MQ=255) tACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGNGACCTATTCCGCAc < 1:868279/76‑1 (MQ=255) cGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCTAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCg < 1:199110/76‑1 (MQ=255) cGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCg > 1:930238/1‑76 (MQ=255) gTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGc > 1:128687/1‑76 (MQ=255) cAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTg > 1:1586436/1‑76 (MQ=255) cAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTg < 1:524542/76‑1 (MQ=255) gCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGGGACCTTTTCCGCACCGCTgcg > 2:271903/1‑76 (MQ=255) cTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTgcgc < 2:615982/76‑1 (MQ=255) tACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCTAATCCGTGGCGTGACCTATTCCGCACCGCTGcgcg < 2:1448635/76‑1 (MQ=255) cGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTa < 2:1411745/76‑1 (MQ=255) cGCCTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTa < 2:685225/76‑1 (MQ=255) gCCTTGGCGAACCGGTGTTTGACGTCCAGTGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa < 2:466218/76‑1 (MQ=255) ccTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaaa < 1:240071/76‑1 (MQ=255) ccTTGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaaa < 2:52298/76‑1 (MQ=255) tGGCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTg < 1:461195/76‑1 (MQ=255) ggCGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGc > 1:1328243/1‑76 (MQ=255) cGAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGt < 1:69352/76‑1 (MQ=255) gAACCGGTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTc < 2:1504397/76‑1 (MQ=255) cGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGt > 1:1551580/1‑76 (MQ=255) ggTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg < 2:65956/76‑1 (MQ=255) ggTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg > 1:508811/1‑76 (MQ=255) ggTGTTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg < 1:574235/76‑1 (MQ=255) tgtTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGAt > 1:1105137/1‑76 (MQ=255) gtTTGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc < 2:70361/76‑1 (MQ=255) tttGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCt > 1:131148/1‑76 (MQ=255) tttGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCt < 2:960977/76‑1 (MQ=255) ttGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTa > 2:627835/1‑76 (MQ=255) ttGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTa > 1:124594/1‑76 (MQ=255) tGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTAt > 2:42270/1‑76 (MQ=255) tGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTAt > 2:914623/1‑76 (MQ=255) tGACGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGGGATCTAt > 2:625410/1‑76 (MQ=255) aCGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGTGTTAAACTGCGTCTGGTGATCTATGa < 2:783023/76‑1 (MQ=255) aCGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGa < 2:704044/76‑1 (MQ=255) cGTCCAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAg < 2:1076197/76‑1 (MQ=255) ccAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAgcgc > 2:542433/1‑76 (MQ=255) ccAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAgcgc > 1:912549/1‑76 (MQ=255) cAGGGATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGcgcg < 1:1196666/76‑1 (MQ=255) | GCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCG > REL606/4161037‑4161177 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |