Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,544,489 | 0 | C | A | 25.0% | 45.4 / 13.4 | 28 | intergenic (+48/+66) | mdfA/ampG | multidrug MFS transporter/muropeptide permease AmpG |
Reads supporting (aligned to +/- strand): ref base C (9/12); new base A (3/4); total (12/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG > NC_002947/1544359‑1544635 | gggggggggCGGGGCCGGTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt < 5:162158/148‑1 (MQ=25) gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTg > 6:199818/1‑147 (MQ=255) gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTAGCGCGTAAAACCCCTCCTAAACCCGGCCCCCCCGCTTTGTa > 6:295260/1‑149 (MQ=255) gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg < 6:82896/149‑1 (MQ=37) ggggggCGTGGCCGCGGGGGCCGCGGTGGTTTGTGTGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg < 8:81367/148‑1 (MQ=37) gcgcTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTTAAACCCCTCCGACCCACGGGCGCGCCGCGTTGTAGGAGCGGGGTTAACcccg > 7:27635/1‑148 (MQ=255) ggtgCCCGCGGGGGTTGGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAagag < 3:117624/148‑1 (MQ=25) tgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca < 7:345610/106‑1 (MQ=255) tgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca > 8:345610/1‑106 (MQ=255) gCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCATCGCGGGTAAACCCGCGCCAACACACGGCCCCGCCGCATTGTAAGAGCGGGGTTAACCGCGAAGAGGCCa > 2:292734/1‑149 (MQ=255) ggggTTTGTGGGGGGGGGGCGGCCTGGGTACGTGGCTGTCGCGCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAg < 3:257438/147‑1 (MQ=25) gggCGGCCTGGGTATGGGGGTGTTGCTCGATATTGGTTTTTGCGGGGGCCGGCCTCTTGGCGGTTAAACCCGCTCCTACACTGGGCCGCGCGGGTTTTTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCa < 8:27635/149‑1 (MQ=25) aagcgccaTAATCGGTATGCCCGGGGCCGGCCTCTTCGCGGGAAACCCCGCTCCAACCCCGGCCCGCGCCGCTTTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg < 5:37790/142‑1 (MQ=37) gATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGt < 7:208482/141‑1 (MQ=255) gCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTc < 5:295260/148‑1 (MQ=255) ctgcgcgGGAAACCCCGCTCCAACACCGGCCCGCGCCGCTTTTTCGGCGCGGTTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTcc < 1:161789/146‑1 (MQ=255) ttCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCt > 3:39507/1‑149 (MQ=255) ttCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCt > 7:308974/1‑149 (MQ=255) cgcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCt < 5:199818/147‑1 (MQ=255) gcgGGTAAACCCGCTCCTACCCAGGGCCGCGCCGCTTTGTAGGAGCTGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGAt < 8:308974/149‑1 (MQ=255) gggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGACACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGc < 1:292734/149‑1 (MQ=255) gTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCa > 3:107198/1‑149 (MQ=255) cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg > 3:107884/1‑147 (MQ=255) cgcgCCTACAACGGGGGGCGGCCCTGTGTTGTAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACCAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg < 5:181749/149‑1 (MQ=21) cGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg > 2:282266/1‑149 (MQ=255) cGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg > 2:108945/1‑149 (MQ=37) cGCTCATACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGGTTACCCGCGAAGAGGCCAGCACAGGCAATCCAAAGCTTCAGGGCCGCTCAGCCCCCGAGGGGTCGT‑CCGCCCCCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg > 5:274238/1‑149 (MQ=25) gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg > 6:126185/1‑147 (MQ=25) gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg < 5:126185/147‑1 (MQ=25) cTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg < 1:282266/149‑1 (MQ=255) acaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag > 6:264436/1‑148 (MQ=255) acaGGCCCGCGCCGCTTTGTAGGAGCGGGTTTACACGCGAAGAGGCCAGCACAGGCAATACAAATCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCAACTGCACGGCgataag > 2:128099/1‑148 (MQ=39) | CGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGT‑CCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG > NC_002947/1544359‑1544635 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |