| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,634,961 | 0 | A | C | 22.9% | 85.5 / 21.4 | 35 | I296S (ATC→AGC) | stfQ | Qin prophage; putative side tail fibre assembly protein |
| Reads supporting (aligned to +/- strand): ref base A (16/11); new base C (0/8); total (16/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.76e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01 | |||||||||||
CATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAG > NC_000913/1634871‑1635058 | cATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgc < 1:1196553/100‑6 (MQ=1)cATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgc < 1:2568546/100‑6 (MQ=1) aTTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgcg < 2:2339454/100‑7 (MQ=1) aTTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgcg < 2:2339395/100‑7 (MQ=1) aTTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgcg < 1:2788028/100‑7 (MQ=1) aTTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgcg < 1:2788642/100‑7 (MQ=1) aTTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgcg < 1:551638/100‑7 (MQ=1) ttCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCATCACCGTTAACGCTGCTGgtaacgcgg < 2:2646583/100‑8 (MQ=1) atacAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGc > 1:1591470/4‑100 (MQ=1) gtgtTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCAt > 1:879965/1‑100 (MQ=35) tgtTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATg > 2:1968450/1‑100 (MQ=35) tttCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATgcgc > 1:553892/1‑100 (MQ=40) tttCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATgcgc > 1:2820645/1‑100 (MQ=40) ttCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATgcgcg > 1:1490077/1‑100 (MQ=40) cgcgTTACCAGCAGTGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCt > 2:1520579/1‑100 (MQ=255) gcgTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTg > 1:1319113/1‑100 (MQ=255) tACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCTACAGTATGTGCATGCGCGCCTgcgc > 1:962867/1‑99 (MQ=255) tACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCt < 2:1591453/100‑1 (MQ=255) ccAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTg > 1:1360007/1‑100 (MQ=255) cagcagCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgc > 2:3171304/1‑100 (MQ=255) cagcagCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTg < 2:2830827/99‑1 (MQ=255) agcagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgca < 1:3226618/100‑1 (MQ=255) agcagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgca > 2:1538000/1‑100 (MQ=255) agcagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgca > 2:2249541/1‑100 (MQ=255) agcagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgca > 2:2761540/1‑100 (MQ=255) gcagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcag < 1:856123/100‑1 (MQ=255) gcagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcag < 1:2837429/100‑1 (MQ=255) cagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcagc < 2:879954/100‑1 (MQ=255) cagcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcagc < 1:905782/100‑1 (MQ=255) cGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGt > 2:2329447/1‑100 (MQ=255) aaCGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCg < 1:137912/100‑1 (MQ=255) aCGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCgg < 2:2274211/100‑1 (MQ=255) aCGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCgg < 1:3250093/100‑1 (MQ=255) cGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGa < 2:2975346/100‑1 (MQ=255) tGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAg > 2:590185/1‑100 (MQ=255) | CATTCTGAATGCCATTATGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAG > NC_000913/1634871‑1635058 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |