Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 982,560 | T→A | 23.6% | intergenic (‑122/+481) | ompF ← / ← asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 982,560 | 0 | T | A | 23.6% | 62.8 / 17.8 | 34 | intergenic (‑122/+481) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base T (13/13); new base A (3/5); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.93e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01 |
GCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTAC > CP009273/982416‑982704 | gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAt > 2:498724/1‑149 (MQ=255) ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAtt > 2:472358/1‑149 (MQ=255) ttGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTGCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCa > 1:493500/1‑143 (MQ=255) tggcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCCGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTaa < 2:109289/147‑1 (MQ=255) gcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTTTTTTTTTTGCCACTGCCGTCAATAATTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCcaa > 1:99532/1‑143 (MQ=255) ttCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTaa < 1:31522/143‑1 (MQ=255) ttCATCATGATTTATTACCCTCCTGGTTTTTTTGATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTc < 2:28886/149‑1 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAAGGTTTCg > 2:167601/1‑149 (MQ=255) ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCg > 2:239007/1‑149 (MQ=255) ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCAGCTTTCCATTCAAACTAACGTTTCGCTACCGTGaaa < 1:352260/143‑1 (MQ=255) aTGGTTTTTTTTAGGACACCTGCCACTGCTGTCCATAAGTTCTGTCAAGCAAAATTTCCGGAACTATTCATGCGAGTTTGGTGGCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACCTTTCGCTACCGTGAAAGTGc < 1:458516/143‑1 (MQ=255) tttttcATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaa < 1:366829/143‑1 (MQ=255) tttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaag < 1:447787/143‑1 (MQ=255) ccTGCCACTGCCGTCAATAAGTTCTGTCAATATAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACGAACGTTTCGCTACCGTGAAAGTGCTACAAAGATATAGATTTGGTTTCa < 2:433928/149‑1 (MQ=255) cTGCCACTGCCGTCAATAAGATCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTgg > 1:206557/1‑143 (MQ=255) aCTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTAGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTGCACAGATAAAGATTTGGGTTGAAAAAgg > 2:438093/1‑148 (MQ=255) tGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaa > 1:206818/1‑143 (MQ=255) ccGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaaa < 1:524557/143‑1 (MQ=255) cGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaaaa < 2:547978/149‑1 (MQ=255) aGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAata > 1:147476/1‑143 (MQ=255) aaaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAg > 1:411252/1‑143 (MQ=255) aaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAtt > 2:14959/1‑149 (MQ=255) aTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTaa < 2:395008/149‑1 (MQ=255) cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAAc > 2:575193/1‑149 (MQ=255) aaCTATTGATGAGAGTTTGGTGTCTTTAGGGGGCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTCCCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACttt < 2:67842/149‑1 (MQ=255) aaCTATTGAGGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACttt < 2:39891/149‑1 (MQ=255) tATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAAc > 1:453328/1‑143 (MQ=255) gTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac < 1:573993/143‑1 (MQ=255) gTTTGGTGTCTTTATGTGCCTGCAGGCATCTTACCATTCAAACTTTCGTGTCTCTACCGTGAAAGTGATACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac < 1:428148/143‑1 (MQ=255) tCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCAt < 2:459652/149‑1 (MQ=255) tATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTcc < 1:423527/143‑1 (MQ=255) gtgtCTGCAGGCATCTTTTCATTCAAACTAAAGTTTTGCTCCCGTGAAAGTGCTACAAAGATAAAGGTTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAAACTTGTGAACACCGTCATATTTcacta > 1:412620/1‑139 (MQ=255) gCATCTTTCCATTCAAACGCACGTTTCGCGACCGTGAAAGTGCGACAAAGATAAAGATTTGGTTTAAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt < 2:515772/149‑1 (MQ=255) tCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTAc > 2:170340/1‑149 (MQ=255) | GCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTAC > CP009273/982416‑982704 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |