Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,340,075 | A→T | 23.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 | 23.4% | 34.0 / 8.3 | 28 | 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 (12/8); new base T (4/2); total (17/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01 |
GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTT > NC_002947/4339929‑4340216 | gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGa < 5:170573/147‑1 (MQ=255) gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGa < 7:263811/147‑1 (MQ=255) gccgccGCTGACATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGGTCGACCCGGTCATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGt < 3:306110/149‑1 (MQ=255) ccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGCTGCTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATCCGAGGGAGaca > 1:337004/1‑149 (MQ=255) cTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGAcacca < 8:93747/149‑1 (MQ=255) gAAATAGGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTGGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGAcaccac < 2:337004/148‑1 (MQ=255) gAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTTACCCGCCAGGGCGGGGGTGaccacc > 1:156739/1‑145 (MQ=255) aTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATTGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGAc < 2:227451/149‑1 (MQ=255) ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc < 4:265414/149‑1 (MQ=255) gCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTTCCCGCGAATGCGAGGGTGACACCACCGCCGCCTTCGCGGGGAAACCCGCT‑CCTACAGGGCCCTTg > 8:294491/1‑149 (MQ=255) cATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTTACCGCGAATGCGAGGGTGTCACCACCGACGCATTAGCGGGTAAAACCCCG‑CCGAAAGGGCCCATGGTGcgt > 7:170997/1‑147 (MQ=38) aTCCCGTTAGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTAACCCCGAAAGCGAGGGTGAAAACACCGAAGCATTAGCGGGTAAAACCCCT‑CCTACACGGCCCATGGTGCttt > 1:35346/1‑149 (MQ=37) ccGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCCCGAATGCGAGGGTGACACCACCGACGCAATCGCGGGTAAACCCCCT‑CCTACACGGCCCTTGGTGCTTTTGACAGGGCg > 5:170980/1‑147 (MQ=255) ggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTTACCGCGAAATCGAGGGTGACACCACCCACGCATTAGCGGGTAAAACCGCTACGTAAAGGG‑CCGTGGTGCTTTTTAAAAGGGAGAGTGATTGc > 2:247151/1‑149 (MQ=255) gaatttGATTTCGCGGGGCCGTGTTGGCGGGGG‑‑TTTCCCCGCTATGGCGTGGGTGACCCCCCCGACCCATGCGGGGGAAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGGGTGATGGCCGGTCGCCCCATACCTCTTg < 8:170997/145‑1 (MQ=255) ccgctcccccTGTAGGCGCGGT‑‑TTACCCCGCTAATGCGCGGGTGACACCCCCGACGCATGCGCGGTAAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 5:2849/142‑1 (MQ=255) gCGGGGCCGTTTTGGTGGGGG‑‑TTTACCCGCTAATGCTAGGGTTACACCACCTACGAATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 2:35346/148‑1 (MQ=255) cTGTAGGAGCGGG‑‑TTTAACCGCGAATTCGAGGGTGACACCACCGACGCATTAGCGGGTTAACCCGCG‑ACTAAACGGGCCTTGGTGGTTTTTACAGGGCGGTGTGATGGCCGGTCGGCCCAAACCTCTTTCCAAAATCTTGGAATACCa > 2:60686/1‑148 (MQ=255) tGTAGGAGCGGG‑‑TTTAAACGCGAAAGCGAGGGTGGCACCACCCACGCATTAGCGGGTATACCCGCT‑CCTACACGGCCCGTGGTGCGTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCCGACCTCATGGCAAAATCTTGGAATAACaa > 6:289579/1‑147 (MQ=255) attcgCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCGAAAATCTTGGAATACCataata > 5:231067/5‑149 (MQ=255) ggAGCGGT‑‑TTTACCCGCGATTTCGAGGGTGACACAACCGACGCATTCGGGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGAAAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataata < 1:13363/148‑1 (MQ=39) g‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATcc > 2:24714/1‑149 (MQ=255) cccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATAACCCGTa < 7:294491/148‑1 (MQ=255) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 3:76886/1‑149 (MQ=25) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 3:73137/1‑149 (MQ=25) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAAACttc > 4:25888/1‑121 (MQ=18) cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa > 2:124726/1‑146 (MQ=25) cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa < 1:124726/146‑1 (MQ=25) cgAATGCTAGGGTGACACCAACGACGCATTCGCGGGTAAACCCGAT‑CATACAGGGCCGTTGGTGCTTATGAAAGGGCTGACTGATGGCCGCTCGCCCCATACCTCTTTACAAAATCATGGAAGACCATAATATGATATCCAGTCGAcat > 8:138576/1‑147 (MQ=255) gAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt > 6:166442/1‑149 (MQ=255) | GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCT‑CCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTT > NC_002947/4339929‑4340216 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |