Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 2,372,557 | C→A | 18.4% | G407* (GGA→TGA) | pbuG ← | Guanine/hypoxanthine permease PbuG |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 2,372,557 | 0 | C | A | 18.4% | 122.9 / 16.1 | 62 | G407* (GGA→TGA) | pbuG | Guanine/hypoxanthine permease PbuG |
Reads supporting (aligned to +/- strand): ref base C (28/21); new base A (4/8); total (32/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.00e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01 |
CAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGG > USA300TCH1516_ALE/2372419‑2372694 | cAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCca > 1:598912/1‑140 (MQ=255) cAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAAGATGATGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCca > 1:224888/1‑140 (MQ=255) gttgTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCcaca < 2:30525/137‑1 (MQ=255) ttttttGTGTATAAAATTAACAATCATCCCATCCGAATTACAGAATGAATAGTACAAAAAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCa < 2:354138/140‑1 (MQ=255) gtCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATCCCAGt > 1:321949/1‑140 (MQ=255) gtCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGt < 2:528410/140‑1 (MQ=255) tCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGtt < 1:44762/140‑1 (MQ=255) tCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTATTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGAAATAAAAGtt > 1:65/1‑140 (MQ=255) cTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTg < 2:159212/140‑1 (MQ=255) ggATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCa < 1:460121/140‑1 (MQ=255) aaCCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGaa > 1:741161/1‑140 (MQ=255) aTGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATACGCATTGTAATCGGATAGAAAATAAATCCACACGCAATCCCAGTTGCAATTGAATATg > 2:46461/1‑140 (MQ=255) cAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATCAATCAACACGACATATCAGTTGCAATTGAATATGAGaaa > 2:475446/1‑140 (MQ=255) aaaTATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAg < 2:62490/140‑1 (MQ=255) aaaTATGAAGTACAGAATGAATAGTACCAATAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAg < 2:464484/140‑1 (MQ=255) aaTATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAgg < 2:414995/140‑1 (MQ=255) atatGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGc > 1:730846/1‑140 (MQ=255) gTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATaa < 1:190580/140‑1 (MQ=255) aGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCACTTGAATATGATAAAGGCATCATAATGAt > 1:683994/1‑140 (MQ=255) aatgaatAGTACCGTTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTg < 2:646422/140‑1 (MQ=255) aaagaatTGTACCATCAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGATATAAGCATTGTAATTGGATAGAAAATAAAGCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTg < 2:764733/137‑1 (MQ=255) gTACCATTAATTCATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAg < 2:511514/140‑1 (MQ=255) tACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGACATCATAATGATGGTAATGAAgg > 2:198385/1‑140 (MQ=255) cATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTg > 2:124169/1‑140 (MQ=255) tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAATGCTGGAACtc > 2:210898/1‑139 (MQ=255) tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTg > 1:519518/1‑140 (MQ=255) tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTg > 1:683455/1‑140 (MQ=255) tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAAGTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTg > 2:204499/1‑140 (MQ=255) agatTATCGGTTGTACTTCTTTATGTTTTTTAGATATAAGCATTGTAATTGGATAGAAAATAAATCTACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACt < 2:310644/138‑1 (MQ=255) acatTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACt < 1:356972/140‑1 (MQ=255) atTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTc < 1:547577/140‑1 (MQ=255) ggATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACt > 1:671216/1‑140 (MQ=255) ggATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTCTAATGAATGCCGGCACTTCAACTTCAAAACt > 2:198193/1‑140 (MQ=255) aTGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACgtt > 1:347983/1‑140 (MQ=255) tACTTCTTTATATTTTTGAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTcc < 2:311669/140‑1 (MQ=255) cttcttTATGTTTTTTAGAGATAAGCATTGTAATGGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAg < 2:598912/140‑1 (MQ=255) tcttTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGAACTGAAACTTCAAAACTTTTCCAGtt > 1:735722/1‑140 (MQ=255) cttTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTATTGAAGGCTGGCACTGCAACTTCTAAACTTTTCCAGttt > 2:277274/1‑140 (MQ=255) tATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATACACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATAATTGTAATGAAGGCTGGCAATGGAACTTCAAAAATTGTCCAgatgagt > 2:587403/1‑134 (MQ=255) ttttttAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAg > 1:302534/1‑140 (MQ=255) tttttAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGc > 1:441960/1‑140 (MQ=255) ttttAGAGATAAGCATTGTAATGGGATAGAAAATATATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCa < 2:689861/140‑1 (MQ=255) tatagaTAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg < 1:478135/137‑1 (MQ=255) tAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAACGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg < 2:741161/140‑1 (MQ=255) tAGAGATAAGCATAGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg < 1:216646/140‑1 (MQ=255) gagaTAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGtt < 2:671216/140‑1 (MQ=255) aGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGCTTGTAATGAAGGCTGGCACTGCCACTTCAAAACTCTACCAGCTTATCTCTGCAAAGTCTGCAGc > 1:711925/1‑140 (MQ=255) aGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGc < 1:501710/140‑1 (MQ=255) gCATTGTAATTGTATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGcc < 1:157758/140‑1 (MQ=255) gCATTGTAATTGGATAGAAATTGTTTGCGCACACAACACCAGATGCAGAAAAAGACTATAAAGGCATTATTGTGATGGGAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGcc < 1:181649/140‑1 (MQ=255) gCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGcc < 1:531799/140‑1 (MQ=255) cATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCa < 2:135866/140‑1 (MQ=255) ttGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCACAATGAGTGGAATGACGGCTGGCACTGCAACTTCAAAACTCTCCCAGCTTATCTCAGCAAAGTTCGCAGCCAtt > 2:359791/1‑140 (MQ=255) ttGTAATTGGATAGAAAATAAAGCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCAtt > 2:492470/1‑140 (MQ=255) gTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTaa > 2:222642/1‑140 (MQ=255) ttatgATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCACCTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGAAGCCATtactac > 2:536160/5‑140 (MQ=255) aTTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATac > 1:497077/1‑140 (MQ=255) ggATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACAcc < 1:188938/140‑1 (MQ=255) tAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACCCCTac > 2:617124/1‑140 (MQ=255) tAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTac < 1:176143/140‑1 (MQ=255) tAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTCCCGCCCTTAATCGACCTag > 2:153420/1‑139 (MQ=255) aGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGCTTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCTTCAATGGTTGCAGCCATTACTACACCTccc > 1:146408/1‑140 (MQ=255) aGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACCCCTacc > 1:392369/1‑140 (MQ=255) aaaaTAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCGGCCATTAATACACCTaccac > 1:42155/1‑140 (MQ=255) tCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAgg < 1:124169/140‑1 (MQ=255) | CAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGG > USA300TCH1516_ALE/2372419‑2372694 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |