Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,820,642 | Δ1 bp | 12.9% | intergenic (+400/+112) | PP_2473 → / ← PP_2474 | hypothetical protein/glutathione S‑transferase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,820,642 | 0 | A | . | 12.9% | 84.8 / 8.5 | 31 | intergenic (+400/+112) | PP_2473/PP_2474 | hypothetical protein/glutathione S‑transferase family protein |
Reads supporting (aligned to +/- strand): ref base A (14/12); new base . (2/2); total (16/15) | |||||||||||
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.91e-01 |
ATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTT > NC_002947/2820504‑2820783 | aTGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAGTCAAAAACCCGGGTTCG‑ACCCGGGtttt < 1:222728/150‑1 (MQ=255) cccAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTTCG‑ACCCGGGtttttgc < 2:524164/150‑4 (MQ=255) ttCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCg > 2:462491/1‑137 (MQ=255) ttCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCg < 1:462491/137‑1 (MQ=255) atatatCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTg > 2:650909/1‑150 (MQ=255) tatatCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGa < 1:566903/150‑1 (MQ=255) tAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGCGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGt < 1:650909/150‑1 (MQ=255) gggCTATTTGCTGGGTGATATAGAGAACGTGAGTCTGTAAAAAAGTAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGa < 2:130818/150‑1 (MQ=255) gggCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGa > 1:256889/1‑150 (MQ=255) gCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGctct < 1:571092/125‑1 (MQ=255) gCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGctct > 2:571092/1‑125 (MQ=255) cTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGa > 1:230015/1‑111 (MQ=255) cTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGa < 2:230015/111‑1 (MQ=255) cTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGc > 2:673537/1‑150 (MQ=255) tATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCAAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTgg > 1:272391/1‑98 (MQ=255) tATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCAAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTgg < 2:272391/98‑1 (MQ=255) ggctgATATAGAGAAAGTGCGCCTCAAAAAAAGCAAAAACCCGGG‑TCGAACTCGGGTATTCGGTTGAGTCCCTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCATTGCTGCTCTCAGACCGGCGTTGGTTTCTTGGACGCTCTTTTTGc < 2:256889/150‑1 (MQ=255) tgATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATgctgct > 1:161492/1‑110 (MQ=255) tgATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATgctgct < 2:161492/110‑1 (MQ=255) tgATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGaa < 1:673537/150‑1 (MQ=255) gagaAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGt > 2:422280/1‑150 (MQ=255) agaAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGtt > 2:676596/1‑150 (MQ=255) aaaGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTa > 2:600208/1‑150 (MQ=255) tGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTg < 1:422280/150‑1 (MQ=255) cGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGc > 1:794052/1‑150 (MQ=255) tGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCt > 2:462890/1‑150 (MQ=255) tcaaccgaaaaCCCGGGTTCG‑ACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATc > 1:240755/8‑150 (MQ=255) aGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAg < 1:462890/150‑1 (MQ=255) aCCCGGGTTCG‑ACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCAt > 1:552507/1‑133 (MQ=255) aCCCGGG‑TTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCAt < 2:552507/133‑1 (MQ=255) ccGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTTTTTATTTTGTTTTAATTGACATCTATATCCTATGCTCGCtt > 2:327219/1‑150 (MQ=255) | ATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTT > NC_002947/2820504‑2820783 |
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: ‑ |