Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,340,084 | C→T | 40.0% | intergenic (+69/‑262) | 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,084 | 0 | C | T | 40.0% | 21.5 / 27.4 | 30 | intergenic (+69/‑262) | ribAB‑II/PP_3814 | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base C (9/9); new base T (8/4); total (17/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01 |
GGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGA > NC_002947/4339949‑4340218 | 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 < 4:233364/122‑1 (MQ=21) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg > 3:233364/1‑122 (MQ=21) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 6:298840/1‑149 (MQ=25) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 4:300643/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 > 1:283869/1‑146 (MQ=25) cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa < 2:283869/146‑1 (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) ggTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataata < 4:174485/122‑1 (MQ=255) ggTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataata > 3:174485/1‑122 (MQ=255) accaccGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCACCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt < 6:6901/134‑1 (MQ=255) accaccGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCACCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt > 5:6901/1‑134 (MQ=255) | GGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGG‑GTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGA > NC_002947/4339949‑4340218 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |