Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 181,805 | T→C | 13.2% | intergenic (+19/‑87) | PP_r01 → / → PP_5SB | 5S ribosomal RNA/5S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 181,805 | 0 | T | C | 13.2% | 102.6 / 6.5 | 44 | intergenic (+19/‑87) | PP_r01/PP_5SB | 5S ribosomal RNA/5S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base T (21/17); new base C (2/4); total (23/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.03e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01 |
AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT > NC_002947/181664‑181954 | aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTACCTGcatgc < 2:25761/150‑4 (MQ=11) cgacCATAGAGCATTGGCACCACCTGTTACTATGCCGAACTCAGGAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATACTCAAGATTC‑ATTACGTAAAAACCCCATCTCCACATCCAGATaagc < 2:520749/150‑4 (MQ=11) cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTACCTGCATGTGCAGATAggg < 2:5536/150‑1 (MQ=11) ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtt > 2:725320/1‑150 (MQ=32) cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGttt > 1:830621/1‑150 (MQ=32) gagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTc < 2:830621/150‑1 (MQ=32) ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa < 2:610807/150‑1 (MQ=32) ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATCCAGGTAGGGGTTTTGTCTTTaa > 2:101492/1‑150 (MQ=17) ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa < 2:800106/144‑1 (MQ=32) ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa > 1:800106/1‑144 (MQ=32) accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAAc > 2:328605/1‑150 (MQ=35) cTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTa > 2:602490/1‑150 (MQ=38) ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga < 2:556417/148‑1 (MQ=40) ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga > 1:556417/1‑148 (MQ=40) catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGAt < 1:220597/150‑1 (MQ=255) tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTcc < 1:725320/150‑1 (MQ=255) gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTa < 1:328605/150‑1 (MQ=255) tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAg > 1:686774/1‑150 (MQ=255) tGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGc < 2:310549/150‑1 (MQ=255) aTCGCCGATGGTAGTGTGGGGTTTCCGCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACGCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAg > 2:530835/1‑150 (MQ=255) cGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTTGGTCATCGTCAAGATTC‑AGTACGGAGTACCCCTATGTGCGCATGTAGGGAGGGGTGTTGTATTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGca > 2:470778/1‑150 (MQ=18) gctgATGGTAGAGTCGGGGTTCCCCATGTGAGAGGAGGTCATCGTCAAGATTC‑ATTTCCAAAAACCCCTATCTGCACATGCAGGGAGGGGTTTTGTCTTTAAGTAGTAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGGCAggaa > 2:613011/4‑147 (MQ=18) gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCCAGCTTCTATTTC‑AAAAACCCCTACCGGGAGATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGGAGGCACGTCCTAGGGACGGGCCAGcacaca < 1:101492/150‑1 (MQ=21) gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt < 1:602490/150‑1 (MQ=255) tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAg > 1:802764/1‑74 (MQ=35) tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAg < 2:802764/74‑1 (MQ=35) cccATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACGTTCAGGTAGGGTTTTTGTCTTTAAGTATGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGACAGCGCACAGAATTTCTTGACGACCATa > 2:390132/1‑150 (MQ=25) ccATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCt < 2:499199/138‑1 (MQ=255) ccATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCt > 1:499199/1‑138 (MQ=255) ccATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATag > 2:677925/1‑150 (MQ=255) gagaGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCa < 1:640328/95‑1 (MQ=255) gagaGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCa > 2:640328/1‑95 (MQ=255) gagaCTAGGACATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTGGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCAAACAGCATTTCTTCACGACCATAGAGCAtt < 2:504435/150‑1 (MQ=25) gTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGaa < 1:677925/150‑1 (MQ=255) gTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccac > 1:58657/1‑150 (MQ=255) tCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGt > 1:57347/1‑60 (MQ=38) tCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGt < 2:57347/60‑1 (MQ=38) tCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc > 2:320582/1‑150 (MQ=255) ttC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatccc > 1:346044/1‑150 (MQ=255) aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGTGCATTGGAaccac > 2:637300/1‑127 (MQ=255) aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGTGCATTGGAaccac < 1:637300/127‑1 (MQ=255) ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGa > 2:841976/1‑150 (MQ=37) cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa > 1:97849/1‑138 (MQ=37) cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa < 2:97849/138‑1 (MQ=37) tCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt < 1:320582/150‑1 (MQ=255) | AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTC‑ATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT > NC_002947/181664‑181954 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |