Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 2,372,863 | C→G | 36.7% | G305R (GGA→CGA) | USA300HOU_2234 ← | NCS2 family nucleobase:cation symporter‑2 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 2,372,863 | 0 | C | G | 36.7% | 23.7 / 45.3 | 41 | G305R (GGA→CGA) | USA300HOU_2234 | NCS2 family nucleobase:cation symporter‑2 |
Reads supporting (aligned to +/- strand): ref base C (14/12); new base G (10/5); total (24/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01 |
GATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTATCAAAGAAATCAATAAATAAGAAC > CP000730/2372726‑2372996 | gATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTa < 2:337713/140‑1 (MQ=255) tAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTCACGATACTAGCAAATCCGGTTCTAGCACCTACAGCAACACCTGACGTTGATTCGATATAAGAAGTTGTTGTAGTAGTTCCAAAAACTGATACTCCac > 2:26013/1‑139 (MQ=255) aaCGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAAt > 2:107934/1‑140 (MQ=255) aCGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCGTACAATc < 1:179625/140‑1 (MQ=255) cTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTTATCCTACAATCGTAGc > 2:40999/1‑140 (MQ=255) aaGAACACCGCTAATAAGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTaa < 1:227015/140‑1 (MQ=255) aCAGGGCTAAGAGGAAACAAGAACCTGTTACGATACTAGCAAATCCTGTTCTCGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGaa < 2:104880/140‑1 (MQ=255) gcgcTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCt > 2:332020/1‑140 (MQ=255) cgcAAATAGCAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTACCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTg < 2:96606/140‑1 (MQ=255) aaTAGGAAACAACAACATGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGaa > 2:77412/1‑140 (MQ=255) acaacaacCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGctc > 1:296170/1‑140 (MQ=255) tcaacCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTGGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTAc < 2:152878/139‑1 (MQ=255) aacCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTCGCTAAAGAATCTGAGAACAATGCTCTACCa > 1:331450/1‑140 (MQ=255) aacCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCa > 2:74670/1‑140 (MQ=255) acCTGTTACGATACTAGCAGATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAg > 1:213402/1‑140 (MQ=255) ccTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGc > 2:114188/1‑140 (MQ=255) cTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTc > 1:108608/1‑140 (MQ=255) tGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTACTATCTTTCa < 1:261938/140‑1 (MQ=255) gTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCat > 1:216620/1‑140 (MQ=255) gTTCTAGCACCTACAGCAACACCTGAAGTTGATGCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCat > 1:71683/1‑140 (MQ=255) tttccccTACAGCAACTCCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAATAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCataata < 1:80296/135‑1 (MQ=255) cACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGGTGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTCCCAGCTCGTGCTAATGTGTTCTCTTTCATAATACCt > 1:342411/1‑140 (MQ=255) aCCTACAGCAACACCTGAAGTTGAGTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCGCTTGGTAATTTATTATCTTTCATAATACCTg < 2:164834/140‑1 (MQ=255) aCAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTg > 1:150285/1‑140 (MQ=255) aCAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAATTAATCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTg > 1:187771/1‑140 (MQ=255) gCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGt > 1:141630/1‑140 (MQ=255) ccTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTCGGTAATTTATTATCTTTAATAATACCTGCttacgtacccac > 2:22705/1‑130 (MQ=255) aGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGcta < 1:107934/140‑1 (MQ=255) tCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTCCTCATGt > 2:88020/1‑140 (MQ=255) tataAGAAGTTGTTGTAGTTGTTCCAAAAACTGATTGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAATAGCTACTAATGTACCa < 1:197628/140‑1 (MQ=255) tataAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCa > 2:90025/1‑140 (MQ=255) taAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGc < 1:77412/140‑1 (MQ=255) gttgttGTAGTTGTTCCAAAATCTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTATc < 1:327290/140‑1 (MQ=255) gttgttGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTATc > 2:9673/1‑140 (MQ=255) ctttttgtAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCGCTACCAGCTCTTGGGAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTATc < 2:149092/136‑1 (MQ=255) gtAGTTGTTCCAAAAACTGATCGTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCGCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCACCAGGTACTAATCTACCAGCGGTATCAAAGaa > 2:80271/1‑140 (MQ=255) cAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTATCAAAGAAATCaataaat < 1:181244/140‑1 (MQ=255) aaaaCTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTATCAAAGAAATCAataaataa < 1:74670/140‑1 (MQ=255) aaCTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTATCAAAGAAATCAATAAATAAGa < 2:296170/140‑1 (MQ=255) cTGATCGTACAATCGCAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTGTTAAAGAGATCAATAAATAAGAAc > 1:264134/1‑140 (MQ=255) cTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTCCCAGCTCTAGGTAATTAATTATCTTTCATAATACCAGCTTGAGTACCATCAGCTACTAATGTACCAGCTGTATCAAAGAAATCAATAAATACGaaa > 2:87503/1‑139 (MQ=255) | GATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAGAACAATGCTCTACCAGCTCTTGGTAATTTATTATCTTTCATAATACCTGCTTGAGTAGCAACAGCTACTAATGTACCAGCTGTATCAAAGAAATCAATAAATAAGAAC > CP000730/2372726‑2372996 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |