| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,282,813 | 0 | T | G | 43.8% | 24.2 / 61.6 | 48 | intergenic (‑162/+67) | PP_2010/PP_2011 | cytochrome b561/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (12/15); new base G (9/12); total (21/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCA > NC_002947/2282672‑2282962 | tCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc < 1:1269862/147‑3 (MQ=255)tCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc > 2:1269862/1‑145 (MQ=255)tCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttccaag > 2:423838/1‑145 (MQ=37) gAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCt < 1:1109744/150‑1 (MQ=255) ggCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGc < 2:448169/150‑1 (MQ=255) ggCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTG‑CTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa < 1:291631/150‑1 (MQ=18) gcgcTCTTGGTGGCGCTTATTTTTCTTTTTAATGGCGACCGGAAGCAGTTTTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGCCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa < 2:245751/150‑1 (MQ=37) gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGACCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa < 2:637617/150‑1 (MQ=37) cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCt > 2:971262/1‑144 (MQ=255) cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCt < 1:971262/144‑1 (MQ=255) cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg < 1:423838/150‑1 (MQ=25) ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGGAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt < 2:618493/150‑1 (MQ=255) ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt < 2:889974/150‑1 (MQ=255) ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt < 2:256061/150‑1 (MQ=37) ttGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGa < 1:876202/150‑1 (MQ=255) ttGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGa < 1:397781/150‑1 (MQ=255) ggAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCggg < 1:279051/150‑1 (MQ=255) tCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCttcc > 1:969164/1‑114 (MQ=255) tCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCttcc < 2:969164/116‑3 (MQ=255) ttttAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc < 1:630503/113‑3 (MQ=25) ttttAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc > 2:630503/1‑111 (MQ=25) gcCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg < 2:380795/92‑1 (MQ=255) gcCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg > 1:380795/1‑92 (MQ=255) gcCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt > 1:597258/1‑97 (MQ=18) gcCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt < 2:597258/97‑1 (MQ=18) aaaaTGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACtgtg < 2:359503/79‑1 (MQ=14) aaTGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCt < 1:42849/74‑1 (MQ=255) aaTGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGATCCCACAGTCCAACTGCACAGTCTTGGAAGCCt > 2:42849/1‑74 (MQ=255) cGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCATGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTTcctcct > 2:190797/1‑150 (MQ=255) cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAg > 1:1086218/1‑100 (MQ=255) cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAg < 2:1086218/100‑1 (MQ=255) acaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGAGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACccacca < 1:420915/122‑1 (MQ=255) acaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGAGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACccacca > 2:420915/1‑122 (MQ=255) acaGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGc > 2:629010/1‑150 (MQ=37) tCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCa > 2:1140561/1‑150 (MQ=255) tCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCa > 2:1168556/1‑150 (MQ=21) ccAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTTCCTCCTGCCAGGCACCAt > 1:660725/1‑150 (MQ=255) cAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt < 1:660535/130‑1 (MQ=25) cAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt > 2:660535/1‑130 (MQ=25) aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc > 2:395059/1‑149 (MQ=255) aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc < 1:395059/149‑1 (MQ=255) aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg > 1:419385/1‑150 (MQ=37) tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa > 2:485245/1‑150 (MQ=255) tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACa < 1:819429/74‑1 (MQ=11) tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACa > 2:819429/1‑74 (MQ=11) gCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg < 1:1055908/137‑1 (MQ=25) cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGc > 1:769144/1‑58 (MQ=255) cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGc < 2:769144/58‑1 (MQ=255) tCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTTGTCGGTGACCATAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCa > 2:985239/1‑150 (MQ=38) tCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCa > 2:172491/1‑150 (MQ=255) | TCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCA > NC_002947/2282672‑2282962 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |