| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_017635 | 123,066 | 0 | A | T | 35.9% | 52.9 / 19.0 | 39 | L15L (CTT→CTA) | aroP | aromatic amino acid transporter AroP |
| Reads supporting (aligned to +/- strand): ref base A (9/16); new base T (14/0); total (23/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.56e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGA > NC_017635/123017‑123112 | cGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTATTTa > 2:1770343/1‑50 (MQ=255) ggTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTaa < 2:1483628/50‑1 (MQ=255) ggTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTa > 2:1051589/1‑49 (MQ=255) gTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAg > 1:404593/1‑50 (MQ=255) gTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAg < 2:1732856/50‑1 (MQ=255) tCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAgg < 1:2031142/50‑1 (MQ=255) tCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGATTTTAAgg > 2:271805/1‑50 (MQ=255) cccTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAGGc > 1:393964/1‑50 (MQ=255) cccTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAAGc > 1:1381236/1‑50 (MQ=255) ccTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAgggc > 2:272799/1‑48 (MQ=37) ccTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAAGcc > 1:1131368/1‑50 (MQ=37) ccTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGcc < 2:1526092/50‑1 (MQ=255) cTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAGGCcg > 1:1940260/1‑50 (MQ=255) cTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAAGccc > 2:298067/1‑49 (MQ=25) gcgcCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAATCCGCGGtt > 2:132022/1‑50 (MQ=17) cgcCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCttt < 1:563008/50‑1 (MQ=255) gcCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTa < 1:2080691/50‑1 (MQ=255) ccACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAAGCCGCCCttttg > 1:1453043/1‑48 (MQ=37) cACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCCCCTTaac > 1:1664285/1‑48 (MQ=37) aCCCAGCGCGATAAGCTGAATATGGAGGGTGTTTAAAGCC‑CGCTATAgct > 2:1094267/1‑50 (MQ=255) aCCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGGCCCGCTTTAgct > 2:1828228/1‑50 (MQ=38) aCCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAgct > 1:401688/1‑50 (MQ=255) aCCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAgct < 2:1589729/50‑1 (MQ=255) cAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAGGCCGCGCTTTTgctgct > 2:1856490/1‑50 (MQ=38) cAGCGCGATAAGCTGAATATGG‑CGGTTTTTTAGGCCGCGCCTTAactgct > 2:448765/1‑50 (MQ=37) cAGCGCGATAAGCTGAATATGG‑CGGATTTTTAAGCCGCGCCTTAtctgct > 2:1402650/1‑50 (MQ=25) aGCGCGATAAGCTGAATATGG‑CGGTTTTTTAAGCCGCGCTTTAGCTGCTc > 2:1326240/1‑50 (MQ=38) aGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTc < 2:1929818/50‑1 (MQ=255) aGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTc < 2:1852443/50‑1 (MQ=255) gcgcGATAAGCTGAATATGG‑CGGTTTTTTAGACCGCGCTTTATATGCTCg > 2:385202/1‑50 (MQ=25) cgcgATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTCGc > 1:1628062/1‑50 (MQ=255) cgATAAGCTGAATATGG‑CGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCg > 1:1667456/1‑50 (MQ=255) gATAAGCTGAATATGG‑CGGTTTTTTAGGCCGCGCCTTAGCTGCTCGCCGt > 2:1405786/1‑50 (MQ=38) gATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTcgtcgt < 1:1206208/50‑1 (MQ=255) aGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCt < 1:1955643/49‑1 (MQ=255) aaTATGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTg < 2:134337/48‑1 (MQ=255) atatGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGAcc < 1:1099977/50‑1 (MQ=255) tatGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCt < 1:1312483/50‑1 (MQ=255) tGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTc < 2:1164392/50‑1 (MQ=255) cGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCcat < 1:2033581/50‑1 (MQ=255) gTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCcatca < 2:1821994/50‑1 (MQ=255) ttAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGa > 2:1101978/1‑49 (MQ=255) | CGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGG‑CGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGA > NC_017635/123017‑123112 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |