Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,340,075 | A→T | 41.4% | intergenic (+60/‑271) | ribAB‑II → / → PP_3814 | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,340,075 | 0 | A | T | 41.4% | 16.0 / 30.3 | 29 | intergenic (+60/‑271) | ribAB‑II/PP_3814 | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base A (8/9); new base T (8/4); total (16/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.51e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01 |
GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGA > NC_002947/4339929‑4340218 | gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGa < 4:56768/147‑1 (MQ=255) gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGa < 4:174498/147‑1 (MQ=255) gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAg < 4:187327/148‑1 (MQ=255) ggggCGGGCGGGCTATGAGCTGGAGGGGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTc < 5:418938/149‑1 (MQ=255) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGTCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg < 1:272516/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg < 8:342524/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg < 5:407129/149‑1 (MQ=255) gAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc < 3:25742/129‑1 (MQ=255) gAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGAGGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc > 4:25742/1‑129 (MQ=255) acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg > 1:147723/1‑80 (MQ=255) acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg < 2:147723/80‑1 (MQ=255) gAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCGTTTGACAGGGCGGAGTGATGGCCGGTGGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAatat > 7:293949/1‑148 (MQ=255) agagggTTTCCCGGCGAATGCGAGGCGTACCCCACCGCCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGa < 8:293949/146‑1 (MQ=255) gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAg > 8:362539/1‑90 (MQ=255) gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAg < 7:362539/90‑1 (MQ=255) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg > 3:233364/1‑122 (MQ=21) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 4:233364/122‑1 (MQ=21) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 4:300643/1‑149 (MQ=25) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 6:298840/1‑149 (MQ=25) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACAATAATATGATATCGCGTAGa > 8:228793/1‑149 (MQ=25) ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 4:45374/1‑149 (MQ=255) ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGAGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 3:357417/1‑149 (MQ=255) ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGAGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa < 4:357417/149‑1 (MQ=255) cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa < 2:283869/146‑1 (MQ=25) cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa > 1:283869/1‑146 (MQ=25) cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc > 3:219353/1‑149 (MQ=25) gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc > 3:224520/1‑149 (MQ=25) gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCCTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTTCCAAAATCTTGGAATACAATAATATGATATCCCGTAGAca > 6:349425/1‑148 (MQ=21) gcgAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc > 4:392815/1‑149 (MQ=255) gAATGCGAGGGTGACACCACCGACGCATTCGCGGTGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt > 1:54237/1‑149 (MQ=255) aTGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa < 3:45374/149‑1 (MQ=255) | GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGA > NC_002947/4339929‑4340218 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |