| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,903,692 | 0 | C | A | 51.1% | 34.8 / 74.3 | 45 | intergenic (+318/+1055) | ygcE/queE | putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis |
| Reads supporting (aligned to +/- strand): ref base C (6/16); new base A (23/0); total (29/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.15e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-04 | |||||||||||
TTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTT > NC_000913/2903543‑2903834 | tttttACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAc < 1:466318/150‑1 (MQ=255) gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAACCACCC‑ATTat > 2:412138/1‑149 (MQ=255) gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTC‑ATTAga > 1:351428/1‑150 (MQ=255) gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACAAAC‑ATTaaa > 1:7511/1‑149 (MQ=255) gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAg < 1:703273/150‑1 (MQ=255) gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTC‑ATTaa > 2:110645/1‑149 (MQ=255) gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAAAACTC‑Ctttg > 1:543027/1‑148 (MQ=255) atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAACC‑Tc > 2:426190/1‑150 (MQ=255) aCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAAc < 1:640333/151‑1 (MQ=255) aCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAAc < 2:559932/151‑1 (MQ=255) cTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACTC‑AATTGAC‑CCcaaacaaaaaa > 1:448863/1‑148 (MQ=255) cTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATC‑ATTAACC‑TCAAAACAAAAAAAAAAttt > 1:294220/1‑149 (MQ=255) aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACTC‑ATTTTCC‑TCAAAAAAAAAAAAAACTTCCAGAt > 1:26439/1‑151 (MQ=255) aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGAt < 1:328475/149‑1 (MQ=255) aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGAt < 1:235562/149‑1 (MQ=255) aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGAt < 1:94902/149‑1 (MQ=255) aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGAt < 1:520403/149‑1 (MQ=255) aaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTAc < 2:503072/151‑1 (MQ=255) aGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACCC‑ATTAACC‑TCCAAACAAAAACAAAATTCCAGAttttttttt > 2:363239/1‑148 (MQ=255) aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACt < 1:705226/149‑1 (MQ=255) cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTC‑ATTTAAC‑TCAAAAAAAAAACAAAAATCCCGATTTTTTACTTaa > 1:109644/1‑151 (MQ=255) aTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGaa < 2:118418/151‑1 (MQ=255) aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑a < 1:249708/151‑1 (MQ=255) cAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAACACTC‑ATTATCC‑TCAAAACAAAAAAAAAATTTCAGATTTATTAATTTAACAAAGAATAAC‑AAAATTaa > 1:556574/1‑151 (MQ=255) cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATAc < 2:80312/151‑1 (MQ=255) tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTg < 1:627765/151‑1 (MQ=255) aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGt < 1:131042/151‑1 (MQ=255) gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACACTC‑ATTTGCC‑TCAAAAAAAAAAAAAAATTCCAGATTTATTTCTTTAACAATGAATAAA‑AAACTTAAAAAAAATACTTTTAt > 1:30845/1‑151 (MQ=255) gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAAAACTCAATTAACC‑‑CAAAAAAAAGAAAAACTTTCAGATTTATTTCTTAAACAAAGAAAAAA‑AAACTTAAAAAATATAATGTTAt > 1:289020/1‑151 (MQ=255) aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAGGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTc < 2:260226/151‑1 (MQ=255) ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTC‑ATTATCC‑TCAAAAAAAAAACAAAATTCCCGATTTTTTACTTTACCAATGAATAAC‑AAACTTAAAAAATATAATGTTTTTTTaaaaaa > 2:21253/1‑151 (MQ=255) tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCCCTC‑ATTAGCC‑TCCAAAAAAAAACAAACTTCCCGATTTATTTCTTAAACAATGAATAAC‑AAACTTAAAAAAAATACTGTTATTCTAAAGaaaaaaa > 2:185759/1‑150 (MQ=255) tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAAAAATC‑ATTAGAC‑TCAAAAAAAAGACAAACTTCCAGATTTATTTATTAAACAAAGAATAAA‑AAACTTAAAAAAAATAATGTTATTTTaaaaaaaaaaa > 1:695048/1‑150 (MQ=255) cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACCC‑CTTTGCC‑TCAAAAAAAAAACAAACTTCCACATTTTTTTCTTAAACAAAGAATAAC‑AAACTTAAAAAAAATACTGTTTTTTTAAAAAAAAAAATaaaaag > 2:510540/1‑146 (MQ=255) cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACAC‑ATTTTCC‑TCCAAACAAAAAAAAAATTCTAGAATTTTTACTTAAACAAAGAAAAACAAAGATTAAAAAA‑ATAAAGTTTTTTTTAAAAAAAAATTaaaaag > 1:684904/1‑146 (MQ=255) tAAGAACTATTCATCAACAGCGTATAGAGGCGGGTATGTAAAAACACAC‑AATTTCC‑TCAAAACAAAGAAAAACTTTCAAATTTATTACTTAAACAAAGAAAAAA‑AAACTTAAAAAAAATAATGTTATTTTTAAAAAAAAAATaaaaaacaa > 1:212722/1‑144 (MQ=255) aaCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAAt < 1:227212/151‑1 (MQ=255) cTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACCC‑CTTACCC‑TTAAAACAAAAAAAAACATCCAGATTTATTAATTAAACAATGAATAAC‑AAAATTAAAAAAAATAATTTTTTTTTAAAGAAAAAAATAAAAAGCAATaaatt > 1:414712/1‑148 (MQ=255) aTTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTC‑ATTAGCC‑TCAAAACAAAGAAAAACTTTCAGATTTATTACTTAAACAAAGAATAAA‑AAACTTAAAAAAAATAATGTTTTTTTTAAGAAAAAAATAAAAGGCATTAATTCaa > 2:435008/1‑151 (MQ=255) aTTCATCAACAGCGTATAGAGGCGGTTATGTAAAAAAACCC‑AATTGGC‑TCAAAAAAAAAAAAAAAATCCCGAGTTTTTTCTTTAACAAAGAAAAAA‑AAACTTAAAAAAAATACAGTTTTTTTTAAAAAAAAAATAAAAAGGATTAATTTaa > 1:331012/1‑151 (MQ=255) catcaACAGCGTATAGAGGCGGTTATGTAAAAACACTC‑CTTAGCC‑TCCAAAAAAAAAAAAAACTCCCGATTTATTTCTTTACCAAAGAAAAAA‑AAACTTAAAAAAAAAAATGTTATTTTTAAGAAAAAAATAAAAAGCATTAATTTAatat > 1:675467/1‑151 (MQ=255) catcaACAGCGTATAGAGGCGGTTATGTAAAAACACTC‑ATTAGCC‑TCCAAAAAAAAACAAAACTCCAGATTTTTTACTTTAACAAAGAATAAA‑AAACTTAAAAAAAATAATGTTATTTTAAAAAAAAAAATAAAAAGGATTAATTTAatat > 2:217953/1‑151 (MQ=255) aCAGCGTATAGAGGCGGTTTTGTAAAAACACCC‑CATAGCCTTCAAAAAAAAAAAAAAATTCC‑GATTTATTTCCTTAACAAAGAAAAAA‑AAAAATAAAAAAAAAAATTTTTTTTTAAAAAAAAAAATAAAAAGCAATTATTCaaaaattaaa > 1:338887/1‑143 (MQ=255) aCAGCGTATAGAGGCGGTTATGTAAAAACAATC‑ATTTACC‑TCCAAAAAAAAAAAAAATTTCAGATTTTTTAATTTAACAAAGAAAAAA‑AAAAATAAAAAATAAAATATTTTTTTTAAAAAAAAAATAAAAAACATTAAATTAAAATTtaaa > 1:339026/1‑151 (MQ=255) aTGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt > 2:546918/1‑151 (MQ=255) aTGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt > 2:646476/1‑151 (MQ=255) | TTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCC‑TCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTT > NC_000913/2903543‑2903834 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |