| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 367,644 | G→T | 100% | intergenic (‑96/‑20) | nanA ← / → USA300HOU_0336 | N‑acetylneuraminate lyase/possible transcriptional regulator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 367,644 | 0 | G | T | 95.7% | 82.8 / ‑5.5 | 29 | intergenic (‑96/‑20) | nanA/USA300HOU_0336 | N‑acetylneuraminate lyase/possible transcriptional regulator |
| Reads supporting (aligned to +/- strand): ref base G (0/2); new base T (12/15); total (12/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.98e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGTT > CP000730/367512‑367779 | gTAACGCTGCATATCAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAGACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATAttt > 2:191969/1‑141 (MQ=255) actcttcATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATg < 2:237313/135‑1 (MQ=255) cGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGt < 2:14279/141‑1 (MQ=255) cTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGaaaa < 1:231295/141‑1 (MQ=255) aaaTCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATAttt > 2:133437/1‑141 (MQ=255) aTCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTAATATTTATGTATAATTAATGAAAATGttttt > 1:105379/1‑141 (MQ=255) tCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTc < 2:245958/141‑1 (MQ=255) aTTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAAt < 2:10826/141‑1 (MQ=255) aTTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAAt < 1:45737/141‑1 (MQ=255) tttaaGTCACTCCAAATATCATTTGATAATATAGTTTATAGCTATAATTGTAAACGCTATTATAAAATGTAACAATATCTTTTTGACAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGa < 2:35866/138‑1 (MQ=255) cccTCCAAAAATCATTTGATAATATTGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgag < 1:252992/141‑1 (MQ=255) aaataTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTTATATTTATGTATAATTGATGAAAATGTTTTTCAAAATCAATAGAAATTGAGTGAGTa < 2:54829/141‑1 (MQ=255) aaaaaTCATTTGATAATATAGTTTACAACTATAATTGTAAACGCTATCATAAAATGTAACAATATTTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTa < 1:159173/141‑1 (MQ=255) aaaTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAg > 1:67899/1‑141 (MQ=255) aTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTTTATATTTTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGGGGGAGTAAggg > 2:240479/1‑140 (MQ=255) ttGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTAtt < 1:194870/141‑1 (MQ=255) taATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAgg > 1:251259/1‑141 (MQ=255) tAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCa > 1:29021/1‑141 (MQ=255) tAAAAGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAAATCATGAAAATGTTTTTTAAAATCAATAGAAATTGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCa < 2:42528/141‑1 (MQ=255) aTCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTaaca > 1:138349/1‑139 (MQ=255) aTAAAATGTAACAATATCTTTTTGAAAATTGTATTCATGTTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCgg > 2:170213/1‑141 (MQ=255) aaaaaTGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGc < 2:129248/140‑1 (MQ=255) ataaCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACCTCGCAATCGATCTTGCAGGCCATCACCTTAAATCGGCAGTTATTGATaa < 2:163977/141‑1 (MQ=255) tatCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAg > 2:59545/1‑141 (MQ=255) tCTTTTTGAAAATTGTATTTATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCa < 1:234160/141‑1 (MQ=255) tttttGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAAGCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGAAGTAATTTATACGCact < 2:129767/141‑3 (MQ=255) gAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAat < 1:129408/141‑1 (MQ=255) aTTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGt > 2:237595/1‑141 (MQ=255) ttGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGtt > 1:208778/1‑141 (MQ=255) | GTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGTT > CP000730/367512‑367779 |
| 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: ‑ |