Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913_3_pae_tpiA | 4,171,489 | G→A | intergenic (+29/‑64) | rrlB → / → rrfB | 23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_pae_tpiA | 4,171,489 | 0 | G | A | 95.7% | 86.9 / ‑2.8 | 46 | intergenic (+29/‑64) | rrlB/rrfB | 23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon |
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (20/24); total (21/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 |
GATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTA > NC_000913_3_pae_tpiA/4171415‑4171630 | gATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa < 4:467891/150‑1 (MQ=11) aTGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg > 2:763187/1‑150 (MQ=14) aTGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg > 2:262678/1‑150 (MQ=14) aTGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg < 1:741257/150‑1 (MQ=14) tGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgcgc < 2:158876/150‑1 (MQ=17) gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg < 3:98205/150‑1 (MQ=32) gCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggt < 1:630212/150‑1 (MQ=17) cTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg > 1:582596/1‑150 (MQ=17) ccGGTACTGATGAAGCGGGAGGCTTAACCTTACAACGCCGATGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATGTGCCTGGCGGCAGTAGCGCGGTGGTcc < 2:63794/150‑1 (MQ=11) tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGa < 4:202667/150‑1 (MQ=17) tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGa > 1:428213/1‑150 (MQ=17) tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGa > 1:303612/1‑150 (MQ=17) tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGCAGCGCGGTGGTCCCACCTGa > 1:428412/1‑150 (MQ=17) aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAcc < 3:210784/150‑1 (MQ=17) tGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAccc < 3:387416/150‑1 (MQ=17) gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc < 4:354540/138‑1 (MQ=18) gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc > 3:354540/1‑138 (MQ=18) cGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGc < 4:271989/150‑1 (MQ=17) gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGcc > 1:6173/1‑150 (MQ=17) tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGTCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCg < 3:120245/150‑1 (MQ=17) tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAACCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCg > 1:490185/1‑150 (MQ=17) cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTc > 1:154042/1‑150 (MQ=17) cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTc > 1:154015/1‑150 (MQ=17) cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAGGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACTCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCCAACTc > 4:52578/1‑150 (MQ=12) ttAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCa > 1:535045/1‑150 (MQ=17) ttAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGCCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCa < 2:572413/150‑1 (MQ=17) aaCCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAAGTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGa > 2:518109/1‑150 (MQ=17) ccTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAg < 2:31385/150‑1 (MQ=17) cTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGt < 4:214696/150‑1 (MQ=17) cTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGACTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGt > 2:565575/1‑150 (MQ=12) ttACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTg < 1:763187/150‑1 (MQ=17) aCAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa > 2:306316/1‑150 (MQ=17) aCAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa > 3:16455/1‑150 (MQ=17) cAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa < 1:726974/150‑1 (MQ=17) aaCGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc > 1:116723/1‑150 (MQ=17) ccGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCg < 1:341417/150‑1 (MQ=17) cTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGTCCCCATGCCGAACTCAGAAGTGAAACGCCGTAgcgc < 1:777321/150‑1 (MQ=17) tGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGcc > 2:670089/1‑150 (MQ=17) ttttGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGa > 1:822477/1‑150 (MQ=17) tttGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGAt < 1:718697/150‑1 (MQ=17) ttGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCCGGCGGCAGTAGCGCGGTGGTCCCACCTGCCCCCATGCCTGACTCAGAAGTGAAACGCCGTAGCGCCGATg < 2:448190/150‑1 (MQ=11) ttGGCGGATGAAAGAAGATTTTCAGCATGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATg < 3:160080/150‑1 (MQ=17) ggCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGt > 2:461653/1‑150 (MQ=17) gCGGATGAAAGAAGATTTTGCGCCGGATACAGATTAAATCAGAACGCAGGAGCGGTCTGATAAAACAGAATTCGCCTGGCGGCAGTAGCGCGGTGGTCCCATCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTa < 2:428412/150‑1 (MQ=11) gCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTa < 2:428213/150‑1 (MQ=17) gCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTa < 1:262678/150‑1 (MQ=17) | GATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTA > NC_000913_3_pae_tpiA/4171415‑4171630 |
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: ‑ |