Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 445,105 | 0 | A | T | 24.5% | 95.8 / 17.5 | 53 | intergenic (+78/‑210) | PP_0367/PP_0368 | hypothetical protein/acyl‑CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (21/19); new base T (4/9); total (25/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.13e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCA > NC_002947/444957‑445246 | tCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGAcaacaa > 7:128652/1‑149 (MQ=255) tCGAGCTCGGATCAAAGGTCGTCAAGTCATCAGACGAAATGCTTGGCCGGATAATCGACACATACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGA‑GGGCTTGCCCGCGAATGCGTCCGATCAGAcaacaa > 7:122317/1‑149 (MQ=255) cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAt > 5:138843/1‑149 (MQ=255) ggtCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt‑t‑t < 1:163820/147‑1 (MQ=255) gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGCGGGGGTTGCCC‑CGCATTCGGGCGCTCAGCCAACAAat‑t‑tttttg > 1:172746/1‑144 (MQ=255) aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCggg < 6:372069/149‑1 (MQ=255) gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGc < 5:209500/145‑1 (MQ=255) gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGc > 6:209500/1‑145 (MQ=255) cAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTTCCCCCGCATGAGTCCCAAAAGACAACAAAT‑T‑TTTGTGGTgggggcacatc > 3:168199/1‑142 (MQ=255) aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGGTTTCCCGCGAATGCGTCCCATCAGACCACAAAT‑T‑TTTCTGGTCGGGCCTATTagc > 5:152048/1‑149 (MQ=255) aGGCATCCGACGAAAAGCTTGGCCGGATCATCGAGACCTACGCCCGAATACTTGCGACATCACCTAATTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGGGAATGCGTCCGATCAGACAACAATT‑A‑TGTCTGGTCGGGCCTCTTcgc < 6:179947/149‑1 (MQ=255) gCATCCGACTAAATGCTTGGCCGGATCATCGACACCTACGCCTGACTACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGA‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCgg < 8:122317/149‑1 (MQ=255) gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCgg < 7:351435/149‑1 (MQ=255) gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCgg < 8:128652/149‑1 (MQ=255) gCATCCGAAGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCACCGCTACTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTcgc < 1:367943/149‑1 (MQ=255) cATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCggg < 4:208648/149‑1 (MQ=255) aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTTAATACTTGCGTCATCCCCTATTTCAGAGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGAAAACAATT‑T‑TGTATGGTCGGGcc > 4:67905/1‑138 (MQ=255) aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTTCCCGCGAATGCGTCCGATAAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTAGCGGGc > 5:365112/1‑149 (MQ=255) aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACATTT‑T‑TGTCTGGTCGGGcc < 3:67905/138‑1 (MQ=255) aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTG‑T‑TGTCTGATCGGACGCATTCGCGGGc < 7:109178/149‑1 (MQ=18) ccGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTG‑T‑TGTCTGATc < 6:378816/131‑1 (MQ=37) ccGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTG‑T‑TGTCTGATc > 5:378816/1‑131 (MQ=37) gacgaATTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGc < 5:412740/149‑1 (MQ=255) aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCct > 3:419940/1‑130 (MQ=255) aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCct < 4:419940/130‑1 (MQ=255) ccTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGC‑‑CGGCCTTGTCGGCGC‑GGGCTGGCCGGCTAATGCGCCCGATCACACAACATTT‑T‑TGGCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatggc < 6:365112/149‑1 (MQ=37) gCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGAc < 1:73142/149‑1 (MQ=255) tGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGa > 7:356328/1‑149 (MQ=255) aaCTTCCTTCAACCCCTATTTCACCGC‑‑AGGCGTTGAAGCTGC‑GGGCTGGCCGGCTAATGCTTCCGATAATACAAAATTTCTATGTCTGGGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAtt < 1:140668/148‑1 (MQ=255) aTCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTc < 8:356328/148‑1 (MQ=255) tCCCCTATTGCCGCGC‑‑CGGCGTTGTCGGCGG‑GGGCTTGCCCGCTAATGCTTCCTATCATACAACAATT‑T‑TGTCTGGTGGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt < 6:152048/149‑1 (MQ=255) cccatttttCCGCGC‑‑CTGCGTTTTTGGTGG‑GGGTTTGCCCGGTATTGCTTCCTATCATAAAAAATTT‑T‑TGTTGGGTGGGGCCTCTTCGGGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg < 8:25661/143‑1 (MQ=255) cccTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTGGGGCCTCTTCGCGGGCAAGCCCGCACCTACAAGGGTTTGTGGCATGGCATGACTGATTTGGCATTTTtggtgg > 5:228839/1‑149 (MQ=255) ccgcgc‑‑CGCCCTTGTCGGCGC‑GGGCTTGCCCGCTAATGCGTCCGATCCCACAAAATTT‑T‑TGCCGGGGGGGGCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt < 4:412225/146‑1 (MQ=255) gcgg‑‑CTGCCTTGTTGGAGG‑GGGCTTGCCCGGGAATGCGTCCTATCAGACAACATTT‑T‑TGTCTGGTCGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc < 6:228839/149‑1 (MQ=38) tgcctTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGc < 7:47387/85‑1 (MQ=255) cctTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGc < 1:225593/63‑1 (MQ=255) cctTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGc > 2:225593/1‑63 (MQ=255) ctTTTAGGAGG‑GGGCTTGCCCGCTAATGCGTCCGAACACAAAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt < 7:20686/149‑1 (MQ=39) c‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa > 5:454054/1‑149 (MQ=255) gCTTGCCCGCGAATGCGTCCGATCAGAGAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGAGTGATTTGGCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa > 5:72414/1‑146 (MQ=255) ccGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt > 5:167159/1‑149 (MQ=255) ccGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt > 2:95798/1‑149 (MQ=255) gcgAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt > 1:257686/1‑147 (MQ=255) cgAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg < 1:95798/149‑1 (MQ=255) ccGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt > 3:172205/1‑148 (MQ=255) ccGACCAGACAAAATTG‑T‑TGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt > 6:48655/1‑148 (MQ=38) ccGACCAGACAAAATTG‑T‑TGTCTGATCGGACCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACa > 2:305720/1‑130 (MQ=37) ccGACCAGACAAAATTG‑T‑TGTCTGATCGGACCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACa < 1:305720/130‑1 (MQ=37) aTCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt > 3:172074/1‑149 (MQ=255) cAGACAAAATTG‑T‑TGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt > 5:3983/1‑148 (MQ=255) gACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc > 3:308158/1‑148 (MQ=255) acaacaATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt < 2:33879/99‑1 (MQ=255) acaacaATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt > 1:33879/1‑99 (MQ=255) acaacaATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATCACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa > 1:367141/1‑148 (MQ=255) aacaaTT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAATGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTTGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc < 5:433128/136‑1 (MQ=255) aacaaTT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc > 6:433128/1‑136 (MQ=255) caaTT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa > 7:301667/1‑138 (MQ=255) caaTT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa < 8:301667/138‑1 (MQ=255) | TCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGC‑‑CTGCCTTGTAGGAGC‑GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATT‑T‑TGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCA > NC_002947/444957‑445246 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |