Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AP009048 | 4,371,271 | Δ2 bp | 100% | intergenic (‑6/+299) | dcuA ← / ← aspA | C4‑dicarboxylate antiporter/aspartate ammonia‑lyase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AP009048 | 4,371,271 | 0 | A | . | 98.2% | 224.4 / ‑1.2 | 56 | intergenic (‑6/+300) | dcuA/aspA | C4‑dicarboxylate antiporter/aspartate ammonia‑lyase |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base . (25/30); total (25/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
* | AP009048 | 4,371,272 | 0 | A | . | 98.2% | 222.4 / ‑1.2 | 56 | intergenic (‑7/+299) | dcuA/aspA | C4‑dicarboxylate antiporter/aspartate ammonia‑lyase |
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (25/30); minor base G (0/1); total (25/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTA > AP009048/4371203‑4371334 || gTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa < 1:5010206/71‑1 (MQ=255) ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat < 1:2965827/71‑2 (MQ=255) ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat < 1:2430263/71‑2 (MQ=255) ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa > 1:4218641/1‑70 (MQ=255) cGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa < 1:2425780/67‑1 (MQ=255) cagcagCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaa < 1:485685/61‑1 (MQ=255) cagcagCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg < 1:3162383/70‑1 (MQ=255) agcagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGa < 1:1500726/70‑1 (MQ=255) agcagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt < 1:222372/71‑1 (MQ=255) agcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAc < 1:703562/71‑1 (MQ=255) gcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCCACCTGCATGGCAG‑‑AAATAGCGGCGATAAcc < 1:2830235/71‑1 (MQ=255) cTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCg > 1:4917176/1‑71 (MQ=255) ttttCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGAGGCGATAACCGc > 1:2510783/1‑71 (MQ=255) tttCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGc > 1:2001963/1‑70 (MQ=255) tctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGAt > 1:4896991/1‑70 (MQ=255) tctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATa < 1:1752987/71‑1 (MQ=255) ctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa > 1:2076693/1‑71 (MQ=255) tgatgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa > 1:83617/1‑70 (MQ=255) atgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATg > 1:3757946/1‑70 (MQ=255) atgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCAATGATAATg > 1:3280303/1‑70 (MQ=255) tgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGga < 1:301849/71‑1 (MQ=255) tgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCCATAACCGCCATGATAATGga < 1:3694452/71‑1 (MQ=255) gaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGga > 1:237792/1‑70 (MQ=255) aaCCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg < 1:5359522/49‑1 (MQ=255) aaCCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGgag < 1:930103/70‑1 (MQ=255) aCCAGATAGTCCAGACCGCCAGCAACCTGCATGG‑‑GAGAAATAGCGGCGATAACCGCCATGATAATGGAGAt < 1:1596109/71‑1 (MQ=255) cAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGa < 1:569851/36‑1 (MQ=255) cAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGa < 1:2608270/71‑1 (MQ=255) cAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGa < 1:2224253/71‑1 (MQ=255) aGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaa > 1:3033125/1‑36 (MQ=255) aGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGa > 1:1445194/1‑70 (MQ=255) gATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGa > 1:478540/1‑57 (MQ=255) aTAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa < 1:4957307/59‑1 (MQ=255) aTAGTCCAGACAGCCAGCAACCTGCATGGCAGaaa > 1:4436122/1‑35 (MQ=255) aGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGa > 1:4157694/1‑54 (MQ=255) aGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa > 1:4434280/1‑57 (MQ=255) tCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCg < 1:2624536/69‑1 (MQ=255) cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAAc < 1:767584/70‑1 (MQ=255) aGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGa < 1:4028718/67‑1 (MQ=255) gACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg > 1:4161092/1‑70 (MQ=255) gACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg > 1:1030668/1‑70 (MQ=255) aCCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGa > 1:199372/1‑47 (MQ=38) ccGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt > 1:64943/1‑56 (MQ=255) cGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt < 1:2575728/55‑1 (MQ=38) gCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGt < 1:4641476/71‑1 (MQ=255) ccAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGt > 1:2604381/1‑70 (MQ=255) cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt < 1:5344189/52‑1 (MQ=38) cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt < 1:3090219/52‑1 (MQ=38) cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt < 1:5078904/52‑1 (MQ=38) cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGtt < 1:173872/70‑1 (MQ=255) aGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACggg > 1:4197969/1‑64 (MQ=255) aGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg > 1:3740688/1‑63 (MQ=255) gCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTAcc > 1:4752602/1‑71 (MQ=255) gCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTAcc > 1:4912293/1‑71 (MQ=255) aaCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGgag < 1:5325712/46‑1 (MQ=37) ccTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGt > 1:2496330/1‑71 (MQ=255) ccTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGt < 1:4707231/71‑1 (MQ=255) cTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATg < 1:1859394/40‑1 (MQ=37) tGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGa > 1:2746987/1‑58 (MQ=255) tGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGtt > 1:3566527/1‑70 (MQ=255) cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt < 1:4355146/69‑1 (MQ=255) cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt < 1:2530449/69‑1 (MQ=255) cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTa < 1:522370/70‑1 (MQ=255) cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTa < 1:206374/70‑1 (MQ=255) cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTa < 1:1184676/70‑1 (MQ=255) || GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTA > AP009048/4371203‑4371334 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |