| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,466,210 | 0 | G | A | 27.8% | 45.1 / 12.6 | 36 | pseudogene (819/2518 nt) | ydbA | putative outer membrane protein, N‑terminal fragment |
| Reads supporting (aligned to +/- strand): ref base G (11/15); new base A (8/2); total (19/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.52e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACG > NC_000913/1465962‑1466448 | gaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCTCAGATGGCTCAACGAAGTGTCAATATAATCACAACGGTGAGCGTGTGATTACGGGTGATAATGCCTCAGTAAACAACAATGGCAAATCGACCGTTGACGGTAAAAATTCCACCGGTACGGAAATCAGCGGGAATAACGGGAAATTGATGCAGGACTGCGAGCTGGATGTCAGCGGCGGTTGGCACGGTATTGATATCACCGGTGACATCGCGAc < 1:155299/251‑1 (MQ=255) cTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGCAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGCCAGCGCGACGGGGGATAACCAGGGCACCAt > 2:37374/1‑251 (MQ=255) aGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCAAAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCCGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACCGCGCAACGGGGGATAACAAGGGCCCCATGACCg > 2:200708/1‑251 (MQ=255) gAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGCCCGTCACCGATCCGGAGTCCATGGGt > 1:283799/1‑250 (MQ=255) aaaaataATCACAACGGTGGGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGAt < 2:93096/247‑1 (MQ=255) aaTATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATc < 2:283799/251‑1 (MQ=255) aTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGt < 2:260934/251‑1 (MQ=255) acaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGAACGGAAATCACCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGGGGATAACAAGGGCACCATGACCGTCCACGGTCCGGAGTCCATGGTTATCCAGGTTGACGGTGAc > 2:249903/1‑251 (MQ=17) aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACa > 1:295103/1‑249 (MQ=255) cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGGGACCAgg > 2:95758/1‑251 (MQ=255) aGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATc < 1:142145/251‑1 (MQ=255) cTTGTTATTACGGGTGATAATGCAACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGt > 1:127224/1‑251 (MQ=255) attGTTATTAAGGGTGATAAAGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGt < 2:282486/250‑1 (MQ=11) ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGtta > 1:327953/1‑247 (MQ=255) tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACAGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCCCGGAATTGATATCACCGGCGACAGCGCAACGGTGGCTAACAAGGGAACCATGACCGGCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGCCAAGGCCACCGTCAAt > 2:49857/1‑250 (MQ=14) ttaCGGGAGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGCTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGTGACAAGGCCATTGTCAATaa > 1:322147/1‑250 (MQ=255) cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACCAGGCCGTCGTCAATAACGa > 1:131800/1‑250 (MQ=255) acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTCATAACGGGAACGTCCTTCAGGACGGCGCTCTGGATGTCAGCCGCGGAGGACCCGCTATTGATATCACCGGTGACAGCGCGACGGTAGCTAAAAAGGGCACCATAACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAAGCGAGAGCACCAt > 1:342146/1‑251 (MQ=255) ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg < 1:153774/251‑1 (MQ=255) aaaGATTCCACCGGTACGGAAAGCAACGGGAATAACGGGAAAGTGATGCAGGACGGCGATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGAACCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGTTGGCACCGGCACCCAGATTAATGGtg < 2:104439/250‑1 (MQ=18) ggAAATCCACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac > 1:117469/1‑250 (MQ=255) gTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGCTGTCAGCGGCGGCGGTCACGGTATTGATATAACCGGTGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCCTGGGTATCCAGATCGACGGCGACAAGGCCATCGTCAATAACGAAGCCGAGAGCACCATCACCAACGGTGGCCCCGGCAACCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa > 1:314428/1‑250 (MQ=17) aataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac < 1:162303/251‑1 (MQ=255) aataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTAAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAAGCCACGGCAAACAACAACGGCAAAac > 1:340903/1‑251 (MQ=255) cGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCCCGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGTCACCGGCACCCAGATGAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg < 2:127224/251‑1 (MQ=255) gAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGGTATCACCGGGGGCAGGGCGACGGTGGAGAACTAGGGGGCCATGACCGCCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGTCAAGGCCATCGTCAATAACGAAGGCGAGCGCACCATCACCAACGGTGGCACAGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGttt > 2:11550/1‑250 (MQ=18) tGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGAGACCGCACCGCAGTCCATCAGTATCCAGCTCGACGGTCACAAGCCCACCGTCAATAACGAAGACGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc > 1:49717/1‑251 (MQ=14) aTTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa < 1:37374/251‑1 (MQ=255) gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGTTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCGCCATGACCGTCGCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAGGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTc < 2:364701/251‑1 (MQ=255) agggcAGCGGCGGCGGCCACGGTACTGATAGCACCAGTGACAGGGAGACGGGGGAGAACAGGGGCACCATTACCGACACCGATCAGGAGTCCATGGGTAACAAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACTGAGATTAATGGTGATGAGGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGTAGTCCACCGTTACGGaa < 2:151903/247‑1 (MQ=17) gcggctgTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACggg > 1:328516/1‑250 (MQ=255) cggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCCCCATGACCGTCACCGATCCGGAGTCCCTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTACCCCGGTACGGAAATCAACGGtaatc > 1:223842/1‑250 (MQ=255) cggTCACGGACTTGATATCCCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGACCCCGAGCCGGAGGCCATGGGTATCCAGATTGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata < 2:71605/251‑1 (MQ=18) atCACCGGTGACAGCGCGACGGTGAATAACAAGGGCACCATGACCGACACCGATCCGGAGGCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc < 1:118133/250‑1 (MQ=255) cccGGTGACAGCGCGAAGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAAGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGa < 1:22206/250‑1 (MQ=255) cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCATATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGTAAATCAACGGTAATAACGGAAAAGTTATTCAGGACg < 1:108411/251‑1 (MQ=255) | GAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACG > NC_000913/1465962‑1466448 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |