Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 5,161,325 | 0 | G | T | 46.7% | 24.8 / 15.5 | 45 | intergenic (‑127/‑377) | CJ019_04854/puuA_3 | hypothetical protein/Gamma‑glutamylputrescine synthetase PuuA |
Reads supporting (aligned to +/- strand): ref base G (18/6); new base T (0/21); total (18/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.92e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGTCATGGGGCTTTTTAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTAT > ACB143/5161192‑5161444 | cagtcaTGGGGCTTTTTAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCATTTGTagcgg > 1:162295/1‑134 (MQ=255) tGGGGGTTTTTTGTGATGCTGGGGAGCCTGCTTTTTCATGATCGGAATTCGGAAAGTCGGGCGATTGTCGCTGATTTGCGGATTCTCTGCTGCCTGCACGCTCCGTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGc < 8:371653/135‑1 (MQ=255) cTTTTTAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCaa > 7:92129/1‑135 (MQ=255) cTTTTTAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCaa > 4:124147/1‑135 (MQ=255) cTTTTTAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTAGCGATAGGGCCGCaa > 2:337965/1‑135 (MQ=255) tAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCg > 1:115583/1‑135 (MQ=255) tAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGAAAGGCCg > 2:19589/1‑135 (MQ=255) tAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTAGCGCTAGGGCCGCAAGGccc > 7:330072/1‑134 (MQ=255) tAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCGGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCg > 3:214294/1‑135 (MQ=255) gCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACt > 6:169883/1‑136 (MQ=255) gTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATGGGCGCAGATTTGCGGATTCTCTGCGGCCGGCACGCTCCCTGTGG‑AAGCGGCCTTGTGGCGCGATAGGGCCGCAAGGCCGCTTCCACTTTg < 1:337965/135‑1 (MQ=255) cgCCTGGTTTTTCATGATCGGAATTCGTAAAGGCGGGCGATTGTCGCAGATTTGCGTATTCTCTGCTGCCTGCCCGCGCCGTGTGG‑ACGCGGCCTTGTGGCGCTATAGGGCCGAAAGCCCGCTTCCACTTTGAcc < 4:313086/134‑1 (MQ=255) aGCCTGTTTTTTCATGAGCGGCATTGGTAACGTGGGGCGTTTGTCGCATTTTTGCGTTTTCTCTGCGGCCTGCACGCCCCCTGGGT‑AAGGGGCTTTGTGGCGCTATTGGCCGGAATGCCCGCTTCCACTTTGACCg < 6:128191/136‑1 (MQ=255) gCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑CAGCGGCCTTGTGTAGCGATAGGGCCGCAAGGCCGCTTGCACTTTGACCgtggt > 1:167356/1‑135 (MQ=255) tttCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACGTTGACCGTGGTGCttt > 6:172326/1‑136 (MQ=255) tttCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAAGCCGCTTACACGTTTACCGTGGTGCtt > 5:246399/1‑135 (MQ=255) ggCATGCGGAACGGGGGGCGATTGTCGCATTTTTGCGGATTCTCTGCTGCCTCCACGCCCCCTGTGC‑AAGCGGCCTTGTGGCGCTATAGGGCCGCATGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAAt < 6:397837/136‑1 (MQ=255) aaGTCTGGCGATTGGCGCAGATTTGCGGAGTCTCTGCTGCCGGCACGCTCCCTGTGG‑ACGCGGCTTTGTGTCGCTATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCaa < 6:303693/136‑1 (MQ=255) aaGTCGGGCGTTTGTCGCTTATTTGCGGTTTCTCTGCTGCGGGCACGCTCCCTGTGG‑AAGCGGCGTTGTGTCGGTTTAGGGCCGCATGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCaa < 5:313647/136‑1 (MQ=255) cgTCGGGCGTTTGTCGCATATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGGCGCTATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCaa < 1:383487/134‑1 (MQ=255) atccGGGCTATTTGCGCAGATTTGCGGAGACTCTGCGGCCTGCACGCTCCCTGTGG‑ACGCGGCCTTGTGTCGCGATAGGGCCGAAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCaa < 1:341465/132‑1 (MQ=255) aggcGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCTATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCaa < 1:325645/132‑1 (MQ=255) aGTGGGGCTATTGTCGCAGATTTGCGGTTGCTCTGCGGCCTGCACGCGCCCTGTGG‑ACGCGGCGTGGTGGCGCGATAGGGCCGAAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCaa < 1:26287/135‑1 (MQ=255) aGTCGGGCTATTGTCGCAGATTTGCTTATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGTCCTTGTGTCGCTATAGGGCCGAAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCaa < 7:280034/135‑1 (MQ=255) aGTCGGGCGATTGTCGCTGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCTATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCaa < 1:53801/135‑1 (MQ=255) gggCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATATGGCCGCAAGGCCGCATACACTTTGACAGTGGTGCTTTCAACAAGAATGTTTAGCCAAGTTCg > 5:117526/1‑136 (MQ=255) tCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCATGTGTCGCGATAGGGCCGCAAGGCCGCTTACACGTTGCCCGTGGTGGGTTCAACACGGATTGTTATCCAAAGGGGATGTaaat > 1:203241/1‑134 (MQ=255) ttGCGTATTCTGTGCGGCCTGAACGGCCCCTGTGG‑AAGCGGCCTTGTGTCGCTATAGGGCCGCATGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCAt < 7:322818/135‑1 (MQ=255) ttGCGTATTCTCTGCTGCCTGCACGCTCCGTGTGG‑AAGCGGCCTTGTGTCGCTATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCAt < 5:201393/135‑1 (MQ=255) tGCGGTTGCTCTGCTGCCTGCACGCGCCCTGTGT‑AAGGGCCGTTTTGTCGCTATAGGCCCGAATGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCAt < 1:345973/134‑1 (MQ=255) tGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCTATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCAt < 1:39808/134‑1 (MQ=255) gCGTATTCTCTGCTCCCGGCACGCTCCCTGTGG‑ATGGGGCCTTGTGTCGCTATAGGGCCGAATGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCa < 6:140065/136‑1 (MQ=255) ggATTCTCTGCTGCCTGAACGCTCCCTGTGG‑ACGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCAc < 1:237601/135‑1 (MQ=255) gATTCTCTGCTGCCTGCACGAGCCCTGTGG‑AAGCGGCATTGTGTCGCTATAGGGCCGAAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACt < 7:74462/135‑1 (MQ=255) tgctgcCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTa < 5:370310/134‑1 (MQ=255) gccggCTCCCTGTGT‑ACGGGGCGTTTTGTCGCTATTGGGCGGAATGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCtt < 3:160083/131‑1 (MQ=255) gcACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTa > 5:379665/1‑136 (MQ=255) gcACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACGTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGGTTTAATCAACCt > 1:187016/1‑134 (MQ=255) gcACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGAAAGGGCCCCAAGGCAGCTTACACTTTGAACGTGGTGCTTTCAACACGAATGTTTAGCCAAGTTCCATGTAAAGGTTGCCAATGACTGCTTTAATAAAccgt > 7:186794/1‑133 (MQ=255) gcACGCGCCCTGTGG‑AAGCGGCCTTGTGTCGCTATCGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCtt < 4:128480/135‑1 (MQ=255) gaaCGCGCCCTGTGC‑AAGCGGCCTTGTGCCGCTAATGGGCCGAATGCCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCtt < 7:156660/133‑1 (MQ=255) cccggccccTGTGGTACGGGGCGTTGTGGGGCTATCGGCCCGAAAGCCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCtt < 8:345640/129‑1 (MQ=255) tgtgG‑AAGAGCCCTTGTGTCGCGAAAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGGGCTTTCAACACCAATGTTTAGCCAAGTGCGATGTAAATGGTGCCAATACCAGCGTTAATCAACCGTACGCCGCAtgg > 6:98772/1‑135 (MQ=255) tgtgG‑AACCGGCCTTTTGTCGCGATAGGGCCGCAAGGCCGCTTACACGTTGAACGGGGTGGATTAAACACGAAGGTTTATCCAACTGCCATTTAAAGGTTTCCATTCACATCTTTAATCAACCTAACGCCGCAtg > 2:181612/1‑135 (MQ=255) tgcggcgTTGTGGCGCTATAGGGCGGCACGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTAt < 5:349213/135‑1 (MQ=255) gCGGCCTTGTGTCGCTATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTAt < 7:380878/135‑1 (MQ=255) | CAGTCATGGGGCTTTTTAGTGATGCAGGTGAGCCTGCTTTTTCATGATCGGCATTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGG‑AAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTAT > ACB143/5161192‑5161444 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |