Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 982,566 | C→A | 23.2% | intergenic (‑128/+475) | 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,566 | 0 | C | A | 23.2% | 57.6 / 11.4 | 35 | intergenic (‑128/+475) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base C (14/12); new base A (3/5); total (17/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.88e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTT > CP009273/982420‑982712 | gAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCaa < 1:652539/149‑1 (MQ=255) atcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCt > 1:691314/1‑149 (MQ=255) tcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTa > 2:363981/1‑149 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAGCGTTTCGCTACCGTGaaa < 1:560532/149‑1 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGaaa < 2:292169/149‑1 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGc < 2:552425/149‑1 (MQ=255) tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACa > 2:212816/1‑149 (MQ=255) gTTTTGTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGCTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagat < 2:558627/149‑1 (MQ=255) ttttttATGACTCCTGCCACTGCCGTCAATACGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTATGCTGGCATCTTGCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa < 2:386657/149‑1 (MQ=255) ttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa < 1:84304/149‑1 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAtt > 1:542817/1‑149 (MQ=255) tGACACCTGCCACTGCCGTCAAGAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTgg < 2:781639/149‑1 (MQ=255) acCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc > 2:94923/1‑149 (MQ=255) cTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGGTAGGGGTATTTATGTGTCTGCAGGCGTCTTTCCTTTCAAACTAACGTTTCGCTACCGTGAAAGGGCTACCAAGATACAGGATTAGTTTCaa > 2:141433/1‑149 (MQ=255) tgcGTCAATAAGTTCTGTTGATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaa < 1:301821/147‑1 (MQ=255) ggcGTCAATAAGTTTTGTTAATAATAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaa < 2:295751/147‑1 (MQ=255) ccGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTCCAAAGATAAAGATTTGGTTTAAAAAAGaaaaa > 1:59253/1‑149 (MQ=255) gTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAat < 2:658616/149‑1 (MQ=255) tCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCTAAAAGAAAAAatt > 2:430161/1‑148 (MQ=255) cAATAAGTTCTGACAATAAAAAATCACGAAACTATTAATGAGAGTTTGAAGTCTTTATATGTCTGCAGGCAACTTTCCAATTAAACAAACCTTACGCTACCGTGAAAGTTCTACAAAGATAAAGATTTGGAAAAAAAAAAAAAAAatat < 1:281014/149‑1 (MQ=255) aaGTTCAGTCAAAAAAAATTGACGCAACTATTGATGAGATTTTTGTTGCTTTATGTTTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAACTGCTACAAAGATAAAAACTTGCTTTCAAAAAGAATATATATGTaa < 1:236498/149‑1 (MQ=255) cTGTTAATAAAAATTTACTGAACTATTGATGGGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCaaa < 1:283490/149‑1 (MQ=255) ttGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa > 2:806593/1‑149 (MQ=255) ttGATGAGAGTTTGGGGTTTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa < 1:123164/149‑1 (MQ=255) ttagagGTTGGTGTCTGTATTTGGCTGCAGGCATTTTTCCATTAAAACTAACGTTTTGCTACCGTGAAAGGGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTATAATTTAAGGAACTTTGTGATCAcc < 1:166592/147‑1 (MQ=255) tttGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGCACACCGTCata > 1:286112/1‑149 (MQ=255) tttGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCata > 1:310353/1‑149 (MQ=255) ggTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAttt < 1:212816/149‑1 (MQ=255) ggCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGatat < 1:363981/149‑1 (MQ=255) gCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt > 2:631388/1‑149 (MQ=255) gCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt > 1:786832/1‑149 (MQ=255) gCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTAAAAAGCTGAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt > 2:817206/1‑149 (MQ=255) ttCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCGAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGAAGAGATGATATTTACAGCAAtta > 2:840594/1‑148 (MQ=255) ttCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAAttt > 2:773090/1‑149 (MQ=255) ttCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAAttt > 1:52951/1‑149 (MQ=255) | GAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTT > CP009273/982420‑982712 |
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: ‑ |