Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,458,741 | 0 | G | A | 41.7% | 10.0 / 24.2 | 24 | A160T (GCC→ACC) | paaG | putative ring 1,2‑epoxyphenylacetyl‑CoA isomerase (oxepin‑CoA forming) |
Reads supporting (aligned to +/- strand): ref base G (5/9); new base A (4/6); total (9/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGGATGATCTGGCAGGTTGTTGATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACA > NC_000913/1458621‑1458872 | gggggCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGAGCGCCATGGGGCTgg > 2:392582/1‑133 (MQ=255) gggCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCa < 2:54021/133‑1 (MQ=255) gggCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCa < 2:537171/133‑1 (MQ=255) aCATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCActgct < 1:348223/133‑1 (MQ=255) tgctgcCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCaa < 2:426464/133‑1 (MQ=255) tgctgcCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCaa < 1:608751/133‑1 (MQ=255) ctgcCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACt < 2:577468/133‑1 (MQ=255) ctgcCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACt < 2:342234/133‑1 (MQ=255) ctgcCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCACTGCTGGGGAATCAACt < 1:266623/133‑1 (MQ=255) cccGTTCATCAAAATTCGCCATAGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAg < 1:176079/133‑1 (MQ=255) cccGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGGGCGCGCACAATGGGCCTGGCACTGCCGGGGAATCAACTGAg > 2:266623/1‑133 (MQ=255) cGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTCAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTg > 1:500635/1‑133 (MQ=255) aaaaTTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAg > 2:365473/1‑133 (MQ=255) aaaTTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAgg < 1:431141/133‑1 (MQ=255) aaaTTCGTCACGGCGCTTAGTAAGGTAGGCTTTATACCCGACTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAgg < 2:203035/133‑1 (MQ=255) tttAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGga < 2:156224/133‑1 (MQ=255) tttAGTAAGTTAGGCGTAATACCCCATTGCGGAGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGga < 2:241281/133‑1 (MQ=255) gTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGgatgat > 1:26518/1‑133 (MQ=255) tAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGga < 2:347002/123‑1 (MQ=255) tAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGga > 1:347002/1‑123 (MQ=255) tAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGAACGAATGGGGGATTATCTGGCa > 1:241281/1‑133 (MQ=255) cGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGGATGATCTGGCAGGTTGTtgatgat > 1:316545/1‑133 (MQ=255) cGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGGATGATCTGGCAGGTTGTtgatgat > 2:563379/1‑133 (MQ=255) cGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGGATGATCTGGCAGGTTGTTGATGATGAAACg > 1:564758/1‑133 (MQ=255) cGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCACCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGGATGATCTGGCAGGTTGTTGATGATGAAACg < 1:365473/133‑1 (MQ=255) cgcCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGGATGATCTGGCAGGTTGTTGATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGaca > 2:475315/1‑133 (MQ=255) | GGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGCTGGCACTGCTGGGGAATCAACTGAGTGCTGAACAGGCGCACGAATGGGGGATGATCTGGCAGGTTGTTGATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACA > NC_000913/1458621‑1458872 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |