| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 69,900 | 0 | A | C | 52.6% | 53.8 / 73.8 | 57 | V370G (GTG→GGG) | thiP | thiamine/thiamine pyrophosphate ABC transporter permease ThiP |
| Reads supporting (aligned to +/- strand): ref base A (17/10); new base C (0/30); total (17/40) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.04e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-07 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATA > NZ_CP009273/69811‑69981 | ggCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCa > 2:862655/1‑90 (MQ=255)ggCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCa > 1:32927/1‑90 (MQ=255)ggCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCa > 2:227986/1‑90 (MQ=255)ggCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCa > 1:580287/1‑90 (MQ=255) cGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCt < 2:279326/90‑1 (MQ=255) aTAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGa > 1:116609/1‑90 (MQ=255) tttgtggAGCAGAAAAAAGAAGCCGGAAGCCAGGCCAATCCCCGGCATGGCGAGGATCAACATGCCGCTCACCTCCAGCACCTGACCCGc < 2:440707/89‑1 (MQ=255) ttgttgAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGcc > 2:562862/1‑90 (MQ=255) gttgAGCAGAAAAAAGAAGCCGGTAGCCAGCAAAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAg < 2:494398/90‑1 (MQ=255) aaaaaGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGc > 1:665269/1‑90 (MQ=255) aagaagCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCcgc > 1:682503/1‑90 (MQ=255) agCCGGTAGCCAGCACAACCCCCGGCATGGGGAGGATCAACATGCCGCTCATCCCCACCCCCTGACCCGCCAGCATTTTCTGCCGCGccc < 2:989887/90‑1 (MQ=255) gCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCg > 2:345029/1‑90 (MQ=255) ccGGTGCCCAGCACAACCCCCGGCAGGGCGAGGAACAACACGCCCCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGc < 1:562862/90‑1 (MQ=255) ccGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGc > 2:930592/1‑90 (MQ=255) ccGGTAGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGc < 2:40680/90‑1 (MQ=255) ccGGTACCCACCACAACCCCCGGCAGGGGGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGc < 1:760706/90‑1 (MQ=255) gCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGAc > 1:884182/1‑90 (MQ=255) cacaACCCCCGGCATGGCGGGGACCAACATGCCCCCCACCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 1:279848/90‑1 (MQ=255) cacaACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 2:4154/90‑1 (MQ=255) cacaACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 2:888837/90‑1 (MQ=255) cacaACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 2:109669/90‑1 (MQ=255) cacaACCCCCGGCATGGCGAGGATCAACATCCCCCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 1:749530/90‑1 (MQ=255) cacaACCCCCGGCAGGGCGAGGATCAACATCCCGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 2:841601/90‑1 (MQ=255) aTCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGTCTGCTc > 1:468123/1‑90 (MQ=255) cccccGGCAGGGCGAGGATCAACAGGCCCCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTcc < 2:975138/89‑1 (MQ=255) cGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATag > 2:232354/1‑90 (MQ=255) ggCATGGGGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 2:797819/90‑1 (MQ=255) ggCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 2:32927/90‑1 (MQ=255) ggCAGGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 2:213654/90‑1 (MQ=255) ggCAGGGCGAGGACCAACATCCCCCTCACCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 1:104660/90‑1 (MQ=255) ggCAGGGCGAGAATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 2:712089/90‑1 (MQ=255) tagggCGAGAACCACCATGCCCCTCATCCCCGGCCCTTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagcag < 2:715660/87‑1 (MQ=255) cATGGGGGGGATCAACATGCCGCTCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagcag < 2:645169/90‑1 (MQ=255) tGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 2:201730/76‑1 (MQ=255) tGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt > 1:201730/1‑76 (MQ=255) cGAGGACCACCATGCCGCTCATCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 1:59754/90‑1 (MQ=255) cGAGAATCAACATGCCCCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 2:244653/90‑1 (MQ=255) accAACAGGCCCCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc < 1:544660/88‑1 (MQ=255) aTCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc < 2:884182/90‑1 (MQ=255) aTCAACATGCCCCTCACCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc < 2:121869/90‑1 (MQ=255) cAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg < 2:53292/90‑1 (MQ=255) cAACATGCCCCTCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg < 2:116609/90‑1 (MQ=255) aaaCATCCCCCTCATCCCCACCACCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg < 2:31847/89‑1 (MQ=255) aaaCAGGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg < 2:852566/89‑1 (MQ=255) accgtGCCGCCCCTCCCCAGCCCCTCCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc < 1:700078/86‑1 (MQ=255) cttcccGCTCACCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCAc < 2:735698/86‑1 (MQ=255) tcccGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTa < 2:905217/88‑1 (MQ=255) ccGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACg > 1:745811/1‑90 (MQ=255) ccGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACg < 2:250807/90‑1 (MQ=255) cccTCACCCCCAGCACCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:359534/88‑1 (MQ=255) cGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc > 2:789349/1‑90 (MQ=255) cGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:889946/90‑1 (MQ=255) cGCTCATCCCCAGCCCCGACCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 1:224337/90‑1 (MQ=255) gCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGca > 1:905217/1‑90 (MQ=255) aTCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGc > 1:210561/1‑42 (MQ=255) aTCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGc < 2:210561/42‑1 (MQ=255) tctcCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATa > 2:254057/1‑90 (MQ=255) gCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca > 1:859302/1‑63 (MQ=255) gCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 2:859302/63‑1 (MQ=255) gCACCTGACCCGCCAGCATTTTCTGCCGAGCCCGc > 2:464214/1‑35 (MQ=38) aCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcgcg < 2:730564/41‑1 (MQ=255) aCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcgcg > 1:730564/1‑41 (MQ=255) | GGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATA > NZ_CP009273/69811‑69981 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |