| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,029,274 | 0 | C | T | 44.9% | 36.1 / 25.3 | 49 | K31K (AAG→AAA) | yedI | DUF808 family inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base C (11/16); new base T (22/0); total (33/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.29e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATGATCAGCGCCAGCGGCACCAGAATCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGCAAAAACGGCAGCTCTGAT > NC_000913/2029125‑2029424 | atgatCAGCGCCAGCGGCACCAGAATCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttct < 2:535138/151‑1 (MQ=255) tgatCAGCGCCAGCGGCACCAGAATCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttctt > 2:134156/1‑151 (MQ=255) atCAGCGCCAGCGGCACCAGAATCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTcg < 2:445991/151‑1 (MQ=255) gcCAGCGGCACCAGAATCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTgccgcc > 2:119249/1‑150 (MQ=255) cAGAATCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTGCCGCCCG‑TTTGccc > 1:672546/1‑148 (MQ=255) cAGAATCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTTGCCGCCCG‑TTTGccc > 2:664362/1‑148 (MQ=255) aTCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACgg < 1:129168/151‑1 (MQ=255) cTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCCGCCG‑TTTGTCCATCAc > 1:423793/1‑145 (MQ=255) attaATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTCCCCCCCG‑TTTGCCCATCACCGGGata > 1:188629/1‑148 (MQ=255) attaATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATc < 2:662393/150‑1 (MQ=255) aTCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTGCCGCCCG‑GTTTCCCCTCACCGGGATATc > 1:231342/1‑146 (MQ=255) tCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCCA‑TTTGCCCCTCACGGCGATATCg > 2:556565/1‑146 (MQ=255) cAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCCTCACGGAGATATCGTCCagt > 1:192530/1‑151 (MQ=255) gATCCTTTCGCCACGCCCCAGACCACGGGAAGGTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGCTTTTTTCCGCGGCCG‑TTTTTCCCTCACCGAGCTAtcctcccgcag > 1:12455/1‑144 (MQ=255) tCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTCCCGCCCG‑GTTGCCCCTCACCGGGATATCGTCCCGGAGggtg > 1:579582/1‑151 (MQ=255) cACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTTGCCCCCCGTTTTGCCCA‑CACTGAGATATCTTACAGTAGTGTTGCGat > 2:845222/1‑147 (MQ=255) ccccAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCCCCAG‑TTTGCCCCTCACGGAGATATCGTCCCGTAGGGTGGCggaa > 1:70708/1‑141 (MQ=255) ccccAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGa < 2:287787/151‑1 (MQ=255) cccAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCAG‑TTTGCCCCTCACGGAGATATCGTCCAGGAGTGTGGCGATat‑cat > 2:408041/1‑147 (MQ=255) ccAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCAG‑TTTGCCCCTCACCGGGAAATCGGCCCGGAGTGTGGGGGTat‑catc > 1:62759/1‑147 (MQ=255) aCGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCCCCCG‑TTTGCCCATCACCGAGGTATCGTCCCGTAGTGTGGCGGTAT‑CATCGAGCAACggc > 1:25014/1‑148 (MQ=255) ccGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTGCCGCCCG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑TATCGAGCAACGTCAGTAAacta > 2:524112/1‑148 (MQ=255) gTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCCG‑TTTGCCCATCACCGAGATATCGTCCCGTAGTGTGGGGGTAT‑CATCCAGCAACGTCAGTAAAATACTccc > 1:333601/1‑150 (MQ=255) gTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGGTAT‑CATCGAGCAACGTCAGTAAACTACTccc > 2:667474/1‑150 (MQ=255) cGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCaaaa < 2:300008/151‑1 (MQ=255) gCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTGCCGCCAG‑TTTGCCCATCACGCAGATATTGTCCAGATGTGGGGCACTAT‑CATCGAGCCACGTCAGGAAACAATTCCCCGCCcaaa > 1:746795/1‑150 (MQ=255) cTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCCG‑TTTGGCCATCACCGAGATATCCTCCCGTAGTGTGGCGGTAT‑TATCGAGCAACGGCAGTAAACTACTCCCGGCCAAAAACATTcc > 1:177506/1‑151 (MQ=255) gAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCct < 1:117921/151‑1 (MQ=255) aaaCTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTccc < 1:560663/149‑1 (MQ=255) tgttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCCG‑GTTGGCCATCACGGAGAAAACTTCCCGTAGTTTGGGGGTAT‑CATCCAGCAACGTCAGTAAACTAATCCCGGCCAAAAGAAATCCCCCT‑CCtt > 1:876562/1‑151 (MQ=255) tgttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTTGCCGCCAG‑GTTGCCCATCACGGAGATATCGTCCAGGAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTTCCTCT‑CAtt > 2:124106/1‑151 (MQ=255) cgcATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTGCCGCCCG‑ATTGCCCCTTACCGAGATATTGGCCCGTAGTGTGGCCGTATCCTTCGA‑CAACCTCCGTAAACTAATCCCCGCCAAAAACAATTCCTCC‑CAttttttt > 2:811987/1‑149 (MQ=255) aTTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCCGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CCtttttttttt > 2:97222/1‑149 (MQ=255) aTTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTAttgt < 1:151044/151‑1 (MQ=255) tGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTGCCGCCCA‑TTTGGCCCTCAACGAGATATCGGTCAGTGGTGTGGGGGTAT‑CATCGAGCAACGGAAATAAAATACTCCCGGCCAAAAGCATTACCTCTCCAtttttttttt > 2:558160/1‑148 (MQ=255) gAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCCG‑TTTTCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CCTTTTTAttttt > 2:731920/1‑149 (MQ=255) gAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCCG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCCAAAACATTCCCTCT‑CATTTTTTttgtt > 2:756305/1‑149 (MQ=255) aGCGATAAGTCATCCCCTAATAGACCGAAGGTATTCTTCGCCGCCCG‑TATGCCCATCACGGAGATATCGTCCAGTAGTGTTTCGTTAT‑CATAGTGCGTCGTTAGTGACCTACTCCCGGCCAAAAGCATTCACTCA‑AATTCTTCTGGTTGa < 1:291147/150‑1 (MQ=255) aGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTTCCCCTCACGGAGATATCGTCCCGTAGTGTGGCGATAT‑CATCCAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATtttttttttt > 2:779579/1‑148 (MQ=255) aaGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCAG‑TTTGCCCCACACGGAGAAAACGTCCAGTAGTGGGGGGGTAT‑CATCGGGCAACGTCAGAAAACTACTCCCGGCCCAAAACATTCCCCCT‑CATTTTTATTGTTTTAggggg > 2:706478/1‑148 (MQ=255) aaGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGgtgt < 1:254311/151‑1 (MQ=255) aGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCCG‑GTTGCCCACCACGAAGATATAATCCAGTAGTTTCGCGCTAT‑TATCGAGCCACAGCCGCAAAATAAACCCCGCCAAAAGCAATTCCTCT‑CATTTTTATTGTTGAAggg > 1:220740/1‑147 (MQ=255) gTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTaa < 1:756231/151‑1 (MQ=255) tCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCCG‑TTTTCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCCGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTTAAggg > 2:321903/1‑145 (MQ=255) tCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTa < 2:653224/150‑1 (MQ=255) ccccTAATACACCGGCGGTTTTCTTCGCCGCCCG‑GTTTCCCATCACGGAGATATCGTCCAGTAGTGTGGCGGTAT‑CATCGAGCCACGTCAGTAAACTACTCCCGGCCAAAAGAATTCCCTCT‑CATTTTTTTTGTTTAAGGGTGGAAg > 1:194260/1‑147 (MQ=255) cccTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGa < 2:674311/151‑1 (MQ=255) aaTACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGc < 1:714661/150‑1 (MQ=255) cacCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGCaaaaa < 2:434924/151‑1 (MQ=255) ccGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGCAAAAAc > 1:433669/1‑150 (MQ=255) gcggTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGCAAAAACGGCa < 2:795523/151‑1 (MQ=255) ttcttcGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGCAAAAACGGCAGctct < 1:407984/150‑1 (MQ=255) ttcttcGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGCAAAAACGGCAGCTCTg > 1:868291/1‑151 (MQ=255) cttcGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGCAAAAACGGCAGCTCTGAt > 2:64831/1‑151 (MQ=255) | ATGATCAGCGCCAGCGGCACCAGAATCACTTTATTAATCAGCGATCCTTTCGCCACGCCCCAGACCACGGGAAGTTCCCGGTTGGCCCGCACGCCTGAAACTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAG‑TTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCGATAT‑CATCGAGCAACGTCAGTAAACTACTCCCGGCCAAAAGCATTCCCTCT‑CATTTTTATTGTTGAAGGGTGTAAGTATGAAGCAAAAACGGCAGCTCTGAT > NC_000913/2029125‑2029424 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |