| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,518,992 | 0 | C | . | 42.1% | 152.9 / 118.9 | 76 | pseudogene (272/504 nt) | insB1 | IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related |
| Reads supporting (aligned to +/- strand): ref base C (22/22); new base . (15/17); total (37/39) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.20e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACA > NC_000913/4518860‑4519124 | ggggCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCtt > 1:868513/1‑139 (MQ=255) aCGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGt < 2:505595/139‑1 (MQ=255) cGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc > 2:941241/1‑139 (MQ=255) cGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTACGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc > 2:851029/1‑139 (MQ=255) gTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg < 1:941241/139‑1 (MQ=255) gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg < 1:245904/139‑1 (MQ=255) gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg < 1:573511/139‑1 (MQ=255) gcgTATGACAGGATACTGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAAt < 2:344796/139‑1 (MQ=255) gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTAAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAAt < 2:868513/139‑1 (MQ=255) gACAGGCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCt < 1:481529/139‑1 (MQ=11) aGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCCGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGa > 1:82900/1‑139 (MQ=255) gctccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg > 1:1136930/5‑139 (MQ=11) ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaaa > 2:842576/1‑139 (MQ=255) ggaaggcGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc > 1:183670/7‑111 (MQ=11) ggaagaCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc < 2:183670/111‑1 (MQ=11) aCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg > 1:209669/1‑107 (MQ=14) aCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg < 2:209669/107‑1 (MQ=14) cGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGta < 1:139410/85‑1 (MQ=11) cGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGta > 2:139410/1‑85 (MQ=11) tgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg > 2:543585/1‑101 (MQ=255) tgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg < 1:543585/101‑1 (MQ=255) gcgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAg > 2:1330406/1‑139 (MQ=14) cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc < 1:255237/139‑1 (MQ=14) gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc < 2:159301/139‑1 (MQ=255) cACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa > 1:476452/1‑139 (MQ=14) cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg > 1:201431/1‑68 (MQ=255) cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAg > 1:1001820/1‑139 (MQ=255) cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAg > 2:196416/1‑139 (MQ=255) cGTCTTCGGTGAACGCAATCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg < 2:201431/68‑1 (MQ=255) ttCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg < 2:1029621/116‑1 (MQ=17) ttCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg > 1:1029621/1‑116 (MQ=17) tCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCgatat < 2:476452/139‑5 (MQ=14) cGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATata > 1:807673/1‑137 (MQ=11) aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTgg > 2:221387/1‑59 (MQ=17) aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTgg < 1:221387/59‑1 (MQ=14) cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa < 1:328569/113‑1 (MQ=255) cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa > 2:328569/1‑113 (MQ=255) cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAgc > 1:314320/1‑139 (MQ=255) cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg > 1:32630/1‑80 (MQ=17) cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg < 2:32630/80‑1 (MQ=17) atggcgacgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa < 1:1026898/130‑1 (MQ=11) tggcgacgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg < 2:1136930/131‑1 (MQ=11) tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg < 1:1035550/139‑1 (MQ=11) ggCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAg < 1:842576/139‑1 (MQ=255) ggCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAg < 2:314320/139‑1 (MQ=255) ggCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAg > 2:1067379/1‑139 (MQ=255) cgacgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa > 1:507317/6‑135 (MQ=14) cgacgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa < 2:507317/130‑1 (MQ=11) gcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg < 1:620406/65‑1 (MQ=37) gcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg > 2:620406/2‑66 (MQ=25) cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggta < 2:543865/136‑4 (MQ=14) cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggta > 1:543865/1‑133 (MQ=14) ggAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc < 2:1113598/105‑1 (MQ=18) ggAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc > 1:1113598/1‑105 (MQ=21) aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTg > 1:20288/1‑38 (MQ=255) aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTg < 2:20288/38‑1 (MQ=255) aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc > 1:1018757/1‑60 (MQ=255) aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc < 2:1018757/60‑1 (MQ=255) aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCaa < 2:545866/106‑1 (MQ=255) aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCaa > 1:545866/1‑106 (MQ=255) tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa > 1:1040223/1‑129 (MQ=255) tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa < 2:1040223/129‑1 (MQ=255) cttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc > 2:1105713/1‑99 (MQ=255) cttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc < 1:1105713/99‑1 (MQ=255) cttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatc > 1:1260557/1‑139 (MQ=255) ttctGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGc > 1:50843/1‑101 (MQ=17) ttctGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGc < 2:50843/101‑1 (MQ=17) tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc < 2:663405/53‑1 (MQ=255) tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc > 1:663405/1‑53 (MQ=255) gAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTgaga < 1:851029/139‑1 (MQ=255) aGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa > 2:1300984/1‑75 (MQ=255) aGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa < 1:1300984/75‑1 (MQ=255) aGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGac < 2:82900/139‑1 (MQ=255) tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca > 2:779433/1‑139 (MQ=255) gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGc > 1:1232409/1‑61 (MQ=255) gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGc < 2:1232409/61‑1 (MQ=255) | GGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACA > NC_000913/4518860‑4519124 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |