| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,465,739 | T→C | pseudogene (348/2513 nt) | ydbA → | pseudogene, autotransporter homolog; interrupted by IS2 and IS30 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,465,739 | 0 | T | C | 98.0% | 169.5 / ‑4.4 | 50 | pseudogene (348/2513 nt) | ydbA | pseudogene, autotransporter homolog; interrupted by IS2 and IS30 |
| Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (14/35); total (14/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 | |||||||||||
TGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATAC > NC_000913/1465597‑1465885 | tGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAAc > 2:148212/1‑151 (MQ=255) aTGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCt < 2:427684/151‑1 (MQ=255) ttCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCgg > 2:1207/1‑150 (MQ=255) ccACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGgaga < 1:260652/151‑1 (MQ=255) cacaGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGgagat < 1:526096/151‑1 (MQ=255) gctgcGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGataata < 2:402886/150‑1 (MQ=255) gctgcGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGataat < 1:205052/149‑1 (MQ=255) gcGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCgg < 2:529626/151‑1 (MQ=255) aaTCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGAt > 2:634950/1‑151 (MQ=255) aGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACg < 1:106403/151‑1 (MQ=255) gTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTacac < 2:723666/151‑1 (MQ=255) aaCGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAAc < 1:747038/151‑1 (MQ=255) acgacgATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAAcc < 1:738769/151‑1 (MQ=255) cgacgaTGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCg < 1:342447/151‑1 (MQ=255) cgacgaTGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCg < 1:8880/151‑1 (MQ=255) acgaTGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGt > 1:301368/1‑151 (MQ=255) tgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCg > 1:704695/1‑151 (MQ=255) gatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCg < 1:261132/150‑1 (MQ=255) gatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGt < 1:473802/151‑1 (MQ=255) gatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGttt < 2:184650/150‑1 (MQ=255) catcatcatcGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaat < 1:732629/150‑1 (MQ=255) cGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGtt < 1:283076/150‑1 (MQ=255) cGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCAAAACCGGTTTCGTTTAATAATGACGTTa < 2:5455/151‑1 (MQ=255) gCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTAt > 1:655899/1‑151 (MQ=255) gCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTAt > 1:608113/1‑151 (MQ=255) gCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTAt < 2:386396/151‑1 (MQ=255) caacaaTTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTAtt < 2:643910/151‑1 (MQ=255) aacaaTTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATtc < 2:63368/151‑1 (MQ=255) aaTTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCg < 2:457616/151‑1 (MQ=255) aTTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCg < 2:598950/150‑1 (MQ=255) tctcCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGAt > 1:766888/1‑150 (MQ=255) ccACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATaaaa < 2:129051/151‑1 (MQ=255) tCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAAc > 1:230296/1‑148 (MQ=255) cccACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAg < 1:109983/149‑1 (MQ=255) cccACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGaa < 2:450226/151‑1 (MQ=255) cccACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGaa < 1:416657/151‑1 (MQ=255) ccACCTACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGaaa < 1:7963/151‑1 (MQ=255) ccACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGaaa < 1:597693/151‑1 (MQ=255) ccACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGaaa < 1:420980/151‑1 (MQ=255) cACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGaaa > 2:176921/1‑150 (MQ=255) aCCGACACCCCCTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGaaaa > 1:135629/1‑150 (MQ=255) acacCCCCTGATGATGAATCAGACCACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCCAACCGGTTTCGTTTCATAATGACGTTATTATCGATAAACCAGAAAAAACGt > 2:313226/1‑151 (MQ=255) cccccTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTaa < 1:371498/151‑1 (MQ=255) ccccTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTaa < 1:719466/150‑1 (MQ=255) cccTGATGATGAATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACt > 2:734378/1‑151 (MQ=255) gatgatgaATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATTc < 2:702083/151‑1 (MQ=255) gatgaATCAGACGACACTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATTcgcg < 1:240674/151‑1 (MQ=255) cacTCCAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACt < 1:216408/151‑1 (MQ=255) cAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATAc > 1:162463/1‑151 (MQ=255) cAGTCCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATAc < 2:693994/151‑1 (MQ=255) | TGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATAC > NC_000913/1465597‑1465885 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |