Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 169,748 | T→A | 12.2% | intergenic (+21/‑30) | fhuA → / → fhuC | ferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 169,748 | 0 | T | A | 12.2% | 149.7 / 11.2 | 49 | intergenic (+21/‑30) | fhuA/fhuC | ferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase |
Reads supporting (aligned to +/- strand): ref base T (23/20); new base A (4/2); total (27/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01 |
GCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGT > NC_000913/169600‑169893 | gctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAt < 2:429097/149‑1 (MQ=255) gTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCTAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGc > 3:284223/1‑149 (MQ=255) gTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt < 2:476210/135‑1 (MQ=255) gTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt > 1:476210/1‑135 (MQ=255) ttAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt < 4:813126/134‑1 (MQ=255) ttAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt > 3:813126/1‑134 (MQ=255) ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACaa < 3:1341138/136‑1 (MQ=255) ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACaa > 4:1341138/1‑136 (MQ=255) aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc < 3:331771/117‑1 (MQ=255) aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc > 4:331771/1‑117 (MQ=255) gctaGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCAt < 4:284223/146‑1 (MQ=255) gctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGc < 2:1224353/100‑1 (MQ=255) gctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGc > 1:1224353/1‑100 (MQ=255) acacTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAttt < 2:22114/100‑1 (MQ=255) acacTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAttt > 1:22114/1‑100 (MQ=255) aTGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAg < 1:922749/115‑1 (MQ=255) aTGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAg > 2:922749/1‑115 (MQ=255) aaCGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTcc < 3:331426/112‑1 (MQ=255) aaCGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTcc > 4:331426/1‑112 (MQ=255) aGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGTATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCgtgc < 4:876640/144‑2 (MQ=255) aGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCgtgc > 3:876640/1‑143 (MQ=255) gTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGTCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCggg > 1:502291/1‑149 (MQ=255) ttgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACCAGTTGGCt < 2:870631/65‑1 (MQ=255) ttgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCt > 1:870631/3‑67 (MQ=255) aaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATac < 4:709323/78‑1 (MQ=255) aaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATac > 3:709323/1‑78 (MQ=255) cTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATcc < 2:502291/149‑1 (MQ=255) ttCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATc > 1:1043665/1‑147 (MQ=255) ttCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATc < 2:1043665/147‑1 (MQ=255) tttCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCg < 1:1332697/144‑1 (MQ=255) tttCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCg > 2:1332697/1‑144 (MQ=255) ttCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGc < 3:734830/45‑1 (MQ=255) ttCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGc > 4:734830/1‑45 (MQ=255) ttttGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTc > 1:515276/1‑67 (MQ=255) ttttGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTc < 2:515276/67‑1 (MQ=255) gCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCg < 2:35257/146‑1 (MQ=255) gCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCg > 1:35257/1‑146 (MQ=255) cACGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGaa > 1:411857/1‑149 (MQ=255) atggATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt > 1:686680/3‑115 (MQ=255) aCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt > 2:304119/1‑115 (MQ=255) aCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt < 2:686680/115‑1 (MQ=255) aCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt < 1:304119/115‑1 (MQ=255) aCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGaaa > 3:987088/1‑149 (MQ=255) aCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGCGCCCGGGCACACgaaa > 4:1402859/1‑110 (MQ=255) aCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGCGCCCGGGCACACgaaa < 3:1402859/113‑4 (MQ=255) aCGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATGCACGa < 4:1411056/51‑1 (MQ=255) aCGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATGCACGa > 3:1411056/1‑51 (MQ=255) aCGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTa > 2:678767/1‑81 (MQ=255) ggATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt > 3:866056/1‑149 (MQ=255) ggATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt > 3:237868/1‑149 (MQ=255) | GCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGT > NC_000913/169600‑169893 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |