Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,282,779 | 0 | C | A | 28.3% | 55.7 / 29.9 | 46 | intergenic (‑171/+58) | PP_2010/PP_2011 | cytochrome b561/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (15/18); new base A (7/6); total (22/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.46e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAC > NC_002947/2282632‑2282920 | gAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGcc > 8:262250/1‑149 (MQ=255) tGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgt < 7:83572/149‑1 (MQ=255) cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCTCCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa < 2:336446/149‑1 (MQ=37) cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTCTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc < 1:272746/149‑1 (MQ=255) cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa < 6:157180/149‑1 (MQ=255) cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCT‑TCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa < 5:68998/149‑1 (MQ=25) ttGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt < 4:116239/147‑1 (MQ=255) ttGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt < 8:79794/147‑1 (MQ=255) gCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTg < 7:293568/146‑1 (MQ=255) gcgcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCTCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACAGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg < 2:302284/149‑1 (MQ=255) gAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCggg > 3:218562/1‑149 (MQ=255) aTTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTGTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGc > 6:200292/1‑149 (MQ=255) tAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTg < 6:155326/125‑1 (MQ=255) tAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTg > 5:155326/1‑125 (MQ=255) gCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAgg < 4:218562/149‑1 (MQ=255) gAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCg > 1:132159/1‑147 (MQ=255) gAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCg > 4:103696/1‑147 (MQ=255) tGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACAt < 6:137189/148‑1 (MQ=255) gttaaaatTGTGTTCCCTCCAACGGGCCTCT‑TCGCTGGAAGGCCCGCTCCCAAAGACAAACTGCACTTTTTTGTACGCGTTTGCTTTTGTAGTGGGGTTGGGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGc < 2:132159/145‑1 (MQ=255) tgtgTTGCCAGCCCCCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTCCACCGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGc < 5:200292/149‑1 (MQ=255) tgtgTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCGTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGc > 7:172186/1‑149 (MQ=255) tgtgTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACAGCACAGTCGTGGGAGCCGGGGGAGTTGGACAGTGGGAGCGGGCGTGCCCGCGCAGGGGCCGG‑AGGAGGCGACATCAGGTAACCGGCAGATAggtc > 8:264030/1‑147 (MQ=25) ttGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGGCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAg < 7:262250/148‑1 (MQ=255) cgccaCCGGCCTCT‑TCGCTGGCACGCCCGCACCCACAGCCCACCTGCACAGTCTTGGCAGCCTCTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGtcatca < 6:159652/148‑1 (MQ=37) aCCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcaccca > 6:162329/1‑149 (MQ=255) ggCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGGGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcacccac > 2:121636/1‑147 (MQ=255) tGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCg < 8:158272/88‑1 (MQ=255) tGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCg > 7:158272/1‑88 (MQ=255) cACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGACGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt < 8:172186/149‑1 (MQ=255) aCGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTg > 4:318456/1‑149 (MQ=255) gCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcc > 1:188073/1‑149 (MQ=255) cccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGcc < 7:264030/149‑1 (MQ=255) ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGc > 6:84816/1‑147 (MQ=255) cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg > 3:216449/1‑149 (MQ=255) aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg < 4:163388/141‑1 (MQ=255) aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg > 3:163388/1‑141 (MQ=255) aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc > 7:287232/1‑149 (MQ=255) aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc > 1:217071/1‑149 (MQ=25) aCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg < 3:318456/149‑1 (MQ=255) aCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg > 5:343103/1‑149 (MQ=25) cTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCAACATCAGTTACCGGGCAGAGAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg > 3:306420/1‑149 (MQ=18) tGCACAGGATTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGCGGTGCCCGCGAAGGGGCCCG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCgg > 4:219980/1‑148 (MQ=37) cacaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt > 6:114292/1‑149 (MQ=255) acaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt > 6:263622/1‑148 (MQ=255) aGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAg < 5:84816/149‑1 (MQ=255) ggcttccaaGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAgg < 2:217071/142‑1 (MQ=37) ccaaGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCAACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAc < 4:306420/146‑1 (MQ=25) | GAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCT‑TCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGG‑AGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAC > NC_002947/2282632‑2282920 |
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: ‑ |