| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,544,489 | 0 | C | A | 32.0% | 22.0 / 18.9 | 25 | intergenic (+48/+66) | mdfA/ampG | multidrug MFS transporter/muropeptide permease AmpG |
| Reads supporting (aligned to +/- strand): ref base C (9/8); new base A (6/2); total (15/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.02e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGC > NC_002947/1544353‑1544629 | tGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTTAACCCGCTCCTAAACAGGgccgcgccg > 5:99228/1‑149 (MQ=255) ggggggggggggggCGTGGGCGCGGGGGCCGCTGTGGTTTGGGGGGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa < 7:135619/149‑1 (MQ=38) ggggggggggggggTTGGCGGTGGGCCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAgg < 4:296762/149‑1 (MQ=21) ggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCCCATAGCCGGGTTAACCC‑CACCAACACAAGGCCCCGCCGCGTTGTTggg > 5:293579/1‑148 (MQ=255) ggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTa < 2:131221/149‑1 (MQ=255) gtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAACCCCCAACAACACACGGCCGCGCCGCTTTGTATGAGCGGGTTTTCACGCGAAGAgg > 7:189114/1‑149 (MQ=38) cGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAAAACCCAGCTAAACAGGGCCGCGCCGCGTTGTAGGAGCGGGGTTTCCCCCGCAGAGGCCAg > 3:63637/1‑148 (MQ=255) gttttGTGGTGGGTGGGCGGCCTGGGTATGGTGGTGTAGGTCGATAATGTGTATTGCCGGGGCGGGCGTTTTCG‑GGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGc < 6:195399/147‑1 (MQ=25) tgggtgggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTGCTAAACAGGGCCGCGCCGCGTTGTAGGAGCGGGTTTACACCCGAAGAGGCCAGCAAACGCAATAAAAAg > 7:181283/1‑149 (MQ=37) tgggGGGCCTGGGTATGGGGCTGTAGCTCGATAATGGGTATTGCCGGGGCCGGCCTCTTCG‑GGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTc < 6:99228/149‑1 (MQ=37) gttgggcGTGGGTAAGGGGTGGTTGTTCTATATTCGGTATGGCGGGGGCGGGCTTGTTCG‑GGGTTAAACCCGCTCCAACACCGGGCGGGGCGTTTTTTTTGTAGGGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCa < 6:293579/145‑1 (MQ=255) tAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTATACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg < 2:280455/149‑1 (MQ=255) ttCG‑CGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt > 7:155698/1‑149 (MQ=255) cg‑cgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc > 1:125162/1‑145 (MQ=255) aaCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg > 1:88113/1‑149 (MQ=255) aaCCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACCAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg > 1:200313/1‑149 (MQ=18) cGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg > 6:4229/1‑149 (MQ=21) cGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg > 7:141177/1‑149 (MQ=37) cGCTCATACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg < 4:63637/149‑1 (MQ=255) gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg < 5:68431/149‑1 (MQ=21) gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg > 6:68431/1‑149 (MQ=21) tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCGCTGCAGGGc > 5:301127/1‑148 (MQ=255) tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc > 5:277810/1‑148 (MQ=21) tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt > 1:146331/1‑132 (MQ=25) tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCg < 7:127716/79‑1 (MQ=255) | TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGC > NC_002947/1544353‑1544629 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |