Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,631,772 | 0 | A | C | 35.5% | 47.0 / 13.8 | 31 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
Reads supporting (aligned to +/- strand): ref base A (17/3); new base C (0/11); total (17/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.30e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG > REL606/2631697‑2631845 | tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTcca > 2:320741/1‑76 (MQ=255) cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTtct > 2:1220652/1‑76 (MQ=255) cATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTctct > 2:53719/1‑76 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCtttt > 1:795012/1‑76 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCC‑TTTTATTTCACC‑ACCTGATCCTCCACCAGCTAGTAACTTCTCttttt < 1:199075/76‑1 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑CCCTAATCCCCCCCCAGCCAGTAACTTCTCttttt < 1:199059/76‑1 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCCCCAGCCAGTAACTTCTCttttt < 2:929104/76‑1 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCttttt > 1:401265/1‑76 (MQ=255) tAACGGCCTGAAGAGTACCCCACCGGCCTTTTTATTCCCCC‑CCCTAACCCCCCCCCAGCCAGTAACTTCTCttttt < 1:199056/76‑1 (MQ=255) aaCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtc > 2:611233/1‑76 (MQ=255) gAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCgccgcc > 2:622005/1‑76 (MQ=255) agagTACCCCAACGCCCTTTTTATTTCCCC‑CCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 2:1137350/76‑1 (MQ=255) agagTAACCCACCGGCTTTTTTATTTCCCC‑CCCTATTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 1:143879/76‑1 (MQ=255) agagTAACCCACCGGCTTTTTTATTTCCCC‑ACCTAACCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 2:1137344/76‑1 (MQ=255) gagTACCCCACCGGCCTTTTTATTCCCCC‑CCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg < 2:250505/76‑1 (MQ=255) tAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt > 1:424181/1‑76 (MQ=255) ccccAACGCCCTTTTTATTCCCCC‑CCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCa < 2:951063/75‑1 (MQ=255) ggCCTTTTTATTTCACC‑ACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt < 2:597861/76‑1 (MQ=255) gCCTTTTTATTCCCCC‑CCCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 1:1121564/76‑1 (MQ=255) gCCTTTTTATTCCCCC‑CCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:241203/76‑1 (MQ=255) cttttttATTTCCCC‑CCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 2:724622/74‑1 (MQ=255) ccTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa > 2:147167/1‑76 (MQ=255) tttttttttTCCCCAACCAAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 2:724634/75‑1 (MQ=255) ttttttATCTCCCC‑CCCAACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg < 2:812286/75‑1 (MQ=255) cTTTTTATTCCCCC‑CCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg < 2:816126/76‑1 (MQ=255) cTTTTTATTCCCCC‑ACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg < 2:812297/76‑1 (MQ=255) ttATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc > 2:519991/1‑76 (MQ=255) ttATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc > 2:668921/1‑76 (MQ=255) aTTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTg > 1:1109096/1‑76 (MQ=255) tCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTaa > 2:902577/1‑76 (MQ=255) acc‑accTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCCCTGTAAcg > 1:1012117/1‑75 (MQ=255) acc‑accTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACa > 1:1012118/1‑76 (MQ=255) c‑accTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACaaa > 1:696142/1‑76 (MQ=255) cccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAat < 2:1146207/75‑1 (MQ=255) ccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAata > 2:357576/1‑76 (MQ=255) tAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt > 1:858796/1‑76 (MQ=255) tAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt > 1:621140/1‑76 (MQ=255) aTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaa > 1:453527/1‑76 (MQ=255) cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaaa > 1:1178552/1‑76 (MQ=255) ccaccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAg > 2:959254/1‑76 (MQ=255) | TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACC‑ACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG > REL606/2631697‑2631845 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |