Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,423,294 | 0 | T | C | 44.4% | 39.6 / 12.6 | 36 | intergenic (+100/+129) | yhdZ/rrfF | putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon |
Reads supporting (aligned to +/- strand): ref base T (3/17); new base C (16/0); total (19/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.35e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC > NC_000913/3423147‑3423443 | ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCGAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg < 1:137717/151‑1 (MQ=255) aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTGCCCGCCCCCCtaccgtc > 2:422607/1‑145 (MQ=255) aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGTCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCTTCCTGCCCGCCACCCTCCGGGCCg < 2:488277/151‑1 (MQ=255) cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCggt > 1:490717/1‑149 (MQ=37) gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt < 1:411370/151‑1 (MQ=255) aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCACCCCCCCGGCCGGTGGTTT‑GCAACAtt > 1:565765/1‑148 (MQ=255) ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTT‑GCAACGTTCaaa > 1:708242/1‑150 (MQ=255) cGCAGGTAATCCATTAATTGAATTTTAGTTCGAAAAGAAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCt < 2:176541/150‑1 (MQ=255) tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTTTCCCACACCCCCCAGGCCCGTGATTC‑GCTACGCTATAATCTGCTGCCcgg > 1:544845/1‑150 (MQ=255) cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGTTTCTTT‑GCCAACTTTAAAACCCCC‑CCCCGCGGGttt > 2:375797/1‑151 (MQ=255) tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCTCCCGGCCGTTTCTTC‑GCAAAGTTTAAAACCCCC‑CCccgcggggtt > 1:126513/1‑144 (MQ=255) aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGCTGGTTCTCCAACATTCAAATCTCCC‑CCcgggggtgt > 2:493710/1‑144 (MQ=255) tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGGCGTTGCTTC‑CCCACGTTCAAATCCCCC‑CCCCGCGGGttttt > 1:60524/1‑144 (MQ=18) ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCggg < 1:568577/150‑1 (MQ=255) gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTC‑GCAAAGTTCAAAACCCCC‑CCCCGCGGGGTTTTTCTACTCggg > 1:426826/1‑147 (MQ=17) tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCACCCCCCCGGGCGGTGCTTC‑CCAACGGTCAAAATCGCT‑CCccgccgggttttccctc > 1:285645/1‑127 (MQ=255) cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCCGCCGTTGCTTC‑GCCACGTTCCAATCCGCC‑CCCCGCGGGTTTTTTCTTCTCcgg > 2:624078/1‑144 (MQ=255) aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAcc < 1:159967/151‑1 (MQ=255) gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCg < 1:667363/151‑1 (MQ=255) cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCGCCCCCCCCCCCCCCGGCCCGTGCTTC‑CCAACGTTCAAAACCCCC‑CCCCGGGGGTTTTTTCTACTcccgggggggttcac > 1:161424/1‑134 (MQ=9) aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGac < 2:181257/151‑1 (MQ=255) aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa < 1:10450/151‑1 (MQ=255) aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa < 1:234260/151‑1 (MQ=255) aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac < 1:432542/151‑1 (MQ=255) tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTCGGACTGCCCCCCCCCCCCCGGCCCCGTTCTTT‑CCCACGCTTCAATCCGCT‑CCCTGCGGGTTTTTTCTTCTCCGGGGAGCGGTGACCGACAAaaaacaaa > 1:247696/1‑146 (MQ=255) ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCTTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaa < 1:533593/150‑1 (MQ=255) cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGCCGGTTCTTT‑GCAACATTCAAATCCCCC‑TCCCGCGGGTTTTTTCCACTCCGCAAAGTGGTCACCGACAAAAAACAGATAaaaaaaaaaac > 1:333243/1‑147 (MQ=9) gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGGCGGTGCTTT‑GCCACGTTCCAATCCGCT‑CCCCGCGGGGTTTTTCTCCTCGGGGGGGGGTTTCCCGCCACACAACAGATAAAAAAAAAAGCGCAggc > 2:291938/1‑149 (MQ=9) gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCTCCTGCCCCCCCCCCCTCCGGCCCTTGCTTC‑GCAAACGTCAAATCCCCT‑CCCGGCGGATTTTTTCCACTCGGGAGAGGGTTCACCGACAAACAACAAATAAAAAAAAAAGCCCAGTc > 1:238718/1‑151 (MQ=255) gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc < 2:67518/151‑1 (MQ=39) tttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc < 2:589558/149‑1 (MQ=38) ttgatGTCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCCCCCCCCGGGGCGGTTCTTC‑GCCAACTTCAAATACCCC‑CCCCGCGGGGTTGTTCTACTCCGGGGGGGGGTTACCGACCAAAAACAGATAAAAAAAAAAGGCCAGCCTTcctcc > 2:534130/1‑148 (MQ=255) tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt < 2:515605/150‑1 (MQ=38) aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa < 2:489646/150‑1 (MQ=38) tGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt > 2:255381/1‑151 (MQ=38) cGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc < 2:66243/151‑1 (MQ=40) cTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc > 2:672817/1‑151 (MQ=40) | CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC > NC_000913/3423147‑3423443 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |