Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,631,772 | 0 | A | C | 27.3% | 142.7 / 27.3 | 66 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
Reads supporting (aligned to +/- strand): ref base A (40/8); new base C (0/18); total (40/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.28e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAA > REL606/2631697‑2631841 | tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca < 1:1417544/76‑1 (MQ=255) aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccac < 1:1040751/76‑1 (MQ=255) tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag > 2:836986/1‑76 (MQ=255) tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag < 2:1944362/76‑1 (MQ=255) tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag < 1:1537594/76‑1 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc > 2:1045965/1‑76 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc < 1:257267/76‑1 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc < 1:1466561/76‑1 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc > 1:1128393/1‑76 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag < 2:46719/76‑1 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag > 2:675820/1‑76 (MQ=255) aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTaa > 2:322389/1‑76 (MQ=255) aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTaa > 1:2237022/1‑76 (MQ=255) aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTaa > 1:2054098/1‑76 (MQ=255) aaaTACCTATTTACCTGCTGGAAGACAACCCCAACGCCCTTTTTATTTCTCCCCCTAATCCCCCCCCAGCCAGTaa < 2:1625009/76‑1 (MQ=255) tGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCCCCTAATCCTCCACCAGCCAGTAACtt < 2:532275/76‑1 (MQ=255) gCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtc > 1:1695750/1‑76 (MQ=255) cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct > 2:441655/1‑76 (MQ=255) cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct > 1:784661/1‑76 (MQ=255) cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct > 2:1986730/1‑76 (MQ=255) cATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTctct > 1:1921813/1‑76 (MQ=255) aTTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtt > 2:1230386/1‑76 (MQ=255) aTTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtt > 1:2096419/1‑76 (MQ=255) tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTTTTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCttt < 2:2132573/76‑1 (MQ=255) tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCCCCACCAGCCAGTAACTTCTCttt < 2:1878111/76‑1 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTCCACCACCAAATCCCCCCCCAGCCAGTAACTTCTCtttt < 2:2048529/76‑1 (MQ=38) taaaCTGCCAGAAGAGTACCCCACCGGCCTTTTTATTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCtttt < 2:776654/74‑1 (MQ=255) aaCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtc > 2:697746/1‑76 (MQ=255) aCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtct > 1:110044/1‑76 (MQ=255) aCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtct > 1:394430/1‑76 (MQ=255) cgCCTGAAGAGTACCCCAACGGCCTTTTTATTTCCCCACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTcg < 1:811305/75‑1 (MQ=255) gCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc > 1:1898896/1‑76 (MQ=255) gCCTGAAGAGTAACCCAACGGCCTTTTTATGTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc > 2:1192994/1‑76 (MQ=255) ccTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgcc > 1:1442049/1‑76 (MQ=255) tGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccgc > 2:2178898/1‑76 (MQ=255) tGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccgc > 2:241875/1‑76 (MQ=255) agagTACCCCAACGGCCTTTTTATTTCCCCCCCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 1:2138503/76‑1 (MQ=255) agagTACCCCAACGGCCTTTTTATTCCACCCCCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 1:1597660/76‑1 (MQ=255) agagTAACCCAACGGCTTTTTTATTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 1:714937/76‑1 (MQ=255) agagTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt > 1:1411036/1‑76 (MQ=255) gagTACCCCACCGGCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg < 1:383527/76‑1 (MQ=255) gagTACCCCAACGGCCTTTTTATTTCACCCCCTAGCCCCCCTCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg < 2:943991/76‑1 (MQ=255) agTACCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGc < 2:1827859/76‑1 (MQ=255) ccaaCCCAGCGGCCTTTTTTTTCCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCg < 1:292334/74‑1 (MQ=255) tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCgc > 1:1412391/1‑75 (MQ=255) tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt > 1:1260669/1‑76 (MQ=255) tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt > 1:698294/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:936904/1‑76 (MQ=255) ccccACCGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCa < 1:1355163/75‑1 (MQ=255) cccACCGGCCTTTTTATTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg < 2:1433123/76‑1 (MQ=255) cccAACGGCCTTTTTATTTCCCCACCTAATCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg < 2:721075/76‑1 (MQ=255) cccAACGCCTTTTTTATTCCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg < 2:1476060/76‑1 (MQ=255) cAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCg > 2:498487/1‑76 (MQ=255) aaCGGCCTTTTTATTCCACCCCCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt < 1:1102085/76‑1 (MQ=255) ggCTTTTTTATTTCCCCACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt < 2:1747652/76‑1 (MQ=255) ggCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt > 1:302680/1‑76 (MQ=255) ggCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt > 2:1248198/1‑76 (MQ=255) gCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTCCTCGCCGCCCTGCGTCAGCGTGttt < 2:104735/76‑1 (MQ=255) ttttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc > 1:423347/1‑76 (MQ=255) ttttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc > 1:2063002/1‑76 (MQ=255) ttttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc > 1:122849/1‑76 (MQ=255) ttttATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc < 1:1096810/76‑1 (MQ=255) tttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCa > 1:1026387/1‑76 (MQ=255) ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc > 1:1928840/1‑76 (MQ=255) ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTCCGTCAGCGCGTTTAGCCAc > 2:793728/1‑76 (MQ=255) tttCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGt < 1:361063/76‑1 (MQ=255) ttCACCACCTAATCCTCCACCAGCCAGTAACCTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa > 1:1588445/1‑76 (MQ=255) cTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat > 2:535426/1‑76 (MQ=255) cTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat < 1:1481256/76‑1 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:222910/1‑76 (MQ=255) tcctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaa > 2:742951/1‑76 (MQ=255) | TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAA > REL606/2631697‑2631841 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |