| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,570,326 | 0 | T | C | 52.6% | 18.4 / 27.8 | 19 | G517G (GGA→GGG) | yddB | TonB‑dependent receptor plug domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base T (7/2); new base C (0/10); total (7/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.14e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGC > NZ_CP009273/1570237‑1570394 | tCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAAt > 1:760528/1‑90 (MQ=255) tataAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCa < 1:505753/90‑1 (MQ=255) ggCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCa < 1:557435/90‑1 (MQ=255) aTCCTGATAACGCGTAAGATTTTTATTCCCTGATCCCGTTCCCGCCCAGCTATTGGGGATATCACGTAACCCAATATAAAAAATATcccc < 2:94435/90‑1 (MQ=255) gcgTCAGTTTTTTTTTCCCTAACCCCGATCCCGCCCAGCTATGGCGGATATCACGTAACCCCATACAAAGAATATTCCCGCCATAGTAAc < 1:26462/90‑1 (MQ=255) gcgTCAGTTTTTTTTTCCCTAAACCCGATTCCGTCCAGTTTTTGGGGATATCACGTAACCCAATAAAAAAAATTTTCCCGCCATAGTAAc < 2:295215/90‑1 (MQ=255) gcgTCAGATTTTTATTCCCTAATCCCGTTCCCGTCCAGTTATTGGGGATATCACGTAACCCCATATAAAAAATATTCCCGCCATAGTAAc < 1:243837/90‑1 (MQ=255) cgTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACg > 1:77952/1‑90 (MQ=255) ttttttCCCTGATCCCGTTTCCGTCCAGCTTTTGGGGATATCACGTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAAcc < 1:170188/90‑1 (MQ=255) ttttttCCCTGATCCCGTTCCCGCCCAGCTATTGCGAATATCACGAAACCCCATAAAAAAAATATTCCCGCCATAGTAACGGTTATAAcc < 2:77952/90‑1 (MQ=255) ttattaCCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACct > 1:368506/1‑90 (MQ=255) aCCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGt > 1:547100/1‑90 (MQ=255) ccTGATCCCGTTCCCGTCCAGTTTTTGCGAATATCACGAAACCCAATATAAAAAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTa < 1:769790/90‑1 (MQ=255) taaTACCCATCCCCCCCAGCTTTTGCGAATATCACATACCCCAAAAAAAAAAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAAt < 2:758059/88‑1 (MQ=255) gATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAg > 1:408095/1‑90 (MQ=255) gATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAg > 2:248685/1‑90 (MQ=255) ccGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGttt > 1:795538/1‑90 (MQ=255) gTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGa > 2:320702/1‑90 (MQ=255) ccAGTTTTTGGGGATATCACTTACCCCCATAAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt < 1:382379/90‑1 (MQ=255) gCTATTGCGAATATCACGAAACCCAATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc < 2:54737/90‑1 (MQ=255) | TCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGC > NZ_CP009273/1570237‑1570394 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |