| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 888,014 | 0 | T | G | 34.7% | 164.5 / 14.3 | 75 | noncoding (44/81 nt) | rybB | small regulatory RNA RybB |
| Reads supporting (aligned to +/- strand): ref base T (15/34); new base G (26/0); total (41/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.49e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAA > NC_000913/887966‑888062 | aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGtt < 1:838208/50‑1 (MQ=255)aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGtt < 1:6379114/50‑1 (MQ=255)aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGtt < 1:342873/50‑1 (MQ=255)aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGt < 1:346061/49‑1 (MQ=255)aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGt < 1:8176263/49‑1 (MQ=255) ttttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGa > 1:4985332/1‑50 (MQ=255) ttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATg > 1:5188706/1‑50 (MQ=255) tAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATg < 1:6625287/49‑1 (MQ=255) tAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATgg > 1:2316721/1‑50 (MQ=255) aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATgg < 1:5513069/49‑1 (MQ=255) aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATggg > 1:2646912/1‑50 (MQ=255) aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATggg > 1:2087463/1‑50 (MQ=255) aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATggg > 1:7799559/1‑50 (MQ=255) gATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGc < 1:3904459/50‑1 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCt < 1:5196389/50‑1 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCt < 1:1965710/50‑1 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATggggt > 1:3527424/1‑48 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGcc > 1:6960556/1‑49 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGca > 1:2961879/1‑49 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCt > 1:223946/1‑50 (MQ=255) tAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTc > 1:2741149/1‑50 (MQ=255) aCAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTcc < 1:2021322/49‑1 (MQ=255) aCAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCca < 1:6629231/50‑1 (MQ=255) aaaaaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcaca < 1:1774881/50‑1 (MQ=255) aaaaaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcaca < 1:966860/50‑1 (MQ=255) aaaaaaCCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCCCcac > 1:1671780/1‑49 (MQ=255) aaaaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACaa < 1:6035396/50‑1 (MQ=255) aaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACaaaa > 1:6958980/1‑50 (MQ=255) aaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACaaaa < 1:2958054/50‑1 (MQ=255) aaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACaaaa < 1:2215842/50‑1 (MQ=255) aaaCCCATCAACCTTGAACCGAAATGGCGGGGGTGGAGGAATCCACaaac > 1:7585554/1‑49 (MQ=255) aaaCCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACaaa > 1:4024895/1‑49 (MQ=255) aaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCCCAAAAt > 1:25041/1‑50 (MQ=255) aaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAt < 1:4506367/50‑1 (MQ=255) aaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCCCCACAAAAt > 1:1409634/1‑50 (MQ=255) aaCCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAt > 1:4212915/1‑50 (MQ=255) cccATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATgg < 1:8073720/50‑1 (MQ=255) cccATCAACCTTGAACCGAAATGGCGGGGGTGATTGGCCCCCCAAAAAgg > 1:6551168/1‑50 (MQ=255) ccATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg < 1:1068371/50‑1 (MQ=255) cATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCCCAAAATtgg > 1:2015809/1‑49 (MQ=255) cATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATggg > 1:6263672/1‑49 (MQ=255) cATCAACCTTGAACCGAAATGGCGGGGGGGAAGGGGTCCCCAAAAAgggg > 1:4510663/1‑50 (MQ=255) aTCAACCTTGAACCGAAATGGCGGGGGTTATGGGCCCCACAAAAAGGGGa > 1:692572/1‑50 (MQ=37) aTCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATggggg > 1:2775123/1‑49 (MQ=255) aTCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGa > 1:2942378/1‑50 (MQ=255) tCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGAc < 1:6229448/50‑1 (MQ=255) tCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGa > 1:5572373/1‑49 (MQ=255) cAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCa > 1:2216825/1‑50 (MQ=255) aaCCTTGAACCGAAATGGCGGGGTTTGTGGGCCCCCCAAAATGGGGACAt > 1:6726078/1‑50 (MQ=255) aaCCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACAt < 1:1084853/50‑1 (MQ=255) aaCCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACAt > 1:936235/1‑50 (MQ=255) aaCCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACAt > 1:5023953/1‑50 (MQ=255) aaCCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGACAt > 1:1077190/1‑50 (MQ=255) aCCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATc < 1:1572167/50‑1 (MQ=255) aCCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATc < 1:5449826/50‑1 (MQ=255) ccTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCa < 1:6150280/50‑1 (MQ=255) ccTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAATGGGGGCATc > 1:3900142/1‑49 (MQ=38) cTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATaaa < 1:3800307/50‑1 (MQ=255) cTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCaa < 1:5650393/50‑1 (MQ=255) cTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCa > 1:2642843/1‑49 (MQ=255) gAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCaaa < 1:688555/48‑1 (MQ=255) gAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAAGGGGACAACAAAg > 1:7130894/1‑49 (MQ=255) aaCCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGaa < 1:5814830/50‑1 (MQ=255) ccGAAATGGCGGGGTTGATGGGCTCCACAAAAAGGGGACATCAaaaaaa > 1:7674770/1‑49 (MQ=255) cGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAg < 1:6767892/50‑1 (MQ=255) cGAAATGGCGGGGTTGATGGGCTCCACAAAAAGGGGGCATCAAAGaaaa > 1:3280712/1‑49 (MQ=255) aaaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCa < 1:1800049/50‑1 (MQ=255) aaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAg < 1:2270755/50‑1 (MQ=255) aTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGt < 1:1189215/50‑1 (MQ=255) tGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTg > 1:7167908/1‑50 (MQ=255) tGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTg > 1:5928300/1‑50 (MQ=255) tGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTg < 1:24633/50‑1 (MQ=255) gCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGc < 1:2532794/50‑1 (MQ=255) cGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa > 1:293492/1‑50 (MQ=255) cGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa < 1:3388082/50‑1 (MQ=255) ggggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCAc < 1:7468323/50‑1 (MQ=255) ggggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCAc < 1:6415102/50‑1 (MQ=255) gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTaa > 1:554637/1‑50 (MQ=255) gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTaa > 1:2812891/1‑50 (MQ=255) gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTaa < 1:1384014/50‑1 (MQ=255) | AGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAA > NC_000913/887966‑888062 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |