| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 360,147 | 0 | G | C | 26.2% | 101.3 / 14.9 | 42 | intergenic (‑78/+208) | lacZ/BW25113_RS23320 | beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base G (0/31); new base C (11/0); total (11/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.34e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAC > NZ_CP009273/360059‑360196 | aGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCgg < 2:401725/90‑1 (MQ=255) ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCcgg > 1:73207/1‑90 (MQ=255) cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt < 2:122222/90‑1 (MQ=255) tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCccgccgct > 2:110626/1‑88 (MQ=255) tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt < 1:120388/90‑1 (MQ=255) tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa < 1:19010/90‑1 (MQ=255) caGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa < 1:147890/89‑1 (MQ=255) aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa < 2:452899/90‑1 (MQ=255) aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa < 1:110626/90‑1 (MQ=255) aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa < 1:318789/90‑1 (MQ=255) aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa < 1:185477/90‑1 (MQ=255) aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGGTGCTTTGCaa > 1:174566/1‑90 (MQ=25) gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg < 2:234512/90‑1 (MQ=255) gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGTCCTCCCCCCCCCCCCCCCCGCCGGTGGTTTGCCACCTTc > 1:490395/1‑90 (MQ=19) gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCCGCCCGTGCTTTGCAACGGTc > 2:125159/1‑90 (MQ=19) aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCCAAAccc > 1:98300/1‑89 (MQ=18) aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCGGCCGGTGGTTTGCCACCTTCCAAAccc > 1:234512/1‑89 (MQ=16) tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc < 2:13590/90‑1 (MQ=255) tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCCCAAACTTCAAAAcccc > 1:362054/1‑88 (MQ=18) tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc < 2:405587/90‑1 (MQ=255) ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGGTCAAAAACGCCCCccgcg > 1:3003/1‑90 (MQ=14) gATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCTGCCCCCCCCCCCCCCCGCCCCTGCTTTCCAACCGTCAAAACCGCCCCccgcggg > 2:323563/1‑89 (MQ=9) gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACCTTCAAAACCCCTCCccgcggg > 1:41761/1‑89 (MQ=14) tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt < 2:305156/90‑1 (MQ=255) tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt < 2:41761/90‑1 (MQ=255) tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt < 1:323563/90‑1 (MQ=255) ccccaAATCCCTGGCCGTTTATGTCGGGCCTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt < 2:282078/89‑1 (MQ=11) ttaaTTCCTGGCAGTTTATGGCGGGGGGCCTGCCCCCCCCCCCCCCGGCCGTTGGTTCCCAAAGATTAAAACCCCTCCCCGCGGGttttt > 1:457221/3‑88 (MQ=2) ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt < 1:225879/90‑1 (MQ=38) ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt < 2:48656/90‑1 (MQ=38) cagaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt < 1:448302/87‑1 (MQ=17) aaTTCCTGGCAGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 1:322187/90‑1 (MQ=25) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:361711/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:348833/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:49200/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 1:309042/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:175616/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:158723/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:141668/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:139106/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 1:376932/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 1:180843/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 1:60064/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGAGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:162075/90‑1 (MQ=25) tCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCCCAACGTTCAAAACCCCTCCCCGCGGGTTTTTCCttc > 1:66676/1‑88 (MQ=11) | AGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAC > NZ_CP009273/360059‑360196 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |