| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_pae_tpiA | 4,110,974 | 0 | C | A | 57.1% | ‑2.9 / 77.2 | 50 | intergenic (+289/+557) | cdh/yiiQ | CDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein |
| Reads supporting (aligned to +/- strand): ref base C (11/10); new base A (13/15); total (25/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.76e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGCTTTAGCCCTGGCAGTTAGTTTAGCTATT > NC_000913_3_pae_tpiA/4110832‑4111118 | cccTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCgg > 1:868623/1‑150 (MQ=255) ccTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGTGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGt < 2:922973/150‑1 (MQ=255) ccTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGt < 1:329472/150‑1 (MQ=255) ccTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGt > 4:482601/1‑150 (MQ=255) cTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTc > 2:270834/1‑150 (MQ=255) tCGTGCCGAGTTTACGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCTCTAGCTGGACAGTTTCCGCCACGGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTcc < 2:390590/150‑1 (MQ=255) tCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTcc < 3:482601/150‑1 (MQ=255) tGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGa < 1:270834/150‑1 (MQ=255) ccGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATc > 2:663656/1‑150 (MQ=255) gAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAg > 1:106973/1‑150 (MQ=255) ttCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAGCAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAAc > 1:850177/1‑150 (MQ=255) cGTCCTCTCCTCTCTCAATACCCTCCCCGGTTATCACGTCGACTTCATTGAATTTCCGCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGa < 1:348991/150‑1 (MQ=255) ccTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTg > 1:369300/1‑150 (MQ=255) ctctCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTTGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTcc < 2:625540/113‑1 (MQ=255) ctctCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTcc > 1:625540/1‑113 (MQ=255) tctcCGCTCTCAATACCCGCCCCAGTNATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCa > 1:613697/1‑150 (MQ=255) cGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCGCTCCCGGTCCCGATCAGCTCggg > 3:388753/1‑127 (MQ=255) cGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACGTGGAGACAGCACTCCCGGTCCCGATCAGCTCggg < 4:388753/127‑1 (MQ=255) cccGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGggcg > 2:927258/1‑150 (MQ=255) ttattaCGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCaa > 1:139929/1‑150 (MQ=255) ttattaCGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCaa < 2:868623/150‑1 (MQ=255) aGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGggcggcggcggc < 2:45579/132‑1 (MQ=255) aGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGggcggcggcggc > 1:45579/1‑132 (MQ=255) aGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATc > 1:832830/1‑150 (MQ=255) ggAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCgg > 2:723714/1‑150 (MQ=255) aaGTGCCTCCCCACTAGCTGGACAGTTTCCACAGCAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCgggg < 2:850177/150‑1 (MQ=255) aaGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCgggg < 2:369300/150‑1 (MQ=255) tGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGgcg < 1:927258/150‑1 (MQ=255) tGCCCCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCGTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGgcg < 4:573949/150‑1 (MQ=255) tCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAg < 1:612579/73‑1 (MQ=255) tCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAg > 2:612579/1‑73 (MQ=255) ccACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATaa < 4:568126/150‑1 (MQ=255) cACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAAc < 2:832830/150‑1 (MQ=255) cACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAAc > 2:729391/1‑150 (MQ=255) aCTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACt < 2:106973/150‑1 (MQ=255) tttCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAg > 2:6162/1‑150 (MQ=255) caacaGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGCtt < 1:788353/150‑1 (MQ=255) aacaGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGc > 1:55870/1‑66 (MQ=255) aacaGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGc < 2:55870/66‑1 (MQ=255) aacaGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCa < 2:831710/123‑1 (MQ=255) aacaGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCa > 1:831710/1‑123 (MQ=255) caGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGc < 2:875424/145‑1 (MQ=255) caGTCTCGGGTTTATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGc > 1:875424/1‑145 (MQ=255) cGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGCTTTAGCCCTgg > 4:543855/1‑150 (MQ=255) tttATACTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGCTTTAGCCCTGGCagt < 2:613697/150‑1 (MQ=255) aCTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGCTTTAGCCCTGGCagttagtt < 1:911454/150‑1 (MQ=255) aCTTGGAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGCTTTAGCCCTGGCagttagtt > 4:476308/1‑150 (MQ=255) ggAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTg < 2:519849/82‑1 (MQ=255) ggAGACAGCACTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTg > 1:519849/1‑82 (MQ=255) aCAGCCCGCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGCTTTAGCCCTGGCAGCTAGTTTAGCTAtt > 4:473547/1‑150 (MQ=255) | CCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGTTATTACGGCGACGTCAGGGAAGTGCCTCCCCACTAGCTGGACAGTTTCCACAACAGTCTCGGGTTTATACTTGGAGACAGCCCTCCCGGTCCCGATCAGCTCGGGGGGCTCAACGGCGACTGCATGCGGGCTTAGGGCGGCGGCGGCAAGGCTGGTCCTCGGATCGGGGGCGCAGATAACTACGTCAAGCCCAGCGGCTTTAGCCCTGGCAGTTAGTTTAGCTATT > NC_000913_3_pae_tpiA/4110832‑4111118 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |