Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 699,725 | 0 | A | G | 39.6% | 26.4 / 41.8 | 53 | intergenic (+115/‑129) | PP_t06/PP_23SD | tRNA‑Ala/23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (19/13); new base G (14/7); total (33/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.73e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCG > NC_002947/699580‑699808 | gTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGa > 1:1149436/1‑150 (MQ=255) tCAGTGGTTGGAGCCCGCTTGGCTCAACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAACTGACCATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGAAATGTGTTAGTTCTAGACTGaa > 2:199496/1‑150 (MQ=11) aGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAaca > 2:650245/1‑150 (MQ=38) gCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTTAacac > 2:1123112/1‑150 (MQ=38) gATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTc > 1:193055/1‑150 (MQ=255) cccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAct < 1:1175213/150‑1 (MQ=39) tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggt < 2:708781/150‑1 (MQ=255) ggCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg > 2:264613/1‑150 (MQ=38) gCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg > 2:82191/1‑150 (MQ=38) accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAg > 2:167408/1‑150 (MQ=38) ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAgg > 1:1020076/1‑150 (MQ=38) ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg > 1:338880/1‑150 (MQ=39) ccacTTTCGCCGTACGCAGTTACTTCTCAGATCTTAGAAATGAACATTC‑AGCTTGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATACGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg < 1:272771/150‑1 (MQ=255) cacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGc < 1:673361/150‑1 (MQ=38) cTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTa < 1:706878/150‑1 (MQ=255) tttGGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAAAATTA‑AGATTGAATGTTGATTTCTGGCTTTGGGCAGATCGTTGTTTAAGAATTCGTATGAGTTATAGATATTGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa > 1:672457/1‑150 (MQ=11) tttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa > 1:767711/1‑150 (MQ=255) gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTa < 1:698303/150‑1 (MQ=39) tACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAt > 2:409059/1‑150 (MQ=38) cGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGTTAAAAttt < 2:1020076/150‑1 (MQ=255) gCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtg < 2:987982/150‑1 (MQ=255) gCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtg > 2:583803/1‑150 (MQ=38) gCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtg > 1:504965/1‑150 (MQ=38) gaGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGGGATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt < 2:825982/149‑1 (MQ=255) cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt < 1:169413/150‑1 (MQ=39) ttACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAg > 2:68913/1‑150 (MQ=37) ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt < 2:173430/150‑1 (MQ=34) gTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc < 2:672457/150‑1 (MQ=21) gTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc < 2:1149436/150‑1 (MQ=25) tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt > 1:580339/1‑150 (MQ=35) aGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGCTATAGACTGAACACCAGTTTCACTGCtgg > 1:671260/1‑113 (MQ=38) aGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtgg < 2:671260/113‑1 (MQ=37) aGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg > 2:1285420/1‑150 (MQ=35) aGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTAGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCGGGTGGATCAGGCTAAGGTAAAGTTTGTGAGTTCTGCTCg > 2:519993/1‑150 (MQ=255) aaCTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaa > 2:1012664/1‑150 (MQ=35) aCTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa < 2:783824/150‑1 (MQ=35) ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga > 2:662876/1‑150 (MQ=21) ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTCCTCGAAaga > 2:232826/1‑150 (MQ=25) aTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACg > 1:458519/1‑150 (MQ=25) aTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga < 2:185419/144‑1 (MQ=35) aTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga > 1:185419/1‑144 (MQ=35) caaCATTCGAGCT‑GAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTa > 1:945830/2‑150 (MQ=2) aaCATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt < 1:1208259/124‑1 (MQ=37) aaCATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt > 2:1208259/1‑124 (MQ=37) cATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGa < 2:567775/150‑1 (MQ=21) aTTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa < 1:1173942/123‑1 (MQ=35) aTTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa > 2:1173942/1‑123 (MQ=35) aTTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGaa > 2:1178650/1‑150 (MQ=35) tC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtt > 1:541300/1‑150 (MQ=18) c‑AGCTCGAATGTTGATTTCTGGCTTTTGGCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGGAAAATTTGTGAGTTCTGCTCGAAGGAGCGACGTACGAAtct > 2:10205/1‑148 (MQ=11) cTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCg > 1:1157467/1‑150 (MQ=11) cTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCg > 2:758457/1‑150 (MQ=11) cTCGAATGTTGAGTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCg < 1:199496/150‑1 (MQ=11) | GTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCG > NC_002947/699580‑699808 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |