Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,992,881 | 0 | A | C | 37.1% | 41.0 / 11.2 | 35 | F238V (TTC→GTC) | PP_3517 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (14/8); new base C (0/13); total (14/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.91e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATCCTGGCCGAACAGCGGGCC > NC_002947/3992747‑3993013 | gAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAg > 1:476842/1‑136 (MQ=255) aGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGcc > 2:794063/1‑136 (MQ=255) gCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCTTGAAGCTCATCACCCGGGTGTCGTCGTAGTTGGCTGG‑CCAGGCGACCATGGCGAAGCCa < 1:204408/136‑1 (MQ=255) cGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCAc > 2:444825/1‑136 (MQ=255) ttttCATCTCCTGCCCCCCCCTGGCGCCCTGCGGCCCGAGGTCTTTCTCCAACACTTGGCCGTCAGGGCTGATCATGAACCACATCACCCCGGTGGCGTCGAACTGGGCCGG‑CCCGGCGACCAGGGCTACGCCACg < 2:584286/135‑1 (MQ=255) tctcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCAt > 1:625030/1‑136 (MQ=255) tctcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCAt > 1:578948/1‑136 (MQ=255) ccAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCCGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGcc < 2:476842/136‑1 (MQ=255) cAGACTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGGTGCCg > 2:445895/1‑136 (MQ=255) ccTGTCGCCTTGCGGCCGGTGGTCTTTCTCTAACACTTGGCGGTCTTGGCTGATCTTGAAGCTCATCACCCGGGTGGGGTGTTACTTGGCGGG‑CCAGGCTACCTTGGCGACCCCACGGCTCATCTGGTTGCCGATc < 2:457670/136‑1 (MQ=255) ccTTTCTCTAACACTTGCCCTTCTTGGCTTATAATCACCCTCTCCCCCCGGGTGTCGGCTTAGTTGGCCCG‑CCCGGCTACCTGGGCCAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCg < 2:814367/135‑1 (MQ=255) acacTTGGCCGTCATGGCTGACCAGCACGCCCACCCCCCCGGGGCCGTCGAACTTGGCCGG‑CCCGGCGACCTGGGCCACGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATgg < 2:525919/136‑1 (MQ=255) cGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGc > 2:415253/1‑136 (MQ=255) caCCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATc > 1:550298/1‑70 (MQ=255) caCCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATc < 2:550298/70‑1 (MQ=255) aCCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCAt > 2:279622/1‑136 (MQ=255) ggTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGACCCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAg < 1:368984/136‑1 (MQ=255) tgtCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCAGTGCATACAGGATGCGGAAGTGGTAGCCATACCAgcc > 1:649173/1‑134 (MQ=255) tcgtcgTACTTGGCCGG‑CCAGGCGACCTTGGCCACGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGt < 1:250007/136‑1 (MQ=255) gggtcgtAC‑TGGCCGGCCCCGGCGACCAGGGCTACGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGt < 1:833035/134‑1 (MQ=255) gcgtcgtACTGGGCCGG‑CCCGGCGACCTTGGCGACGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGt < 1:245471/135‑1 (MQ=255) cgtcgtACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTc > 2:284478/1‑136 (MQ=255) gtACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCg > 2:53939/1‑136 (MQ=255) tACTTGGCAGG‑CCAGGCGACCATGGCGAAGCCACTGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGTGCCCAGTGCATCCATGATGCGGAAGTGTTAGCCATACCAGGCGTCACCGa > 2:208674/1‑136 (MQ=255) ttGGCCGC‑CCCGGCCACCTGGGCCACGCCACGGCTCATCTGGTTGCTGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACca < 1:453410/136‑1 (MQ=255) ggccgg‑ccCGGCGACCATGGCGACCCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACcaca < 2:578948/136‑1 (MQ=255) ggccgg‑ccCGGCGACCATGCCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACcaca < 2:193401/136‑1 (MQ=255) ggccgc‑ccCGGCCACCATGCCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACcaca < 2:625030/136‑1 (MQ=255) ccgc‑ccAGGCGACCAGGCCCACCCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATc < 2:64911/136‑1 (MQ=255) cgg‑ccAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATcc > 2:554391/1‑136 (MQ=255) ctacgttGGCTACGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATCCTGGCCGaa < 1:60353/130‑1 (MQ=255) ccctgggCGACCCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATCCTGGCCGAACa < 1:75364/131‑1 (MQ=255) gCGACGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATCCTGGCCGAACAGCGGGc < 1:444825/136‑1 (MQ=255) gCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATCCTGGCCGAACAGCGGGc < 1:415253/136‑1 (MQ=255) cGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATCCTGGCCGAACAGCGGGcc < 1:236215/136‑1 (MQ=255) | GAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGG‑CCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGCGGAAGTGGTAGCCATACCAGGCGTCACCGACCACATCCTGGCCGAACAGCGGGCC > NC_002947/3992747‑3993013 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |