| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_vch_tpiA | 4,314,263 | 0 | A | T | 20.8% | 38.7 / 10.2 | 24 | intergenic (+36/+111) | yjdP/phnP | putative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase |
| Reads supporting (aligned to +/- strand): ref base A (5/14); new base T (3/2); total (8/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAC > NC_000913_3_vch_tpiA/4314116‑4314409 | ccAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAAt < 1:314208/149‑1 (MQ=255) gTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATggg < 1:22778/149‑1 (MQ=255) gATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACCGGAGGAGGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgcgagaggg > 2:532199/1‑149 (MQ=255) agATCGTCACCGCCAGCTTGATTAGGAGGGTCGGCAACTGGAGGATGAAGAGTGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTc < 2:107921/148‑1 (MQ=255) ccagcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgagg < 1:286765/122‑1 (MQ=255) ccagcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgagg > 2:286765/1‑122 (MQ=255) ccagcaccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATc > 2:204592/1‑149 (MQ=255) gcgccagcTTGATCAGGCGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCTATGAGGTTGGGTATGGAGCAGGCATCGTCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAg < 1:208993/149‑1 (MQ=255) gcgccagcGTGATCAGGAGCAGCGGCAACGGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGATGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAg < 2:473975/149‑1 (MQ=255) cagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGcaca > 2:235823/1‑149 (MQ=255) cTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATc < 2:186090/149‑1 (MQ=255) tGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga < 1:438289/149‑1 (MQ=255) gaggAAGGGGGGAGAATGGAAGCGGAGTATGGGCTATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAggg < 2:423341/146‑1 (MQ=255) agagCGGAGAATGGAAGATGGGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGATCGCAATATTCAATATCGAGTACAATCGGTCCACTCGCCCCTCTGGGGAGAGGGTTAgg < 2:4791/149‑1 (MQ=255) gtgaATGGAAGAGTAGTATGGTGTATGAGGTTGTGTATGGAGCAGGCATCGCCCCTCACCCCAATTGGGAGTTTGTTTTTTCGTGCGCAAGATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg < 2:14986/147‑1 (MQ=21) gagaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGAAGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGACGGGGTTAGGTTGAgg > 1:221522/1‑149 (MQ=255) ggggGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCTCAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc < 1:310414/149‑1 (MQ=255) aTCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatc < 1:97246/149‑1 (MQ=255) ccTCTCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTccc < 1:401436/116‑1 (MQ=255) ccTCTCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTccc > 2:401436/1‑116 (MQ=38) ctctCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGATCCCCACGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCGAGCCACCCCAATCTCCATCCCATCaa > 1:160040/1‑149 (MQ=38) tCACCCCAATGGGGAGAGGGAGAACACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCaaa < 1:509382/149‑1 (MQ=255) aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAcc > 2:515382/1‑149 (MQ=255) cAACGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAc < 1:204592/149‑1 (MQ=255) | CCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAC > NC_000913_3_vch_tpiA/4314116‑4314409 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |