Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 699,722 | 0 | A | G | 37.7% | 37.9 / 51.4 | 55 | intergenic (+115/‑129) | PP_t06/PP_23SD | tRNA‑Ala/23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (13/20); new base G (12/8); total (26/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.68e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTT > NC_002947/699580‑699803 | aGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGaa < 6:99871/148‑1 (MQ=37) ggTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAaca > 2:64289/1‑147 (MQ=37) ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCa > 2:136133/1‑148 (MQ=37) ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCa < 8:108752/148‑1 (MQ=37) cccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCa > 5:138760/1‑148 (MQ=37) cccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa < 5:10983/148‑1 (MQ=37) ttGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgct > 6:114965/1‑148 (MQ=37) cTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggt < 4:62727/147‑1 (MQ=37) cTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt < 7:215/147‑1 (MQ=37) tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggtgg > 8:37446/1‑148 (MQ=37) tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggtgg > 1:100529/1‑148 (MQ=37) accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGa > 7:95601/1‑146 (MQ=37) accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATc < 6:61921/148‑1 (MQ=37) ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa > 8:137057/1‑148 (MQ=37) acTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg < 3:49854/148‑1 (MQ=37) cGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTa < 8:2857/146‑1 (MQ=37) gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg < 8:77641/147‑1 (MQ=37) gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg > 6:85218/1‑147 (MQ=37) gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg < 8:60053/147‑1 (MQ=37) ccGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGa > 6:69895/1‑136 (MQ=37) cGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGt > 3:39262/1‑147 (MQ=37) cGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt < 5:69895/136‑1 (MQ=37) aCGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTaaaa < 4:39262/148‑1 (MQ=37) aGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtg > 2:76885/1‑148 (MQ=37) gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt < 3:97022/148‑1 (MQ=37) gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt < 4:147634/148‑1 (MQ=37) gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt < 2:2251/148‑1 (MQ=37) gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt < 3:10885/148‑1 (MQ=37) tACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa < 6:30443/148‑1 (MQ=35) tACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGGTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa < 6:6901/148‑1 (MQ=25) tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTg > 7:104074/1‑148 (MQ=21) tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTg > 4:113179/1‑148 (MQ=21) tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTg > 4:112932/1‑148 (MQ=34) cAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc < 4:138165/148‑1 (MQ=21) cAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc < 2:134622/148‑1 (MQ=34) cAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc < 4:36857/148‑1 (MQ=34) aGAACTTAGAAATGATCATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt > 5:157425/1‑148 (MQ=21) aGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt > 4:69451/1‑148 (MQ=21) aaCTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg < 1:171783/148‑1 (MQ=34) ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGTCTGAACAACAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg > 6:150640/1‑145 (MQ=21) ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa < 7:94770/148‑1 (MQ=34) ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGGGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa < 7:54980/148‑1 (MQ=21) aaaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag < 1:76885/147‑1 (MQ=21) aaaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag < 7:137057/147‑1 (MQ=34) aaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAGCACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCa > 6:27426/1‑148 (MQ=21) aaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCa > 3:103168/1‑148 (MQ=34) attaaatttC‑ATCTCGAATGTTGATTTCTTTCTTTTGTCAGATCGTTATTTAAAAATTCTGTTATGTGATATATATAGTCTGAACACCATTTTCACTGCTGTTGTATCAGTCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCa < 4:87995/140‑1 (MQ=11) aTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa > 6:29406/1‑148 (MQ=21) aTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa > 6:3401/1‑148 (MQ=34) tcaaCATTCGAGCT‑GAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc > 1:79197/3‑148 (MQ=2) aTTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTAc > 8:17839/1‑147 (MQ=21) tC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGaa > 3:102806/1‑148 (MQ=18) tC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGaa > 2:6679/1‑148 (MQ=33) tC‑AGCTCGAATGTTGATTTCTGGCGTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGaa < 2:15/148‑1 (MQ=18) cTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtttt < 4:102806/148‑1 (MQ=11) | AGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTT > NC_002947/699580‑699803 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |