| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 851,606 | 0 | G | A | 37.1% | 146.0 / 90.1 | 116 | S3F (TCT→TTT) | ybiU | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (61/12); new base A (1/42); total (62/54) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.70e-19 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
Only 100 of 180 total aligned reads displayed.
TGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTATAAATAA > REL606/851531‑851681 | tGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGGCATGTGGCGCTGGTTAAAGC‑Ag > 1:1682147/1‑76 (MQ=255)tGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑Ag < 1:3594660/76‑1 (MQ=255) gcgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGGAATGTGTCGCTGGTAAAaga‑aga > 1:7375657/1‑76 (MQ=255) gcgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGGAATGTGTCGCTGGGAAAAGT‑AGa > 1:6655955/1‑76 (MQ=255) gcgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGa > 2:7511825/1‑76 (MQ=255) gcgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGa < 2:486155/76‑1 (MQ=255) gcgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGa < 1:7341373/76‑1 (MQ=255) gcgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGa > 2:7388422/1‑76 (MQ=255) gcgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAaga‑aga > 2:6392163/1‑76 (MQ=255) cgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCCATGTGTCGCTGGTAAAAAT‑AGaa > 2:5147335/1‑76 (MQ=255) cgcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAAT‑AGaa > 2:6715640/1‑76 (MQ=255) gcCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAg < 2:1233916/76‑1 (MQ=255) ccGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGcc < 1:5978376/76‑1 (MQ=255) ccGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGGGATCGGCAGGGCATGTGTCGCCGGGAAAACG‑AGAAGcc > 2:4920513/1‑76 (MQ=255) gCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCAt < 1:2890493/76‑1 (MQ=255) gCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCAt > 1:3843352/1‑76 (MQ=255) cAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATa < 1:4505724/76‑1 (MQ=255) gcgTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGT‑AGAAGCCATAaccc > 1:1793844/1‑75 (MQ=255) gcgTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGG‑AGAAGCCATAacca > 2:3604324/1‑76 (MQ=255) cgTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGT‑AGAAGACATAaacac > 2:5447724/1‑76 (MQ=255) cgTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGT‑AGAAACCATAaccac > 1:37883/1‑76 (MQ=255) gTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCCTAaccacc > 1:2995000/1‑76 (MQ=255) tcTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCCCctc > 2:4364757/2‑76 (MQ=255) gCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGGAATGTGTCGCTGGTAAAAGT‑AGAAGCCCTAAACACctc > 2:5622058/1‑76 (MQ=255) gCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGT‑AGAAGCCATAAACACctc > 1:4120729/1‑76 (MQ=255) cTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGT‑AGAAGCCATAACCACccct > 2:1568539/1‑76 (MQ=255) tCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCGATGTGTCCCGGGGAAAAGT‑AGAAGCCATAACCACCTCTCa > 2:3819743/1‑76 (MQ=255) tCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCa > 2:2316864/1‑76 (MQ=255) cATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAAT‑AGAAGCCATAACCACCTCTCAt > 2:2308258/1‑76 (MQ=255) cATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTAGTAAAAGT‑AAAAGCCATAACCACCTCTCAt < 1:7714508/76‑1 (MQ=255) aTCTGACGGATAGCTTCTTGGTGAGCGCCGGGAAATGTGCCGCGGGAAAAAAA‑AAAAGCCATAACCACCTCTCAtt < 2:1296113/76‑1 (MQ=255) aTCTGACGGATAGCGGCTTTGTGATCGGCAGGCAATGTGTCGCTGTAAAAAGT‑AAAAGCCATAACCACCTCTCAtt < 1:5690110/76‑1 (MQ=255) tCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGT‑AGAAGGCATAACCACCTCTCATTa > 1:6014239/1‑76 (MQ=255) cggACGGATAGCTGTTTTGTGATCGGCAGGAAATGTGTCGCTGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTaa < 1:641999/74‑1 (MQ=255) cTGGCGAATAGCTGTTTTGTGAGCGGCAGGAATTGTGCCGGTGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTaa < 2:1683761/76‑1 (MQ=255) cTGACGAATAGCTGTTTTGGGATCGGCAGGCAATGTGTCGCTGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTaa < 1:3211961/76‑1 (MQ=255) tGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAAc > 1:2144250/1‑76 (MQ=255) tGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAAc > 1:5225195/1‑76 (MQ=255) gggATGGCTGTTTGGTGAGCGGCGGTCATTGTGTGGCTGAAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACttt < 2:1029266/75‑1 (MQ=255) gggATAGCTTCTTTGTGATCGGCAGGCAATGTGTCGCTGGAAAAAGT‑AAAAGCCATAACCACCTCTCATTAACttt < 1:2455370/75‑1 (MQ=255) ggATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTc > 1:6433200/1‑76 (MQ=255) tAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCAtt > 2:2166766/1‑76 (MQ=255) tAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCAtt > 1:2218469/1‑76 (MQ=255) aGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTc > 2:4010367/1‑76 (MQ=255) gctgttTTGTGATCGGCAGGCAATGTGTCGCTGGAAAAAGT‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCg < 1:5711512/76‑1 (MQ=255) gctgctTTGTGATCGGCAGGCAAGGTGTCGCTGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCg < 1:4734691/76‑1 (MQ=255) gctgctTTGTGATCGGCAGGAAATGTGTCGCTGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCg < 1:2000440/76‑1 (MQ=255) gctTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGaaa > 1:1920309/1‑76 (MQ=255) gctTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGaaa > 1:788253/1‑76 (MQ=255) gctTTGTGATCGACAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGaaa > 2:5460302/1‑76 (MQ=255) ctTTGTGCTCGGCAGGAAATGTGTCTCTGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAAc < 1:3756541/76‑1 (MQ=255) tgtgtgGTCGGCAGAAAATGGGGCGCTGGAAAAAGT‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAAcc < 1:1944656/74‑1 (MQ=255) ttGTGATCGGAGGGCAATGGGTGGTGGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCa < 1:1316883/76‑1 (MQ=255) tgtgATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCAt > 2:7357069/1‑76 (MQ=255) gtgATGGGAAGGCAATTTGCCGCTTGAAAAAGT‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCAtt < 1:7808503/76‑1 (MQ=255) gtgATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCAtt > 1:4033489/1‑76 (MQ=255) tgATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAAACATTc > 2:333235/1‑76 (MQ=255) gagcGGCAGGCAATGTGCCGCTGGAAAAAGA‑AAAAGCCATGACCACCTCTCATTAACTTTCATTCGAAACCATTCt < 2:571785/73‑1 (MQ=255) gATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCt > 2:1113177/1‑76 (MQ=255) gtCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTg < 1:2699247/75‑1 (MQ=255) gtCGGCAGGCAATGGGTTGCTGGAAAAAGTAAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCt < 2:5922006/75‑1 (MQ=255) ggcaggcaATGTGTCGCTGGTAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGAtt < 1:5867729/76‑1 (MQ=255) gcgggcaTTGTGTCGTGGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGAttt < 1:1721821/73‑1 (MQ=255) caggaaTTGTGTCGTGGAAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTg < 1:4464540/76‑1 (MQ=255) aggcaATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGc > 1:2659827/1‑76 (MQ=255) attgtgTCGGTGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgt < 1:5678204/74‑1 (MQ=255) aaTGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgt > 2:3762371/1‑76 (MQ=255) aaTGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgt > 1:3398405/1‑76 (MQ=255) tggtgtCGCTGAAAAAAAA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg < 2:863747/74‑1 (MQ=255) aTGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg > 2:1664846/1‑76 (MQ=255) aTGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg > 2:2939563/1‑76 (MQ=255) aTGTGTCGCTGAAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg < 1:1177695/76‑1 (MQ=255) tgtgTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGa > 1:3418236/1‑76 (MQ=255) tgtgTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGa > 1:6190635/1‑76 (MQ=255) tgtgGCTGTAAAAAGT‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACa < 1:960264/76‑1 (MQ=255) tgccGCTGAAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACa < 2:7688295/73‑1 (MQ=255) gtCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACaa > 1:3022304/1‑76 (MQ=255) tCGTTGGAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAat < 1:5389014/76‑1 (MQ=255) tCGTTGGAAAAAAA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAat < 1:6840842/76‑1 (MQ=255) cGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAata > 1:853691/1‑76 (MQ=255) cGCTGAAAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAata < 1:7114811/76‑1 (MQ=255) cTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAAtata > 1:3385860/1‑76 (MQ=255) cTGGTAAAAGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAAtata < 1:7511825/76‑1 (MQ=255) ggTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAg > 1:2778164/1‑76 (MQ=255) gaaaaaaGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAg < 2:4760907/73‑1 (MQ=255) gaaaaaaGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAg < 1:7126314/73‑1 (MQ=255) gTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGt > 1:6745007/1‑76 (MQ=255) tAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGtt > 2:7303537/1‑76 (MQ=255) aaaaaGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGtt < 1:5544336/75‑1 (MQ=255) aaaaaGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGtt < 2:6910439/75‑1 (MQ=255) aaaGA‑AAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTtat < 2:5953489/76‑1 (MQ=255) aaGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTtata > 2:3823227/1‑76 (MQ=255) aaGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTtata > 1:1363516/1‑76 (MQ=255) aaGG‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTtata < 1:5859514/76‑1 (MQ=255) gT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTataaa > 2:5353583/1‑76 (MQ=255) t‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTataaat > 1:4187246/1‑76 (MQ=255) t‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTataaat > 1:2072417/1‑76 (MQ=255) t‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTataaat > 2:1878963/1‑76 (MQ=255) aaGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTataaat < 1:4560526/75‑1 (MQ=255) gAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTataaataa > 2:5070789/1‑76 (MQ=255) | TGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGT‑AGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTATAAATAA > REL606/851531‑851681 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |