| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,842,941 | 0 | C | . | 51.8% | 87.0 / 129.5 | 56 | intergenic (+40/‑257) | ccoP‑II/PP_4259 | cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein |
| Reads supporting (aligned to +/- strand): ref base C (14/13); new base . (13/16); total (27/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCA > NC_002947/4842796‑4843089 | gCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGAt < 1:336645/150‑1 (MQ=255) aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGa < 2:1013801/148‑1 (MQ=255) aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGa > 1:1013801/1‑148 (MQ=255) aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATcc < 1:174026/150‑1 (MQ=255) aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATcc > 2:467791/1‑150 (MQ=255) aTGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt > 2:293105/1‑150 (MQ=255) cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCTTACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCtgtgta < 1:414994/150‑6 (MQ=38) tgcacAGGAACAACTGCAGGGCAACGCCAAGTTGCACTTTCTGGCGTCTTAAGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAg < 2:744645/146‑1 (MQ=39) ggCCCAGGAACAACTGCAGGTCAACGACAAGGTGCACGTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGCCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAg < 1:274264/150‑1 (MQ=25) cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGa < 1:293105/150‑1 (MQ=255) cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGa < 1:1876235/150‑1 (MQ=255) ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAg < 1:844486/150‑1 (MQ=255) ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAg > 2:844486/1‑150 (MQ=255) ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAg < 1:467791/150‑1 (MQ=38) ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGGAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg < 2:1477493/150‑1 (MQ=255) ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg < 2:1711175/150‑1 (MQ=255) ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg < 1:284969/150‑1 (MQ=255) ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg < 2:674853/150‑1 (MQ=255) ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg < 2:1060039/150‑1 (MQ=255) ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGAAACCTCTGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg < 1:1573569/150‑1 (MQ=37) aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGc > 1:407039/1‑146 (MQ=255) aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGc < 2:407039/146‑1 (MQ=255) aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCggg < 2:553372/135‑1 (MQ=255) aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCggg > 1:553372/1‑135 (MQ=255) ggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAg > 1:205436/1‑105 (MQ=37) ggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAg < 2:205436/105‑1 (MQ=37) cGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg > 2:622177/1‑105 (MQ=255) cGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg < 1:622177/105‑1 (MQ=255) gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg < 1:1649953/104‑1 (MQ=25) gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg > 2:1649953/1‑104 (MQ=25) tCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCtgtgt > 1:1494993/1‑91 (MQ=255) tCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCtgtgt < 2:1494993/95‑5 (MQ=255) tCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc > 1:664503/1‑80 (MQ=255) tCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc < 2:664503/80‑1 (MQ=255) gCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg > 2:1196766/1‑83 (MQ=18) gCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTAGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg < 1:1196766/83‑1 (MQ=18) gCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg < 1:1070347/74‑1 (MQ=17) gCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg > 2:1070347/1‑74 (MQ=17) cAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAg > 1:395215/1‑69 (MQ=255) cAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAg < 2:395215/69‑1 (MQ=255) tctTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGt < 1:166178/96‑1 (MQ=255) tctTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGt > 2:166178/1‑96 (MQ=255) tctTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGCTCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAAcgc > 1:903867/1‑150 (MQ=255) gcgGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGCaga < 1:369965/56‑1 (MQ=255) gcgGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAaga > 2:369965/1‑56 (MQ=255) aCCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGc > 2:1055905/1‑150 (MQ=255) cccGCTCCTACACAGGATCCT‑TGTAGGAACGGCTCTAGCCGCGAAGAGGCCAGTACTGGCACAACATCTCTCGCAGGCAAGTGACCTGGATCAATTGACACCAGCCATAACACGTTTCATACCACGAACCCAACGCGACGAAAGGTGGCa > 1:733882/1‑150 (MQ=255) ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg > 2:1580884/1‑150 (MQ=255) ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg > 1:1245659/1‑150 (MQ=255) ccGCTCCTACA‑AGGGTCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg > 1:1901359/1‑150 (MQ=255) ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACACGACTAAAGGTCGCAg > 2:378389/1‑150 (MQ=255) gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTaaa < 2:1290905/140‑1 (MQ=255) gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTaaa > 1:1290905/1‑140 (MQ=255) gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg > 2:467052/1‑150 (MQ=255) gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg > 1:1664959/1‑150 (MQ=255) tCCTACA‑AGGATCCTGTGTAGGAGCGGGATTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg > 1:1306172/1‑150 (MQ=255) cacaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCa > 2:340618/1‑69 (MQ=255) cacaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCa < 1:340618/69‑1 (MQ=255) acaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa < 2:1143550/55‑1 (MQ=255) acaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa > 1:1143550/1‑55 (MQ=255) acaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa > 2:638308/1‑150 (MQ=255) | GCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCA > NC_002947/4842796‑4843089 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |