Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,340,076 | G→C | 43.0% | intergenic (+61/‑270) | 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,076 | 0 | G | C | 43.0% | 12.3 / 34.7 | 28 | intergenic (+61/‑270) | ribAB‑II/PP_3814 | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base G (8/8); new base C (8/4); total (16/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.59e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01 |
GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGA > NC_002947/4339929‑4340218 | 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 < 2:147723/80‑1 (MQ=255) acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg > 1:147723/1‑80 (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 < 4:357417/149‑1 (MQ=255) ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGAGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 3:357417/1‑149 (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: ‑ |