Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,577,123 | 0 | A | C | 30.0% | 83.8 / 15.9 | 50 | V166G (GTA→GGA) | yddB | putative TonB‑dependent receptor |
Reads supporting (aligned to +/- strand): ref base A (28/7); new base C (0/15); total (28/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.58e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCG > NC_000913/1576986‑1577251 | ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcc > 2:609567/1‑140 (MQ=255) ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcc > 2:193173/1‑140 (MQ=255) ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcc > 1:410839/1‑140 (MQ=255) ccTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcccc > 2:154289/1‑140 (MQ=255) aaGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATACCCC‑AATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCa < 1:193173/140‑1 (MQ=255) agagCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCAt > 1:21585/1‑140 (MQ=255) agagCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCAt < 2:410839/140‑1 (MQ=255) tccTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGa > 1:440969/2‑140 (MQ=255) gCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGa > 2:587442/1‑140 (MQ=255) aTGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt > 1:358619/1‑140 (MQ=255) ctctTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc > 2:396623/1‑140 (MQ=255) ctctTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc > 1:762162/1‑140 (MQ=255) ttATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGcac < 1:362622/140‑1 (MQ=255) aTTCTCATCGTTATGCGATGTTACCCGGTCGGAACGAGTAGTGCAATAACCA‑AATTTCCCCTTGCTATCACCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaa < 2:358619/140‑1 (MQ=255) tctcATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa < 1:561295/140‑1 (MQ=255) tctcATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa < 1:680463/140‑1 (MQ=255) cATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTg > 1:242754/1‑140 (MQ=255) atatGCGATGTTAACCAGCCCGAACGCGTAGTGCGATAACCC‑AATTTCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATa < 2:580501/140‑1 (MQ=255) cGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCg < 2:191765/140‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTca < 2:332678/140‑1 (MQ=255) ttAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca < 1:609567/140‑1 (MQ=255) ttAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTGGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca < 2:620566/140‑1 (MQ=255) ttAACCAGTCCGAACGCGTAGTGCGAAAACCC‑AATTTCACTTTGTTATCATCAGCGTTAAGCCGTTTGATCTTTGCAACAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca < 2:440969/140‑1 (MQ=255) tCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg > 2:563052/1‑140 (MQ=255) tCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg > 1:327927/1‑140 (MQ=255) tCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg > 1:605843/1‑140 (MQ=255) cGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa > 2:455645/1‑140 (MQ=255) cgcgTAGTGGGATAACCC‑AATTCCCCCTTGCTATCATCGGCGTTGAAGCGTTTGATGTTTGCAGAAATTGCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt < 2:568342/140‑1 (MQ=255) cgcgTAGTGCGATAACCC‑AATTTCCCCTTGCTATCATCAGCGTTGAGGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt < 1:729122/140‑1 (MQ=255) aGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTAATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACAt < 2:290463/140‑1 (MQ=255) gTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 1:154289/140‑1 (MQ=255) gTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc > 2:286133/1‑140 (MQ=255) gtgcaaaaCCCAAATTGCTCCTTGCTATCACAAGCGATGAAGCGTTCAACCTTCGCATCACCTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:368701/134‑1 (MQ=255) gCGATACCCC‑AATTTCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa < 2:327927/140‑1 (MQ=255) gATAACCC‑AATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGa < 2:349892/140‑1 (MQ=255) aCCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataat < 1:465243/140‑1 (MQ=255) cc‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAc < 1:563052/140‑1 (MQ=255) c‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt > 1:49558/1‑123 (MQ=255) c‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt < 2:49558/123‑1 (MQ=255) c‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAcc > 2:602965/1‑140 (MQ=255) aTTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGAc < 2:599277/124‑1 (MQ=255) aTTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGAc > 1:599277/1‑124 (MQ=255) cTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACa > 1:706844/1‑140 (MQ=255) cTTCCTATCATCAGCGTTAAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACa < 2:706844/140‑1 (MQ=255) ttGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACAt < 2:762162/140‑1 (MQ=255) gCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATcc < 1:534016/140‑1 (MQ=255) aGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAatat > 1:642083/1‑140 (MQ=255) cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTg > 2:414913/1‑140 (MQ=255) cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTg > 2:129378/1‑140 (MQ=255) gAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTg < 1:587442/140‑1 (MQ=255) aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGc > 1:474968/1‑140 (MQ=255) aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGc > 2:379636/1‑140 (MQ=255) aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGc < 1:602965/140‑1 (MQ=255) tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactac > 1:85555/1‑140 (MQ=255) tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactac > 2:734055/1‑140 (MQ=255) ttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCTTAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCg > 2:419331/1‑140 (MQ=255) | TTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCC‑AATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCG > NC_000913/1576986‑1577251 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |