Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 716,832 | 0 | T | G | 27.5% | 103.7 / 17.6 | 51 | A31A (GCA→GCC) | ECB_00663 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (7/30); new base G (14/0); total (21/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.00e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAAAGAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAT > REL606/716757‑716907 | gAAAAGAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGt < 2:767855/76‑1 (MQ=255) aaaaGAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTg < 1:1143393/76‑1 (MQ=255) aaTGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTa < 1:1293325/76‑1 (MQ=255) aTGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTAc < 2:893755/76‑1 (MQ=255) gAACATTAGAGACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCa < 2:578193/76‑1 (MQ=255) aaCATTAAATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCAt < 2:1380062/76‑1 (MQ=255) cATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTATCAACg > 2:345780/1‑76 (MQ=255) aTTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGt < 2:282028/76‑1 (MQ=255) aTTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGCTTACCATCGt > 2:610811/1‑76 (MQ=255) ttAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTTCCATCgcg > 2:833107/1‑74 (MQ=255) ttAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTACCATCggg > 2:959635/1‑74 (MQ=255) gATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAAc < 2:833142/76‑1 (MQ=255) aaaCGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTAACATCCTGGACa > 1:551626/3‑76 (MQ=255) aTACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTTCCATCGTGAACa > 1:551631/1‑76 (MQ=255) aCGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAAt < 1:1501141/76‑1 (MQ=255) gcgATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAg < 2:1101342/76‑1 (MQ=255) gcgATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAg < 2:734557/76‑1 (MQ=255) cgATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGCTTACCATCGTGGACAATaag > 2:364764/1‑74 (MQ=255) cgATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAAAAATAgg > 2:1866296/1‑76 (MQ=255) gATTCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGa < 1:200436/76‑1 (MQ=255) gATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGCTTACCATCCTGGACAAAAggg > 2:1367925/1‑75 (MQ=255) aTGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGAc > 1:1258683/1‑76 (MQ=255) tGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACg < 2:1491824/76‑1 (MQ=255) cggcAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTa < 2:455118/76‑1 (MQ=255) ggcAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTACCATCGTGGACAATAGGACGCCGTAc > 2:254818/1‑76 (MQ=255) cAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTActg > 2:1336817/1‑76 (MQ=255) ggTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTActgct < 1:959635/76‑1 (MQ=255) tACACGGTTCTCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGa > 1:39092/1‑76 (MQ=255) tACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGa > 2:1597560/1‑76 (MQ=255) tACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTACCATCGTGAACAAAAGGACGCCGTACTGCTGGa > 1:32428/1‑76 (MQ=255) acagGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGaa < 2:1617910/76‑1 (MQ=255) cacGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTActgctg < 2:790215/72‑1 (MQ=255) cacGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTACCATCGTGAACAATAGGACGCCGGActgctg > 1:790215/1‑72 (MQ=255) cacGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGcaa > 2:578421/1‑76 (MQ=255) acGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCCTCCGTAACAATTAGAAGCCGTACCGCTGGCAAt > 1:1194131/1‑76 (MQ=255) ggTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAAttt > 1:1611855/1‑76 (MQ=255) tATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAATTTCtt > 1:506953/1‑76 (MQ=255) tCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTAc < 1:345780/76‑1 (MQ=255) aCTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTTTTACCt > 1:1513989/1‑76 (MQ=255) cTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCtt < 2:372922/76‑1 (MQ=255) tACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCttt < 2:1805506/76‑1 (MQ=255) tACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCttt < 1:1731506/76‑1 (MQ=255) aCGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTa < 2:614552/76‑1 (MQ=255) aCGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTa < 1:1469966/76‑1 (MQ=255) aCGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTa < 2:829356/76‑1 (MQ=255) cGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTAt < 2:1371786/76‑1 (MQ=255) cGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTAt < 1:1490680/76‑1 (MQ=255) tAGCCAATACTGGCTGGGGGATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAATTTCTTACCTTTTTGat > 2:116821/1‑75 (MQ=255) tAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAATTTCTTACCTTTATGAc > 2:1410013/1‑76 (MQ=255) aGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACt < 1:654282/76‑1 (MQ=255) cAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTccc < 1:1369146/76‑1 (MQ=255) aaTACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTcccc < 1:53419/76‑1 (MQ=255) tACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTcccccc > 2:20300/1‑74 (MQ=255) tggctggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGcgcc > 2:1667167/1‑76 (MQ=255) ctggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCgccgcc < 1:666771/76‑1 (MQ=255) ctggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCgccgcc > 2:1550646/1‑76 (MQ=255) tggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGccgccg > 1:774266/1‑76 (MQ=255) ggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTCCCTTTATGACTCCCCGCGCCGCCGt < 2:1790758/76‑1 (MQ=255) tGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAt < 1:1088083/76‑1 (MQ=255) tGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAt < 2:53514/76‑1 (MQ=255) tGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAt < 2:32428/76‑1 (MQ=255) tGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAt > 2:1462956/1‑76 (MQ=255) | GAAAAGAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAT > REL606/716757‑716907 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |