Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 126,279 | 0 | G | A | 26.8% | 96.8 / 34.6 | 56 | A581T (GCG→ACG) | CJ019_00134 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (15/26); new base A (3/12); total (18/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCTGGCTGGGGACAAGGTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGCATTCCGTTCTGGCGAACGCAGCGACCAGGTGCTCAA > KT2440/126145‑126412 | aTCTGGCTGGGGACAAGGTCCTCAAGATCGTGGCCGACTAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATAGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTg > 1:31559/1‑135 (MQ=255) tCTGGCTGGGGACAAGGTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTgc < 1:199664/135‑1 (MQ=255) gctggGGACAAGGTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCaa > 7:282055/1‑135 (MQ=255) ggggACAAGGTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTg > 2:16678/1‑135 (MQ=255) ggggACAAGGTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGGGATTTCATCGCGCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTg < 5:250551/135‑1 (MQ=255) ggACAAGGTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTg < 8:295092/133‑1 (MQ=255) aaGGTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAg < 8:4968/135‑1 (MQ=255) ggTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGgt < 6:69514/135‑1 (MQ=255) ccTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGtt < 3:220030/135‑1 (MQ=255) ccTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGGTTTGGCGGCGCGGAATTCGTCATGCTGCTGCCGCCGACTACGCCCCCAGTCGCTGCGCAAGGGGCCGAGGTGtt > 3:34224/1‑135 (MQ=255) aGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGCTTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCCGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTg > 1:106605/1‑135 (MQ=255) aGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTg < 1:187326/135‑1 (MQ=255) aGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTg < 8:51606/135‑1 (MQ=255) aGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTg < 6:203031/135‑1 (MQ=255) aGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTg < 8:163896/135‑1 (MQ=255) aGATCGTGGCCGACCAGTTGCGCAAGCGGCTGCGGGGCCGTGAGTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTg < 7:174322/135‑1 (MQ=255) cGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCcac < 5:297083/135‑1 (MQ=255) cGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCcac < 3:279871/135‑1 (MQ=255) ggCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGGCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTAAGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTGGCGTGCCACAGt > 7:319156/1‑135 (MQ=255) ccGACCAGTTGCGCAAGCGTCTGTGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCg > 6:26913/1‑135 (MQ=255) ccGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAggc > 6:44212/1‑132 (MQ=255) cGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGa < 3:21311/135‑1 (MQ=255) gactggCGTCGGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCa < 5:191489/130‑1 (MQ=255) cgcAAGCGTCTGCGGGGCCGTGATTTAATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCa < 2:2877/135‑1 (MQ=255) gcAAGCGTCTGTGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCa < 5:26913/134‑1 (MQ=255) gcAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCAAAGCCTGCCCa < 1:16678/134‑1 (MQ=255) ggggCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGCCTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAgg > 5:159784/1‑135 (MQ=255) ggggCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGGAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCa < 1:229820/132‑1 (MQ=255) gTGATTGCAGCGCCCGTTTTGGCGGCGAGGAATTCGTCCTGCTGCTGCCGCAGACTACGCCCGCAGGCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCCCTTCACgggcgcac < 4:129563/135‑6 (MQ=255) tttCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACGTCAAGGGAGCACGggt > 6:319206/1‑135 (MQ=255) cGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGAt < 3:289694/135‑1 (MQ=255) gCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATc > 8:275437/1‑135 (MQ=255) gCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCACTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATc > 1:310172/1‑135 (MQ=255) cccGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATc < 4:193261/134‑1 (MQ=255) ccGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATcac < 7:238661/135‑1 (MQ=255) ccGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATcac < 8:308012/135‑1 (MQ=255) tGGCGTCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCt < 8:101885/134‑1 (MQ=255) tGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCt < 1:326699/92‑1 (MQ=255) tGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCt > 2:326699/1‑92 (MQ=255) tcgAGGAATTGTGCATGCTGCTGCCGATGAATACGCCGGCTGGCGCGGCGAATGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCTTTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCg < 7:126977/134‑1 (MQ=255) gcgAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCg < 1:56153/135‑1 (MQ=255) aGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGcc < 6:136904/135‑1 (MQ=255) ggAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATca < 1:124510/119‑1 (MQ=255) ggAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATca > 2:124510/1‑119 (MQ=255) gAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTg < 3:23176/135‑1 (MQ=255) aaTTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTgg > 7:2136/1‑135 (MQ=255) gTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGc < 4:82046/134‑1 (MQ=255) gTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGc < 4:280603/134‑1 (MQ=255) gTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGc < 8:2136/134‑1 (MQ=255) gctgcCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTg > 4:64197/1‑53 (MQ=255) gctgcCGCAGACTACGCCCGCAGTCGCTACGCAAGTGGCCGAGGTGTTGCGTg < 3:64197/53‑1 (MQ=255) cccGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGGCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGCATTCCGTTCTGGCGAACGCAGCGAc < 1:242195/135‑1 (MQ=255) cccGCAGGCGCTGCGCAAGTGGCCGAGGTGTTGCGGGCCACAGTCGAACCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCAGCGGCCTGGGTGCATTCCGTTCTGGCGAACGCAGCGAc < 4:34224/135‑1 (MQ=255) ccGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGggt > 7:258641/1‑79 (MQ=255) ccGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGggt < 8:258641/79‑1 (MQ=255) ccGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGCATTCCGTTCTGGCGAACGCAGCGa > 2:297360/1‑133 (MQ=255) tCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGCATTCCGTTCTGGCGAACGCAGCGACCAGGTGc > 3:325518/1‑135 (MQ=255) tGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGCATTCCGTTCTGGCGAACGCAGCGACCAGGTGCTCaa < 3:238873/135‑1 (MQ=255) gcgcAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGCATTCCGTTCTGGCGAACGCAGCGACCAGGTGCTCa > 3:111350/1‑133 (MQ=255) | ATCTGGCTGGGGACAAGGTCCTCAAGATCGTGGCCGACCAGTTGCGCAAGCGTCTGCGGGGCCGTGATTTCATCGCCCGTTTTGGCGGCGAGGAATTCGTCATGCTGCTGCCGCAGACTACGCCCGCAGTCGCTGCGCAAGTGGCCGAGGTGTTGCGTGCCACAGTCGAAGCCTGCCCATTCCACTTCAAGGGAGAACGGGTGGTGATCACCACCTCCATCGGCCTGGGTGCATTCCGTTCTGGCGAACGCAGCGACCAGGTGCTCAA > KT2440/126145‑126412 |
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: ‑ |