| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,779,939 | 0 | C | G | 26.5% | 67.8 / 17.8 | 34 | A930P (GCA→CCA) | ypjA | adhesin‑like autotransporter |
| Reads supporting (aligned to +/- strand): ref base C (6/19); new base G (9/0); total (15/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.54e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCC > NC_000913/2779872‑2780002 | cGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGc < 1:125744/68‑1 (MQ=255) ggTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCt < 2:965614/68‑1 (MQ=255) ggTATATTGCCCGCCAGAGATAACATTTGTGGCGGAGTCCTGCCCCAGGTTTTGCAATGCCCCACACt > 1:615664/1‑68 (MQ=255) gTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGGCCAAcccat > 1:954784/1‑66 (MQ=255) gTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCtt < 2:705242/68‑1 (MQ=255) atTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGGCCCACCCTTGCCa > 1:273157/1‑68 (MQ=255) atTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGGCCCACCCTTGCCa > 1:795210/1‑68 (MQ=255) atTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCa > 2:618019/1‑68 (MQ=255) atTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTAGCATTGCCCCACCCTTGCCa > 2:864291/1‑68 (MQ=255) atTGCCCGCCAGAGTTAACCTTTGTGGCGGAGGCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCa < 1:131198/68‑1 (MQ=255) ttGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCaa > 2:250494/1‑68 (MQ=255) tGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGGCCCACCCTTGCCAAc > 1:169159/1‑68 (MQ=255) tGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGGCCCACCCTTGCCAAc > 2:603287/1‑68 (MQ=255) tGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAAc > 2:334923/1‑68 (MQ=255) gCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAAcc > 1:34494/1‑68 (MQ=255) ccTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCt < 1:323463/68‑1 (MQ=255) ccTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCt < 1:334923/68‑1 (MQ=255) ccTTTGTGGCGGAGTCCTGACCCAGGTTTCGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCt < 1:864294/68‑1 (MQ=255) tgtgGCGGAGTCCTGACCCAGGTTCTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCgg < 2:851576/68‑1 (MQ=255) gtgGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGt < 2:99498/68‑1 (MQ=255) gtgGCGGAGTCCTGACCCAGGTTTAGCATTGCCCCACCCTTGCCAACCGTGTAGTCGCGAGCTTCGGt < 1:495783/68‑1 (MQ=255) ggCGGAGTCCTGACCCAGGTTTTGCATTGGCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTAc > 1:932069/1‑68 (MQ=255) ggAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTg < 1:618019/68‑1 (MQ=255) ggAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTg < 2:811078/68‑1 (MQ=255) gAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGc < 1:383918/68‑1 (MQ=255) gAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGc < 2:832463/68‑1 (MQ=255) gAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGc < 1:535609/68‑1 (MQ=255) gccTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTaa < 1:615339/67‑1 (MQ=255) cTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATa < 2:338029/68‑1 (MQ=255) aCCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTTGAGTCGCGAGCTTCTGTACCTGCTAATACCa < 1:835072/68‑1 (MQ=255) aCCCAGGTCTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCa < 1:495523/68‑1 (MQ=255) cccAGGTTTTGCATTGGCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCaa > 1:35118/1‑68 (MQ=255) cATTGGCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGc > 2:459556/1‑68 (MQ=255) cATTGGCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGc > 1:598556/1‑68 (MQ=255) aTTGCCCCACCCTTGCCAACCGTGGATTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGcc < 1:703221/68‑1 (MQ=255) | CGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCC > NC_000913/2779872‑2780002 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |