| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913_3_pae_tpiA | 1,197,502 | C→T | G11R (GGA→AGA) | ymfD ← | e14 prophage; putative SAM‑dependent methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_pae_tpiA | 1,197,502 | 0 | C | T | 97.1% | 102.5 / NA | 34 | G11R (GGA→AGA) | ymfD | e14 prophage; putative SAM‑dependent methyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (20/13); total (21/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCA > NC_000913_3_pae_tpiA/1197362‑1197642 | tAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat‑at > 1:47691/1‑149 (MQ=255)tAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCAAGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAAGATTATCTCTGTGAat‑at > 1:274794/1‑149 (MQ=255) ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCTTTTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAATTTTGTGCTt‑attat > 2:54814/1‑149 (MQ=255) ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAt‑attat > 1:637430/1‑149 (MQ=255) aGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCa > 2:153877/1‑149 (MQ=255) gAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGCTATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAAGATTAACTCTGGACGT‑ATTATAGTTCaa > 2:444315/1‑149 (MQ=255) atCTAAGTTTTACGCAGCCAAAATAAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGATGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCa < 2:152084/149‑1 (MQ=255) tAAGTTTACCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTACTAATATAATCGCAAAGGTGTCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGt > 1:241172/1‑149 (MQ=255) gcAAAATCAAGGGCGGCGCGATTTTTCTCAATGCCTCTAATATACTCGCAAAGATAACTAGAGGGCATCGTATGAGGTTTTGCTGCATTATCTGAGCGAATATTCA‑GCCTTTCATCCTTTCAGGTCAAAGCAAGTACCATTCAATTAAt > 2:319803/2‑149 (MQ=255) cAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATg > 1:499353/1‑149 (MQ=255) cAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAttt > 1:57744/1‑149 (MQ=255) aaGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAtttt < 1:153877/149‑1 (MQ=255) gtgTGGCCTTTTTTCTCAATTCTTCGAATAGATTCGCAAAGAGATCTAGAGGGTATCGTATGAGGTTTTGCTGCATTCTCTTAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGt < 2:198418/147‑1 (MQ=255) caaTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTCTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaaaaa < 1:140592/147‑1 (MQ=255) aTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAAtt > 2:7787/1‑149 (MQ=255) tttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt < 2:499353/149‑1 (MQ=255) aTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttt > 1:658582/1‑149 (MQ=255) cTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTc < 2:179071/149‑1 (MQ=255) tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTATATCGTcc > 2:68568/1‑149 (MQ=255) ataCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATCTTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGa > 2:204462/1‑149 (MQ=255) tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATACTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTCTCGTCCGGGAtttt > 2:187894/1‑149 (MQ=255) gATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGTTTTTTCCgcct > 2:131503/1‑149 (MQ=255) ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAAc < 1:121004/149‑1 (MQ=255) gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa < 1:2487/147‑1 (MQ=255) gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa > 2:2487/1‑147 (MQ=255) gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACt < 1:575751/149‑1 (MQ=255) gtaaaatctcTGAGCAAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACGATTCAATTAATGTTATTTTGAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaa < 1:67002/143‑1 (MQ=255) tgcATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTCTCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag > 1:638659/1‑149 (MQ=255) gAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataat > 1:224478/1‑149 (MQ=255) aGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataata > 2:106497/1‑149 (MQ=255) gCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATc > 2:21317/1‑87 (MQ=255) gCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATc < 1:21317/87‑1 (MQ=255) gCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAc < 1:7787/149‑1 (MQ=255) actAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGAGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca < 1:346466/147‑1 (MQ=255) | TAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCA > NC_000913_3_pae_tpiA/1197362‑1197642 |
| 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: ‑ |