Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 1,070,246 | 2 | . | A | 77.8% | 70.4 / 16.2 | 27 | intergenic (+46/+96) | aroP‑I/PP_0928 | aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger related‑protein |
Reads supporting (aligned to +/- strand): ref base . (0/6); new base A (14/7); total (14/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.80e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACTGGGCGAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGT > AE015451/1069952‑1070490 | aactggtcaaccagcAGCATTACCCACAGCAGATTCCGCAAGGGCAAGCTGGCCGCCGGGATCACCCCCTTCTACAGAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGGCGCCGCCGGTTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑Gagatc < 1:499031/287‑5 (MQ=255) aCTGGGCGAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑Gagatcg < 2:681889/301‑6 (MQ=255) tGGGCGAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑Gagatcgag < 1:66938/301‑8 (MQ=255) ataccgtcaGGCCAAGCTGGCCGCCGTGATCACCCCCTTCTACACAAGCCTGGGGCACTCGCGTACCAACTACCTGTGGCTGACGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg < 1:571284/293‑1 (MQ=255) aTCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCAGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCCCGGTCGGCGCCCGATCCCCCAGCCGATGAAAGCCTAAAGTCAtccgcctg > 1:348478/1‑294 (MQ=255) aaaaGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTACAGGCATGCGTCTGGCGCGCCTGACGCGt > 1:134319/1‑301 (MQ=255) ggcggCACCCACTGACCAACTACCGGTGCCTGGCGTTCATCGTGCTGGTGCTGGTAGGGATGTACATGACGACGCCGATTCGCATCGCGGTAATGCGGATCCCGGAATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt < 2:409506/298‑1 (MQ=255) ccAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTTAGATCGAGCGCCGCGCGGGAGGCGCCCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACTCGTTCCCCTTCCCTGccccctcatt > 2:146644/1‑295 (MQ=255) aaCTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCACATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct > 1:552090/1‑301 (MQ=255) ccTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCACATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa < 2:551904/301‑1 (MQ=255) tGGCGTTCAGCGTGCGTATCCTGGTAGGTATGTACCTGACGCCGCCGATTCGCATCTCGGTAAGTCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc < 2:137080/301‑1 (MQ=255) ggTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAac > 1:819814/1‑301 (MQ=255) gtgcgtccgATTCGAATGTCGGAAATGTTGTGCCCGGCATGGATGGATGTGCTGTGGGTGACTTTCAAGCTGAAGAAGGCTAGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCg < 1:752255/295‑1 (MQ=255) cgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCCCGCACAGTGGTTTACA‑GCAACACCGTACGGAGCCCCGGAAGTAc > 2:729097/1‑301 (MQ=255) cgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCCCCGCCACGCACGCACAGTGGTGT‑CACGCAACACCGTAGCCAGCAcccgccgcac > 2:630403/1‑293 (MQ=255) tCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCCCAATTCTGCCTCCGCCAAGCACGCACAGTGGGGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACCTAg > 1:790203/1‑301 (MQ=255) ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATc > 1:617124/1‑301 (MQ=255) aaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGTC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGGTCCCCTTCCATGCCCCTCAATGCTGCTCCCGCCAAGCACGCACCGTGGGGTACA‑GCAACACCGTAGCGAGCACCGGCAGTATACAGAACAGCATCaaa > 1:827483/1‑300 (MQ=255) ggATCCCGGCATGGATGGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTc < 2:819524/300‑1 (MQ=255) ctggatgggatggcgtgtgtgtgGGTTGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGGTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAggc < 2:616911/282‑1 (MQ=255) tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACgcc > 1:767831/1‑299 (MQ=255) gtgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCGGT‑‑GA‑‑GATCGAGCGCCGCGCGGGAGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATAccc < 1:146695/301‑1 (MQ=255) gtggggttCAAGCCGAAGAAGGCTCGCAAGGCAAAGGGCTGATCGGGTTTGTTCGGAAGGGGGATGCCGTTAGGTTTCCGGCGCC‑‑TGGGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATgg < 1:567668/295‑1 (MQ=255) gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATgg > 1:521295/1‑301 (MQ=255) gggttttCAAGCTGAAGGAGGCTCGCCAGGCCAAGGGGTGATCGGGTTTGTTCGTAGGGGGTATGTTCTAAGGGTTTCTGAGCC‑‑TGCGA‑‑GATCGCGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGACTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGt < 2:437381/297‑1 (MQ=255) ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGAACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGACACACAAAGGCCGTAGGcc > 1:171698/1‑301 (MQ=255) aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCCTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCATGCGGCCAGCCAACGGCATACCCATGGTAACAGCCGCCACATGAAGGCCGTAgtacagga > 1:671768/1‑294 (MQ=255) | AACTGGGCGAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA‑GCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGT > AE015451/1069952‑1070490 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |