Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,562,123 | 0 | A | T | 46.4% | 29.5 / 12.3 | 28 | I189N (ATT→AAT) | glgC | glucose‑1‑phosphate adenylyltransferase |
Reads supporting (aligned to +/- strand): ref base A (10/5); new base T (0/13); total (10/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTAGATACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAA > NZ_CP009273/3562040‑3562203 | cGTAGATACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑GTTTTC‑AACGAATTCGATAATTTTAt > 1:162208/1‑90 (MQ=255) aGATACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGt < 2:590499/90‑1 (MQ=255) aTACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTtc < 1:204875/90‑1 (MQ=255) aTACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTtc < 1:623849/90‑1 (MQ=255) cATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTTCtcatc > 1:407209/1‑90 (MQ=255) ttGCTGGGTTCGTTCGGCATTGGGGGGGGGTTAGCAGT‑TTTTTC‑AACGAATTCGATATTTTTATCTTTCTCATCAACCGCCATAACGCCa < 1:103575/90‑1 (MQ=255) ttGCTCGGTTCGTTCGGCTTGGAGGGGGGGTTAGCAGT‑TTTTTC‑AACAAATTCGATATTTTTATCTTTCTCATCAACCGCCATAACGCCa < 1:69832/90‑1 (MQ=255) ttGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCa > 2:552035/1‑90 (MQ=255) gCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCaaa > 2:519188/1‑90 (MQ=255) gATCGTTCGGCATTGAGGGGGGGTTAGCGGT‑TTTTAA‑AACAAATTCGATATTTTTATCCTTCTCATCAACCGCCATAACGCCAAATGCgg < 2:479149/90‑1 (MQ=255) gATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGt < 1:213203/60‑1 (MQ=255) gATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGt > 2:213203/1‑60 (MQ=255) gATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCgg > 2:164175/1‑90 (MQ=255) ggcgggggTTTAGCAGT‑TTTTTA‑AACGAATTCGATATTTTTATCCTCCTCATCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATc < 2:67871/88‑1 (MQ=255) ggcgggggGTTAGCAGT‑TTTTAC‑AACGAATTCGATATTTTTATCTTTCTCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATc < 1:376135/88‑1 (MQ=255) gACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGt < 2:169006/46‑1 (MQ=255) gACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGt > 1:169006/1‑46 (MQ=255) gcgggggGTTAAGAGGTTTTTTC‑AACAAATTCAATATTTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATc < 2:539386/89‑1 (MQ=255) ggggggTTAGCAGT‑TTTTTC‑AACGAATTCGATATTTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGt < 2:259269/87‑1 (MQ=255) gcggGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTa > 2:353135/1‑90 (MQ=255) ggggTTAGCAGT‑TTTTTA‑AACAAATTCGATATTTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTAc < 2:197416/89‑1 (MQ=255) ggggTAAGCGGT‑TTTTTA‑AACAAATTCGATATTTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTAc < 1:98206/89‑1 (MQ=255) ggTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTg > 2:465099/1‑90 (MQ=255) gttttttA‑AACAAATTCGATATTTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACaa < 1:100852/88‑1 (MQ=255) tttttA‑AACGAATTGGATATTTTTATCGTTCTCATAAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGc < 2:90178/90‑1 (MQ=255) ttttCAAACAAATTCAATATTTTTATCTTTCCCACCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGc < 1:149100/90‑1 (MQ=255) ttttC‑AACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCa > 1:171327/1‑90 (MQ=255) aaTTCAATATTTTTATCGTTCTCATCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCaa < 2:407209/90‑1 (MQ=255) | CGTAGATACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTC‑AACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAA > NZ_CP009273/3562040‑3562203 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |