| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,903,692 | 0 | C | A | 41.7% | 51.1 / 20.0 | 36 | intergenic (+318/+1055) | ygcE/queE | putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis |
| Reads supporting (aligned to +/- strand): ref base C (2/19); new base A (15/0); total (17/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.44e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01 | |||||||||||
ACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAA > NC_000913/2903548‑2903829 | aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTc < 1:479476/150‑1 (MQ=255)aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCa < 2:22727/151‑1 (MQ=255)aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑ATATGTAAAACCCCTc > 1:177712/1‑150 (MQ=255) ttAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAAACACCCATTaaactc > 2:231153/1‑144 (MQ=255) tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTAAAAACAACCATTATcc > 2:454965/1‑146 (MQ=255) aaCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCaaa < 1:2865/150‑1 (MQ=255) acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACa < 1:312817/149‑1 (MQ=255) cTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTTAAAACACTCCTTATCCTCAAAACaaaaa > 1:637007/1‑148 (MQ=255) gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTTAAAACACTCATTaacaaaacgggaaa > 2:381613/1‑132 (MQ=255) gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTTAAAAAACACATTAACCTCaaaaaaaaaaaaa > 2:287005/1‑145 (MQ=255) gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGc < 1:381582/134‑1 (MQ=255) gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTGTTAAAACACTCATTAGCCTCAAAAAAAAAACaaaa > 2:573639/1‑148 (MQ=255) gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACt < 2:36907/150‑1 (MQ=255) cTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTTTAAAAACTCATTATCCTCCAaacaaaaacaaa > 1:600722/1‑146 (MQ=255) tAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAg < 1:104075/150‑1 (MQ=255) gAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑tt < 2:73355/151‑1 (MQ=255) aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑ttt < 1:32385/151‑1 (MQ=255) aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑Tttatta < 2:382840/151‑1 (MQ=255) cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATa < 2:70501/151‑1 (MQ=255) cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATa < 1:189636/151‑1 (MQ=255) gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATaa < 1:112232/151‑1 (MQ=255) cAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTTAAAACAATCCTTAACCTCCAAACAAAAAAAAAATTTCACG‑TTTATTATTTTAACCataaaaaaaaaactt > 2:41346/1‑135 (MQ=255) aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTTAAAACACTCATTTTCCTCAAAACAAAAAAAAACATTCAGA‑ATTATTA‑ATTTAACAAATAAAAACaaaaatta > 1:802323/1‑144 (MQ=255) cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTTAAAAAACACATTTAACTCAAAAAAAAAAAAAAAATCCAGA‑TTTATTA‑TTTTACCAAAGAATAAAAAAATTaaaaaaaa > 1:196904/1‑148 (MQ=255) tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAATACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATAACAAGCTTAAAAAATATACt < 2:427961/150‑1 (MQ=255) tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATAACAAGCTTAAAAAATATAc < 1:594672/149‑1 (MQ=255) tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTGTAAAAACACTCCTTTTCCTCAAAAAAAAAAAAAAAATCCCGA‑TTTTTTT‑TTTAAACAAAGAATAAAAAAATTAAAAAAAAAAATGTtttt > 2:308304/1‑151 (MQ=255) ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTAAAAACACTTCTTAAACCCCAAAAAAAAAAAAAATTTCACA‑TTTATTA‑CTTAAACCAAAAAAAAAAAAAATAAAAAAAATACTattattttcaaa > 2:336795/1‑143 (MQ=255) ggATATCTTTAGTCATAAGAATTATTCGTCAATAGCGTATAGTTGCTG‑TTAGGTAATACCACTCCTTAGCCTCAAAACAAAGACAAACTGCCAGA‑TTTATTT‑CTTAACCAATAAATAACAATCTTAAAAAATATTCTGTTATTCTaagtca < 1:54836/151‑5 (MQ=255) atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTATTCAAAACAACTATTAATCCCCAAAAAAAAAAAAAAC‑TACAGAGTTTTTTA‑TTTAACAAAAAAAAAAAAAAATTAAAAAAAAAAATTTTTTTTTaaaaaaaa > 2:255279/1‑150 (MQ=255) atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTAAAAACACCCATTAACCTCCAAAAAAAAAAAAAAATCCAGA‑TTTATTT‑CTTAAACAATGAATAAAAAACTTAAAAAAAAAAATGTTTTTTTTaaaaaaaa > 2:90766/1‑150 (MQ=255) aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑CTATGTAAAAACCCACAATAGACTCAAAAAAAAAACAAACTTCCCAA‑TTTATTT‑ATTAAACAAAGAAAAAAAAACTTAAAAAAAAAAATGTTATTTTTAAAAAAAAAATaa > 2:35024/1‑149 (MQ=255) aaCAGCGTATAGAGGCGG‑TTATTTAAAAACACCCCTTAACCTCAAAAAAAAAAAAAAAATCCAGA‑TTTTTTT‑TTTTAACAAAAAAAAAAAAACCTAAAAAAAAAAATTTTTTTTTTAAAAAAAAAATAAAAAGCATTAAATTAATATTaa > 1:655888/1‑150 (MQ=255) aGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa < 2:148580/151‑1 (MQ=255) gCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaaa < 1:163628/151‑1 (MQ=255) cGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaa < 2:233373/151‑1 (MQ=255) | ACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGA‑TTTATTA‑CTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAA > NC_000913/2903548‑2903829 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |