| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_017635 | 4,363,963 | 0 | A | G | 46.0% | 11.9 / 121.1 | 74 | V340V (GTA→GTG) | ECW_RS21820 | phage late control D family protein |
| Reads supporting (aligned to +/- strand): ref base A (31/9); new base G (26/8); total (57/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTCTGGCTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCACTAAG > NC_017635/4363914‑4364011 | gTCTGGCTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTa > 1:1135994/1‑50 (MQ=32) ggCTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTGTCAg > 1:243123/1‑50 (MQ=255) ggCTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAg > 2:1208752/1‑50 (MQ=255) gCTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAgg > 2:1079329/1‑50 (MQ=255) cTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGc > 1:997316/1‑50 (MQ=255) cTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGc < 1:1430566/50‑1 (MQ=255) cTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGAATCAGGc > 2:1473836/1‑50 (MQ=255) tACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCt > 1:29946/1‑50 (MQ=255) tACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCt > 2:2034047/1‑50 (MQ=255) tACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCt > 1:498139/1‑50 (MQ=255) cTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCtt > 2:932355/1‑49 (MQ=255) cTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCtt > 2:37392/1‑49 (MQ=255) cTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCttt < 2:1181160/50‑1 (MQ=255) tGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCttt > 2:975515/1‑49 (MQ=255) tGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCttt > 1:1350165/1‑49 (MQ=255) tGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCttt > 2:1376199/1‑49 (MQ=255) ggTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCttt > 1:1352309/1‑48 (MQ=255) tCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAgc > 1:741262/1‑50 (MQ=255) tCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAgc > 1:629434/1‑50 (MQ=255) tCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAgc > 2:444129/1‑50 (MQ=255) cGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAgcg > 2:1473991/1‑50 (MQ=255) cAGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTc < 2:1247478/50‑1 (MQ=255) cAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGt > 2:1564467/1‑49 (MQ=255) aGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCa > 2:1488659/1‑50 (MQ=255) aGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCa > 2:1260215/1‑50 (MQ=255) aGATATTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCa > 1:1289889/1‑50 (MQ=255) gATATTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCa > 1:1232918/1‑49 (MQ=255) atatTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCa > 2:737354/1‑48 (MQ=255) atatTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCa > 2:972872/1‑48 (MQ=255) tatTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAg < 2:2054814/50‑1 (MQ=255) atTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGa > 1:1441445/1‑50 (MQ=255) atTTACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGa > 2:187480/1‑50 (MQ=255) atTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGa > 2:626433/1‑50 (MQ=255) atTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGa > 1:1565004/1‑50 (MQ=255) tttACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGa < 2:616080/49‑1 (MQ=255) tACACGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGAc > 1:1119462/1‑48 (MQ=255) acacGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAg > 2:1061733/1‑50 (MQ=255) acacGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAg > 2:518191/1‑50 (MQ=255) acacGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAg > 2:1512412/1‑50 (MQ=255) cacGGAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGc < 1:1188901/50‑1 (MQ=255) acGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCa < 2:563638/50‑1 (MQ=255) acGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCa < 2:275344/50‑1 (MQ=255) ggAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAgg < 2:1136522/50‑1 (MQ=255) ggAAACGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAgg < 2:355561/50‑1 (MQ=255) gAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGa > 2:795480/1‑50 (MQ=255) aaaCGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGAc > 1:1043457/1‑50 (MQ=255) aaaCGCCGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGAc > 1:219885/1‑50 (MQ=255) aaCGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACt > 1:1418949/1‑50 (MQ=255) aaCGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACt > 2:55732/1‑50 (MQ=255) aaCGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACt < 2:310728/50‑1 (MQ=255) aaCGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACt > 2:1102959/1‑50 (MQ=255) cGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGaaa < 2:1745723/50‑3 (MQ=11) cGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGAc > 1:1326114/1‑49 (MQ=255) cGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGAc > 2:870081/1‑49 (MQ=255) cGGTCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGAc > 1:1297495/1‑49 (MQ=255) ggTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACaa > 2:1326247/1‑50 (MQ=255) ggTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACaa > 2:1850950/1‑50 (MQ=255) ggTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACaa > 1:1055694/1‑50 (MQ=255) gTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAAt > 2:1192274/1‑50 (MQ=255) gTCAAAGTATAAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAAt < 1:356540/50‑1 (MQ=255) tCAAAGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATc < 1:1242714/50‑1 (MQ=255) tCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATc < 1:1704634/50‑1 (MQ=255) cAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCa > 1:2047994/1‑50 (MQ=255) cAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCa > 1:1698688/1‑50 (MQ=255) aaaGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCAc > 2:472918/1‑50 (MQ=255) aaaGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCa < 1:662725/49‑1 (MQ=255) aaaGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCAc > 1:59479/1‑50 (MQ=255) aaaGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCAc > 2:165784/1‑50 (MQ=255) aaaGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCAc > 2:72439/1‑50 (MQ=255) aaGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCAc > 1:2059824/1‑49 (MQ=255) aaGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCAc > 1:878479/1‑49 (MQ=255) aaGTGTCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCACt < 1:1613314/50‑1 (MQ=255) aaGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCAc > 2:830958/1‑49 (MQ=255) gTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCACt > 2:1792016/1‑48 (MQ=38) tATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCACTAAg > 1:1200711/1‑50 (MQ=35) | GTCTGGCTACTGGTCGGGCAGATATTTACACGGAAACGCCGGTCAAAGTATCAGGCTTTAAGCGCGTCATAGACGAGCAGGACTGGACAATCACTAAG > NC_017635/4363914‑4364011 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |