| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,131,440 | 0 | T | C | 50.9% | 66.4 / 90.6 | 55 | G323G (GGA→GGG) | yghQ | putative inner membrane polysaccharide flippase |
| Reads supporting (aligned to +/- strand): ref base T (22/5); new base C (0/28); total (22/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.21e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-12 | |||||||||||
AAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTA > NC_000913/3131292‑3131584 | aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc < 2:21627/151‑1 (MQ=255)aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc > 2:294057/1‑151 (MQ=255)aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc < 1:247993/151‑1 (MQ=255)aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc < 2:715471/151‑1 (MQ=255) atgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgc > 2:444349/1‑151 (MQ=255) gTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAg < 2:943844/150‑1 (MQ=255) tCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGt > 1:500137/1‑150 (MQ=255) aGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTaa > 1:782376/1‑150 (MQ=255) aaaaaTTCCCCCTCTTATCCGCCCCTCCCCAAAATCAGGTCATCCGCTCCGAGATACTTCACGCCAACCCCCAGCAAAAGGGGCGGCTTCCCGCCGCTCACCCCCGGGCCCGCCCCCAGAATCCCGTTCCCGC‑CCGCCAGTAAACCGGAc < 2:42716/147‑1 (MQ=255) agaaTTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACtt > 1:489063/1‑151 (MQ=255) aaTTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTc > 1:104819/1‑150 (MQ=255) tttCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTcac > 2:413006/1‑150 (MQ=255) cccGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCt > 1:289732/1‑150 (MQ=255) cccGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCt > 1:608887/1‑150 (MQ=255) tatcccccccTCCCCAAAACCAGGTCATCCGCTTCGAGATACTAAACGCCAAACCCAAGCAAAAGGGGCGGCTTGCCGACGATCACCCCCGCGAGCGCCACCAGAATCCCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTaacaac < 1:449788/142‑1 (MQ=255) tcccccccTCCCCAAAATCGGTTCATCCGCTTGGAGATACTTACCCCCAACCCCCAGCAAAAGGAGCGGCTTCCCGACGACCACCCCCGCGCGCCCCACAAAAATACCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAAcc < 1:757474/143‑1 (MQ=255) ccccccTCCCCAAAATCGGTTCATCCGCTTCGAGATCCTAACCCCCAACCCCCAGCAAAAGGGGCGCCTTGCCGACAATCAACCCCGCGAGCCCCACAAAAATCCCTATCCCCC‑CCGCCAGTAAACCGGACTTCACACCTaacaac < 1:542356/142‑1 (MQ=255) cACCTCCCCAAAACCGGCTCCTCCGCTTCGAGATACTCACCCCAAAACACCAGCAAAGTGACCGCCCCCCCGCCAATCACCCCCGCGCCCCCCACCAAAACCCCTATCCCCC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTc < 2:886043/151‑1 (MQ=255) cttccccAAAATCAGGTCATCCGCTCCGAAATCTTTAACCCCAACCACCAGCGAAAGGAGCGGCTTGCCGACAATCAACCCCCCGGGCCCCCCAAGAATCCCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCt < 1:175360/146‑1 (MQ=255) ccTCCCCAAAATCAGGTCATACGCTTCGAAATACTAAACCCCAACCCCCAGAAAAAGGCGCGGCTCCCCGACGACCAACCCCGCGACCGCCCCCAAAATCCCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtg < 2:357744/151‑1 (MQ=255) ccTCCCCAAAACCAGGTAATCCGCTCGGAGATACTAACCGCCAAACACCAGCAAAATGAGCGGCTGGCCGACGACCAACCCCGCGAGCGCCACCAAAATACCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtg < 1:456626/151‑1 (MQ=255) cccTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtggtgg > 1:994854/1‑151 (MQ=255) aaTCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTgcgc > 2:932667/1‑151 (MQ=255) caagatCAGACGCTTCGATATTCTAATCACCAAACCCCAGCGAAATGCCGGCCTAGCGCAAGATCACCCCGGCGCGCGCCACAACAATACCACTCCCCC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCgg < 1:531711/146‑1 (MQ=255) ggTCATCCGCTCCGAAATACTTAACCCCAACCACCAGCAAAATGGGCGGCTTGCCGACGATCAACCCCGCGAGCCCCACAAAAATACCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGAt < 1:857241/150‑1 (MQ=255) atccGCTTCGAGATAATAACCCCCAAACCCAAGAAAAATGGGCGCCTCCCCGACGCCCACCCCCGCGAGCCCCACAAAAACACCGACCCCCC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTaaa < 2:664763/148‑1 (MQ=255) tCGAGATCCTCAACCCAAAACCCCAGAGAAGGGGGGGGCTTCCAGACGAACAACCCCGCGGGCGCCACAAAAATCCGGTTCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATc < 1:680479/150‑1 (MQ=255) tCGAGATCCAAAACACCAACCCCCAGCAAAGGGGGCGCCTCCCCGACGAACACCCCCGCGGGCCCCACCAAAACACCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATc < 2:359527/150‑1 (MQ=255) aaaTACTAACCCCCAACCACCAGCAAAAGGGGCGGCTTGCCCACAATCACCCCCGCGACCCCCACCAAAATCCCGACCCCCC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCAt < 2:88553/144‑1 (MQ=255) aaCGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTa > 1:742308/1‑150 (MQ=255) cGCCAACCACCAGCGAAATGGGCGGTTTGCCGACAACCAACCCCGCGACCCCCACAAAAACACCGATCCCGC‑CCGCCAGTAAACCTGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTa < 2:230859/148‑1 (MQ=255) gCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTaaaa > 1:121342/1‑150 (MQ=255) aacccccAGCAAAATGGGCGGCTTGCCGCCAAAAACCCCCGCGCGCCCCACAAGAAACCCGTTCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTggg < 1:94424/137‑1 (MQ=255) accccagaaaaaTGAGCGGCTCCCCGCCAACCACCCCCGCGACCCCCCCCAAAATCCCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAAc < 2:256853/138‑1 (MQ=255) acCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAc > 1:949536/1‑151 (MQ=255) ccAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAcc > 2:159940/1‑151 (MQ=255) agaaaaaTGGGCGGCTTGCCGACGATCAACCCCGCGAGCCCCACAAAAATACCTTTCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt < 2:575444/142‑1 (MQ=255) cGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGc > 1:780357/1‑151 (MQ=255) aagtgggCGGCTTCCCGACGAACAACCCCGCGAGCCCCACAAAAACCCCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCa < 2:242509/141‑1 (MQ=255) aaagggggggCTGGCCGACGACCACCCCCGCGGCCCCCACAAAAACCCCGCTCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCaa < 2:597925/143‑1 (MQ=255) aaaGGAGCGGCTTGCCGACGATCACCCCCGCGGGCCCCACAAGAATCCCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAc < 1:782351/144‑1 (MQ=255) aaaGGAGCGCCTTGCCGACGACCAACCCCGCGAGCCCCACAAGAATACCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCa < 1:259407/147‑1 (MQ=255) aaggAGCGGCTCGCGGACGACCACCACCGGGGGCCCCACCAAAATACCGTTCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCaaa < 2:627902/148‑1 (MQ=255) cgacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTcc > 1:827628/1‑151 (MQ=255) gacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcg > 2:554377/1‑151 (MQ=255) gacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcg > 1:334533/1‑151 (MQ=255) acgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcgg > 1:137557/1‑151 (MQ=255) cgaTCAACACCGCGGCCCCCCCCAAAATACCGATCCCCC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGtt < 2:741353/146‑1 (MQ=255) ccccGCGAGCCCCCCAAGAACCCCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTc < 2:423204/148‑1 (MQ=255) acCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCg > 1:779102/1‑150 (MQ=255) cccaccaAAACACGTATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCa < 1:397971/150‑1 (MQ=255) ccacaaAAATACCGCTCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATg < 2:919139/147‑1 (MQ=255) aaaaTCCCGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCg < 1:126339/148‑1 (MQ=255) aaaccccGATCCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAt < 2:489007/146‑1 (MQ=255) atcactgttCCGCGCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATg < 1:198749/138‑1 (MQ=255) ccGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTa > 1:202203/1‑150 (MQ=255) | AAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGC‑CCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTA > NC_000913/3131292‑3131584 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |