Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,764,783 | 0 | C | A | 44.7% | 53.3 / 16.9 | 38 | intergenic (‑397/+151) | sufA/ydiH | Fe‑S cluster assembly protein/uncharacterized protein |
Reads supporting (aligned to +/- strand): ref base C (13/8); new base A (0/17); total (13/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.20e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGG > NC_000913/1764640‑1764932 | gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGc < 2:461117/151‑1 (MQ=255) gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTa < 2:449456/151‑1 (MQ=255) aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑Cccgcc < 1:648639/151‑1 (MQ=255) aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGt > 1:252252/1‑151 (MQ=255) gCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAac > 1:702127/1‑151 (MQ=255) gCGTACAAAAATAAAAAACCCCTTAAGTCTGTGCGACACGGGTTTAAGGGTTTCTCCCCCATCCGGCGCTTATCCCCGGCCCTCTCAGTGGTTTAGCTCTTGAG‑GGGGCGATAAGAATAATCTCATAAAGCTAC‑CCCGCCGTTTTAAcaca < 1:804224/151‑1 (MQ=255) gTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACaa > 2:328215/1‑150 (MQ=255) aaaaaTAAGAAACCCCTTAAGTCTGTGCGACACAGGTTAAAGGTTTTCTCCCCCACCCGGCGTTTATCTCCGGCCCTCCCAGTGGCTTAGCTTTTAAA‑GGGGCGATAAGAATAATCTAAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTg < 2:213763/148‑1 (MQ=255) aTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGAtt > 1:515778/1‑151 (MQ=255) aaacCCTTTAGGTTTGTGCGACACAGGCTTAGGGGTTTCTCCCCCATCCGGCCCTTATCTCCGGCCCTCTCAGTGGTTTAGCTTTTGAG‑GGGGCAATAAAAAAAATCTAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGta < 1:769120/150‑1 (MQ=255) aaaCCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGtatt > 2:444818/1‑151 (MQ=255) ccttaaaGTTTGTGCGACACAGCCTAAGGGGTTTCTACCCCACCCGGCCTTTATCCCCGCCCCTCTCAGTGGCTGAGCTTTTAAA‑GGGGAAAAAAAAATAATCTAAAAAACCAAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTattatt < 2:186995/146‑1 (MQ=255) acacggaattAGGGGTTTCTCCCCCATCGGGCGTTTATCTCCGGCCCTCGCAGTGGCTTAGCTTTGGAA‑GGGGGAAAAAGAATAATCTAATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAg < 1:229951/142‑1 (MQ=255) caGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGTGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCgg > 2:143797/1‑151 (MQ=255) ggCTTAAGGTTTTCTCCCCCACCCGGCGCTTATCCCCGGCCCTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAAAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc < 1:81962/147‑1 (MQ=255) aGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATa > 2:613835/1‑149 (MQ=255) ttCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGt > 2:441376/1‑150 (MQ=255) cccccACCGGGCTCTTATCCCCGCCCCTCTCAGTGGTTTAGCTTTTGAG‑GGGGCGAAAAAAATAATCTAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa < 1:253420/149‑1 (MQ=255) ccccACCCGTCACTAATCCCCGCCACCCTCAGTGGCTGAGCTTTAGAA‑GGGGCGTAAAAAATAATCAAAAAAAGTTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAACATCAGGCGGTAGATAAGCa < 2:615729/143‑1 (MQ=255) ccccACCCGGCGCTTATCCCCGCCCCCCTCAGTGGTTTAGCTTTGGAG‑GGGGCAATAAGAATAATCTAAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa < 1:198327/149‑1 (MQ=255) ggCGCTTATCCCGGGCACTCTCAGGGGCTTAGTTCTAGAG‑GGGGCGAAAAAAAAAATCTAAAAAACCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc < 1:431937/149‑1 (MQ=255) ttATCTCCGCCCCCCCCGGTGG‑TTACCTTTTGGAGGGGGCAAAAAAAAAAATCAAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCg < 2:774463/148‑1 (MQ=255) cgccGGCTCTCTCAGTGGTTTAGCTCTTGAA‑GGGGCAATAAAAATATTCTCAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTgg < 1:214225/149‑1 (MQ=255) ccccGGCACTCTCAGTGGCTTAGCTTTTGAG‑GGGGCAATAAAAATAATCTAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt < 1:734206/147‑1 (MQ=255) ccccGCCACTCTCAGTGGTTTAGCT‑TTGAAGGGGGCAAAAAAAATAATCTCA‑AAAACTAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGa < 2:520573/143‑1 (MQ=255) ccGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGAc > 2:13847/1‑151 (MQ=255) tctcGGTGGTTTAGCTGGTGAA‑GGGGAAAAAAAAATAATCTAAAAAAGCAAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAg < 2:503847/148‑1 (MQ=255) gTGGTTTAGCTCTTGAA‑GGGTCAAAAAAAATAATCTAATAAACCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc < 2:257994/151‑1 (MQ=255) ggCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc > 2:405694/1‑149 (MQ=255) ttAGCTCTTGAG‑GGGGCAAAAAGAATAATCTAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGt < 1:483134/150‑1 (MQ=255) aGCTCTTGAG‑GGGGAGAAAAGAATACACTCATAACACCAT‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGa < 1:594182/151‑1 (MQ=255) aGCTCTTGAA‑GGGGCAAAAAAAATAATCAAATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg < 1:689160/150‑1 (MQ=255) gAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCg > 2:233220/1‑151 (MQ=255) ggggCGAAAAAAATAATCTAAAAAAGATAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGaa < 2:66666/150‑1 (MQ=255) ggcaaaaaGAATAATCTAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt < 1:258614/145‑1 (MQ=255) ggCAATAAAAATAATCTAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt < 1:142545/151‑1 (MQ=255) aaGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc > 2:803574/1‑150 (MQ=255) aaGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc > 2:185457/1‑150 (MQ=255) tcATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAgg > 1:420307/1‑151 (MQ=255) cATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAgg > 2:233698/1‑150 (MQ=255) | GAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGG > NC_000913/1764640‑1764932 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |