Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,327,127 | 0 | C | T | 38.9% | 27.3 / 35.6 | 37 | intergenic (+97/‑196) | PP_16SE/PP_23SE | 16S ribosomal RNA/23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base C (11/11); new base T (7/7); total (19/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAG > NC_002947/1326986‑1327231 | tAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATAAAATATTTAATAGATGATATAAGCCTCCTGATGATCGCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATg > 1:11745/1‑150 (MQ=11) aCAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGtt > 1:515821/1‑150 (MQ=32) aaGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGa < 2:325586/150‑1 (MQ=32) aaGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGa > 1:255011/1‑150 (MQ=32) cGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGc < 2:405378/150‑1 (MQ=32) atgggAACCTGCGGCTGGATCACCTCCTTAATCGACGGCATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCCTTGTCGAAGAGGTTTTTCCTGGNCAGNGCNTGGAAATGAGCANTCGCATCGAATGTTGATTTCTGGCttt < 1:584713/148‑1 (MQ=11) ccTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTc > 2:79726/1‑147 (MQ=17) ccTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTc < 1:79726/147‑1 (MQ=17) tGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCATATc < 2:310255/150‑1 (MQ=17) tGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCtt > 1:111293/1‑150 (MQ=32) aTCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaa > 1:573222/1‑150 (MQ=32) ctcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAAtt > 2:572261/1‑150 (MQ=32) tcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc < 1:572261/150‑1 (MQ=32) ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTAt < 2:363637/103‑1 (MQ=18) ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCgg > 2:283636/1‑150 (MQ=17) ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCgg > 1:308982/1‑150 (MQ=17) aaTCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGata < 1:283636/150‑1 (MQ=17) tCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtg > 1:550767/1‑150 (MQ=17) tCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtg > 2:258272/1‑150 (MQ=32) gCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAg < 2:515821/150‑1 (MQ=32) ccTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGa < 2:255011/150‑1 (MQ=32) tgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACt > 1:141061/1‑150 (MQ=17) tgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACt < 1:616780/150‑1 (MQ=17) cTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCac < 2:127060/73‑2 (MQ=25) cTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCac > 1:127060/1‑72 (MQ=25) cTCCCACACGAATTGCTTGATTCATTGTCGAAGACGATTTTCATGGTCACAACCTGGAACTGCGCATTGGGATCGAATGGTTATTTCTGGCTTTTGGCATATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGt > 2:70241/1‑150 (MQ=11) ccACACGCATTGCTTGATTACTTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc > 2:6716/1‑150 (MQ=12) ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc > 1:539042/1‑150 (MQ=17) cGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgc < 2:308982/150‑1 (MQ=17) aaTTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgctg > 2:274466/1‑150 (MQ=32) aTTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATATACTGAACACCAGTTTCACTGCtgg < 2:11745/150‑1 (MQ=17) ttGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt < 2:573222/150‑1 (MQ=32) ttGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt > 2:173261/1‑150 (MQ=32) ttGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagata < 2:210214/110‑1 (MQ=18) ttGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagata > 1:210214/1‑110 (MQ=18) tGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg < 2:341674/150‑1 (MQ=17) tGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg < 2:111293/150‑1 (MQ=32) | TAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAG > NC_002947/1326986‑1327231 |
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: ‑ |