Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,539,815 | A→T | 11.1% | intergenic (+109/‑348) | yhgF → / → feoA | putative transcriptional accessory protein/ferrous iron transporter, protein A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,539,815 | 0 | A | T | 11.1% | 143.7 / 9.1 | 45 | intergenic (+109/‑348) | yhgF/feoA | putative transcriptional accessory protein/ferrous iron transporter, protein A |
Reads supporting (aligned to +/- strand): ref base A (21/19); new base T (2/3); total (23/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 |
CAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAC > NC_000913/3539694‑3539958 | cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGATAT‑TCTCCGGc > 3:432748/1‑132 (MQ=255) cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGATAT‑TCTCCGGc < 4:432748/132‑1 (MQ=255) aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCa > 4:1432917/1‑149 (MQ=255) ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCata < 3:1432917/149‑1 (MQ=255) tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATAT‑TCTCCGGCActt > 2:262236/1‑108 (MQ=38) tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATAT‑TCTCCGGCActt < 1:262236/108‑1 (MQ=38) gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaat > 3:592441/1‑149 (MQ=255) aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatt > 1:1025970/1‑122 (MQ=255) aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatt < 2:1025970/122‑1 (MQ=255) aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa < 4:592441/149‑1 (MQ=255) aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAatata < 3:866898/99‑1 (MQ=255) aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAatata > 4:866898/1‑99 (MQ=255) ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata > 3:1097975/1‑109 (MQ=255) ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata < 4:1097975/109‑1 (MQ=255) cACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAg < 2:1343996/145‑1 (MQ=255) cACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAg > 1:1343996/1‑145 (MQ=255) gccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaatt > 1:579693/1‑133 (MQ=255) gccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaatt < 2:579693/133‑1 (MQ=255) cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAc > 1:373401/1‑124 (MQ=255) cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAc < 2:373401/124‑1 (MQ=255) cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa > 3:72224/1‑134 (MQ=255) cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa < 4:72224/134‑1 (MQ=255) tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCAttttt > 1:479627/1‑56 (MQ=255) tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCAttttt < 2:479627/56‑1 (MQ=255) tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAt < 1:561769/97‑1 (MQ=255) tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAt > 2:561769/1‑97 (MQ=255) tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa < 4:1310970/117‑1 (MQ=255) tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa > 3:1310970/1‑117 (MQ=255) tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAAtt > 3:1395075/1‑149 (MQ=255) cAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCAc < 4:1395075/149‑1 (MQ=255) aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCataa > 4:1223470/1‑73 (MQ=255) aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCataa < 3:1223470/73‑1 (MQ=255) aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc > 2:1084388/1‑131 (MQ=255) aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc < 1:1084388/131‑1 (MQ=255) acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAat > 3:1007352/1‑36 (MQ=255) acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAat < 4:1007352/36‑1 (MQ=255) acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa < 2:1237108/123‑1 (MQ=255) acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa > 1:1237108/1‑123 (MQ=255) aacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTc < 4:178594/149‑1 (MQ=255) gCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa > 2:493098/1‑117 (MQ=255) gCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa < 1:493098/117‑1 (MQ=255) gACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAtttt > 4:560889/1‑149 (MQ=255) aatGCCGGAGATAT‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa < 2:224090/98‑1 (MQ=39) aatGCCGGAG‑ATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa > 1:224090/3‑100 (MQ=39) cGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTa > 3:774971/1‑74 (MQ=255) cGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTa < 4:774971/74‑1 (MQ=255) ggAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAc > 1:851499/1‑149 (MQ=255) aaTA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCgg > 3:937732/1‑105 (MQ=255) aaTA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCgg < 4:937732/105‑1 (MQ=255) ata‑tCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTatat > 2:216284/1‑127 (MQ=255) ata‑tCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTatat < 1:216284/127‑1 (MQ=255) | CAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAC > NC_000913/3539694‑3539958 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |