Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 4,261,351 | 0 | A | G | 45.5% | 11.1 / 11.4 | 25 | F51S (TTC→TCC) | nicP_10 | Porin‑like protein NicP |
Reads supporting (aligned to +/- strand): ref base A (4/8); new base G (10/0); total (17/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.69e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGCTTCACGCCGAACAGGGCGATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGC > KT2440/4261217‑4261478 | ccAGCTTCACTCCGAACAGGGCGATACCATCCAGGCCAAACCCGCTCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGa < 6:154525/135‑1 (MQ=255) ccGCCGAACAGGGCGATACCATCCCGGCCAAAACCGATCCGCCCCGGGGTGTCGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATcgcg < 5:134244/133‑1 (MQ=255) aCGCCGAACAGGGCGATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATcgc < 4:42156/133‑1 (MQ=255) aCGCCGAACAGGGCGATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATcgc < 4:190071/133‑1 (MQ=255) gATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGGCTGGCGGGTGTagtaatagtg > 5:117368/1‑134 (MQ=38) cATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCa < 6:140446/135‑1 (MQ=255) aGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTAGCCCTTGGCGGCCCCCGGGTCGCGGCAAACGAGGTTGAACTAGAAGTTGCGCCATGtg > 6:156500/1‑134 (MQ=38) ggCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCa < 7:102816/135‑1 (MQ=255) cGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGAATCGCGGTTGAAGTAGTAGTTTCGCAATGTCAGCCCAccca > 6:79509/1‑134 (MQ=255) cGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATAGCGGTTTAAGGAGTAGTTTCGCAATGTCAGcccccc > 8:129419/1‑133 (MQ=255) gcAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTgg < 1:100156/134‑1 (MQ=255) aGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGCTTAGCGGGTGTAGTTGTAGTTTGGGAATGTCGGGCCCCCCGTGGcg > 3:117303/1‑134 (MQ=255) aGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGCATCGCGGGTGGTGTAGTAGTTTGGGAGTGTCAGCCCCCCCGTGGcc > 1:155429/1‑134 (MQ=255) aGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGCAATAGCGGTTGAAGTAGTAGTTGCGCAATTTAAGCCCACCCTTTGCa > 8:15516/1‑135 (MQ=255) aGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAACTCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCAACCTTGGCa > 8:178816/1‑135 (MQ=255) ccAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGAATGGCGGTTGAAGTATTATTTGCGCAATGTGATCCCCCCCGTGGGAGCTTCCAAGAAcccc > 6:77326/1‑135 (MQ=255) ccAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATAGCGGGTGAAGTAGTAGTTTCGCAATGGAAGCCCACCCTTGGCATCATGCAACAAGccc > 4:110490/1‑135 (MQ=255) tGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGCAATCGCGGTTGAAGTAGTTGTTGCGCCAATTCCGCCCACCCATGGGAACTTCCAAGAAGccccccgcgt > 6:156367/1‑131 (MQ=255) gAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGAATAGCGGTTGTAGTTGTATTTGGCCGATGTAAGCCCCCCCGTTGGGTCGTCCAAAAACccccccgaagg > 8:81567/1‑130 (MQ=255) ccTTGGCGGCCCCCGGGTCGCGGGAATCGCGGGTGAAGTAGTAGTTGGGCAATGTCCGCCCCCCCATGGGATCTTACAAGAAGCCCCCCGCTTGCGCCCACGCCGGCCAACCACGGCCCCACGCCAGCGACACAg > 6:102139/1‑135 (MQ=255) ccTTGGCGGCCCCCGGGTCGCGGAAATCGCGGGGGAAGTAGTAGTTTCGCAATGTCAGACCCCCCGTGGCATCGTCCAAGACGCCCCCCGATTTCGCCCCCGCCGGCCCACCCACGCCCAACCCCCGCGaacaac > 2:120587/1‑131 (MQ=255) ttGGCGGCCCCCGGGTCGCGGGAATCGCGGGTGAAGTAGTTGTTGCGCAATGTGAGCCCCCCCATGGCCTCTTCCACGAAGCCCCCCTCATGCGCCCACGCCGCCCAACCACGGCCCAACCCCACCGAAACAccc > 4:141678/1‑135 (MQ=255) ttGGCGGCCCCCGGGTCGCGGGAATCGCGGGGGGAGGTGTAGTTGGGCAATGTCAGCCCCCCCGTGGCCTCTTCCCAGAAGCCCCCCGCCTGCGCCCCGGCCGGCCCACCCCCGCCCAACGCCCCCGCAACAccc > 4:200011/1‑135 (MQ=255) ggcCCCCGGGTCGCGGGAATCGCGGGTGAAGTAGTAGTTGCGCAATGTCAGCCCCCCCATGGCATCTTACCAGAAGCCCCCCACTTTCGCCCCGGCCCGCCAACCCACGCCCAACGCCCCCGCAACCGCCacgcc > 2:129143/1‑131 (MQ=255) gTCGGGTAAATCGCGTTTTATGTAGTTGTTGCGAAATGTCCGCCCACCCTTGGCATATTACAATAAGCCCCCCTCTTGCGCCCTGGCGGGCAACCAACGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGc < 6:117368/135‑1 (MQ=255) gTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGCCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGc < 2:155429/135‑1 (MQ=255) | CCAGCTTCACGCCGAACAGGGCGATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGC > KT2440/4261217‑4261478 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |