Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,041,386 | 0 | G | C | 27.7% | 96.0 / 24.1 | 47 | G257A (GGA→GCA) | appA | phosphoanhydride phosphorylase |
Reads supporting (aligned to +/- strand): ref base G (25/9); new base C (0/13); total (25/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.54e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATCAATGCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATTTTATTTGCTA > NC_000913/1041323‑1041453 | cATCAATGCTGACGGAGATATATCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAgg < 1:236501/68‑1 (MQ=255) aTCAATGCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGa < 2:778597/68‑1 (MQ=255) tCAATGCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGTGCCGGGGTGGGGAAGGa < 1:970057/67‑1 (MQ=255) tCAATGCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGa > 2:970056/1‑67 (MQ=255) gCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCg > 1:91546/1‑68 (MQ=255) gCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCg > 1:529175/1‑68 (MQ=255) gCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCg > 1:505149/1‑68 (MQ=255) gCTGACGGAGATATTTCTCCTGAAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCg > 1:498752/1‑68 (MQ=255) tGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGAt < 2:505147/68‑1 (MQ=255) tGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGAt < 2:789417/68‑1 (MQ=255) cGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTca > 1:49336/1‑68 (MQ=255) cGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTca > 1:1258838/1‑68 (MQ=255) cGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTca > 2:1245584/1‑68 (MQ=255) ttatTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGTGCAAGGATCACCGATTCACACCa < 1:572291/67‑1 (MQ=255) atatTCCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCa < 1:173415/68‑1 (MQ=255) ttttCTCCTGCAACAAGCACAGGTAATGCCGGAGCCGGGGTGGTGAAGGATCTCCGATTCACACCAGt < 2:473884/67‑1 (MQ=255) atTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGt < 1:957599/68‑1 (MQ=255) atTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGt < 1:1155874/68‑1 (MQ=255) atTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGt < 1:784038/68‑1 (MQ=255) acttCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGGGGGCAAGGATCACCGATTCACACCAgc < 2:1074287/66‑2 (MQ=255) ttCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGAt > 1:667445/1‑56 (MQ=255) ttCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGAt < 2:667443/56‑1 (MQ=255) ttCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTgg > 2:843123/1‑68 (MQ=255) ttCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGTgg < 1:828383/68‑1 (MQ=255) ccTGCGACAAGCGCAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAacg > 2:265770/1‑67 (MQ=255) gCAACAAGCACAGGGAATGCCGGAGCCGGGGGGGGCAATGATCACCGATTCACACCAGTGGAACACCt < 2:91546/68‑1 (MQ=255) aacaaGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGt < 2:1220668/57‑1 (MQ=255) aacaaGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGt > 1:1220669/1‑57 (MQ=255) aacaaGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTg < 1:417141/68‑1 (MQ=255) gCACTGGGAATGCCGGAGCCGGGGCGGGCAAGGATCACCGATTCACACCAGCGGAACACCTTGCTAAg < 2:717492/68‑1 (MQ=255) gCACAGTGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAg < 1:405540/68‑1 (MQ=255) gCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAg > 2:1263505/1‑68 (MQ=255) gCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAg > 1:403571/1‑68 (MQ=255) cacaGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTcaca > 1:886706/1‑46 (MQ=255) cacaGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTcaca < 2:886704/46‑1 (MQ=255) cacaGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGt > 1:803589/1‑68 (MQ=255) acaGGGAATGCCGGAGCCGGGGTGGGGAAGGATCAACGATACACACCAGTGGAACACCTTGCTAAGtt > 1:79352/1‑68 (MQ=255) ggAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCa > 1:335407/1‑68 (MQ=255) ggAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTAACACAAGTGGAACACCTTGCTAAGTATGCa > 2:231537/1‑68 (MQ=255) gAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCAt > 2:39737/1‑68 (MQ=255) gAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCAt > 1:1074545/1‑68 (MQ=255) gAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCAt > 1:776003/1‑68 (MQ=255) aTGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACGCCTTGCTAAGTTTGCATaa > 1:962538/1‑68 (MQ=255) ccGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAAcgc > 1:1149672/1‑68 (MQ=255) ccGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACATTGCTAAGTTTGCATAAcgc > 1:1247229/1‑68 (MQ=255) ggTGGGCAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATTTTAt < 2:1258837/68‑1 (MQ=255) ggTGGGCAAGGAGCACCGATTCACACCAGTGGAACACCTTGCTAAGTCTGCATAACGCGCAATTTTAt < 1:964817/68‑1 (MQ=255) gAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATTTTATTTGCTa > 2:695478/1‑68 (MQ=255) | CATCAATGCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATTTTATTTGCTA > NC_000913/1041323‑1041453 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |