Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,948,683 | T→A | 8.5% | intergenic (‑64/+18) | PP_mr62 ← / ← trxA | Pseudomon‑Rho regulatory region/thioredoxin I |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,948,683 | 0 | T | A | 8.5% | 137.7 / 3.4 | 47 | intergenic (‑64/+18) | PP_mr62/trxA | Pseudomon‑Rho regulatory region/thioredoxin I |
Reads supporting (aligned to +/- strand): ref base T (20/23); new base A (2/2); total (22/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01 |
GATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCGGG > NC_002947/5948541‑5948827 | gATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCggg < 5:235724/149‑1 (MQ=255) aTGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCgggg > 7:146001/1‑149 (MQ=255) aCAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCtttt < 4:245877/149‑1 (MQ=255) cAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGTCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCtttt > 1:58490/1‑148 (MQ=255) cAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTCTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCtttt < 1:277848/149‑1 (MQ=255) gATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTttaca > 2:344956/1‑149 (MQ=255) gATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCAAATGCGGGGCttttctttt < 6:223366/149‑6 (MQ=255) aTCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTttacat < 1:344956/149‑1 (MQ=255) cTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAGGCCCCGCAAATGCGGGGCttttct > 5:31728/1‑135 (MQ=255) cTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAGGCCCCGCAAATGCGGGGCttttct < 6:31728/137‑3 (MQ=255) aTTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTgg < 6:354214/147‑1 (MQ=255) ggATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATg < 8:146001/149‑1 (MQ=255) gaggcgtCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTCCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGc < 8:83648/143‑1 (MQ=255) aTGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCg < 3:229499/149‑1 (MQ=255) aTGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCg < 6:94865/149‑1 (MQ=255) tGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGa < 6:54510/149‑1 (MQ=255) aCTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGcc > 8:215442/1‑148 (MQ=255) gAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGTGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCt > 5:72484/1‑149 (MQ=255) aaaGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg > 7:319391/1‑149 (MQ=255) aGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt < 8:319391/149‑1 (MQ=255) gCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGt > 4:181543/1‑149 (MQ=255) cGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTg < 7:215442/149‑1 (MQ=255) cGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTg < 2:58490/149‑1 (MQ=255) cGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTg < 3:181543/149‑1 (MQ=255) gTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTc > 3:77404/1‑137 (MQ=255) gTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTc < 4:77404/137‑1 (MQ=255) cacTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg > 3:85268/1‑125 (MQ=255) cacTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg < 4:85268/125‑1 (MQ=255) ttttCTGCGTCTAGTGCTCAGAAAAAGAAAAGCCCCGCATTTGCCGGGCCTTTTCACATCACAGCTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCAAACCTTGGTGGCATCGACGTAGCCGTACTTGAACAGCAt > 4:91831/1‑148 (MQ=255) cTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGc < 6:292641/148‑1 (MQ=255) aGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATa > 8:114533/1‑149 (MQ=255) gTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATa > 3:272884/1‑148 (MQ=255) tGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCtt < 4:93151/125‑1 (MQ=255) tGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCtt > 3:93151/1‑125 (MQ=255) ctctGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCa < 7:332229/148‑1 (MQ=255) gAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGca > 5:6673/1‑148 (MQ=255) aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGGTTCGACGTTGCCGTTCTTGAACAGc > 8:226160/1‑125 (MQ=255) aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGc < 7:226160/125‑1 (MQ=255) aaaagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGaa < 2:87269/120‑1 (MQ=255) aaaagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGaa > 1:87269/1‑120 (MQ=255) aagaaaagACCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGAGTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACACAGCACAcct > 8:155960/1‑148 (MQ=255) aaaaggCCCCGCAAATGCGGGGCTTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTg > 3:208922/2‑149 (MQ=255) aaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGc > 7:113760/1‑149 (MQ=255) ccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCg > 8:95484/1‑149 (MQ=255) gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACa > 3:8535/1‑108 (MQ=255) gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACa < 4:8535/108‑1 (MQ=255) gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCggg < 8:113760/149‑1 (MQ=255) | GATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCT‑TTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCGGG > NC_002947/5948541‑5948827 |
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: ‑ |