| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 298,960 | 0 | A | C | 46.8% | 62.9 / 14.9 | 48 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base A (21/4); new base C (0/22); total (21/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.47e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCA > NC_000913/298814‑299102 | ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATAcccc > 2:182923/1‑151 (MQ=255) cgcgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCat > 2:229869/1‑150 (MQ=255) gATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGAtt > 2:866089/1‑150 (MQ=255) gATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGAtt > 2:615698/1‑150 (MQ=255) gATAGCCGCGCTCACGCCGCACAGCCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCCCCTGTTTGGGG‑ATATCCGCATGAC‑CCGGCCCC‑TCATCCCCCGCCATATCGTGATTa < 2:908376/151‑1 (MQ=255) ccccgccaaCGCCACACAGCCCCGGCCCCCCG‑CCCCTTTGGCCCTCCCTGGGGGAGGATAGGGGTTCGTTACCA‑CCAAGAAAAACCCCCCCGTTTTTGGGAATACCCGCGTAAC‑CCGGCCCC‑TAAACCCCCGCCATATCGTGATTAAc < 1:424448/140‑1 (MQ=255) ccGCGCTCACGCCGCACAGGCCCAGCCCCCCGACCCTTTTGC‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCACCTCCTTTTTTGGG‑ATACCCGCATGAC‑CCGGCCCC‑TCATCCCCCGCCATATCGTGATTAAc < 2:478390/148‑1 (MQ=255) ccGCCCTCACGCCGCACGGCCCAAGCCCCCCGCCCCTTTTGG‑CCTTCAGGGGGGAGGATAGGGGTTGGTTATCACCCAGAAAAACCCCCTCCTGTTTGGGA‑TTCCCCGAATGAA‑CCGGCACC‑TCACACCCCGCCATATCGTGATTAACGaa < 1:719470/151‑1 (MQ=255) cgcTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGt > 1:257290/1‑150 (MQ=255) gctacccccgaACAGCCCAAGCACACGAACACTTTTGC‑CCTCAATCGGGGGGAATGGGGGGTCGATATCACCCAGAAAAACCCCCTCCTTTTTGGGA‑ATATCCGCATAAC‑CCGGCCCC‑TCATACCCCGCCATATCGTGATTaa < 1:179194/136‑1 (MQ=255) cTCCCGCCGCACGGGCCCAGCTCACCGCCACTTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCCCCCCCTTTTTTGGG‑ATATCCGCATGAA‑CCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAg < 2:703766/150‑1 (MQ=255) cTCCCGCCGCACAGGCCCAGCTCACCGACCCCTTGGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAC‑CCGGCACC‑TCATCCCCCGCCATATCGTGATTAACGAAGt < 2:922057/148‑1 (MQ=255) cACCCCGCACGGCCCTAGCTCCCCGACGCCTTTGG‑CCTTCACGGGGGGGAATGTGGGGTGGGTATCACCCAGAAAACCCCCCCCCTTTTTGGGA‑AAATCGGAAGGAA‑CCGGACCC‑TAATCCCCCGCCATATCGTGATTAACGAAGTAg < 1:697152/148‑1 (MQ=255) ccgcACAGCCCCAGCTCCCCGACCCCTTTGC‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaaa < 2:419421/150‑1 (MQ=255) ccgcACAGCCCCAGCCCCCCGCCCCTTGTGG‑CCTTCGGGGGGGAGGATTGGGGGTCGGTATCACCAAGAAAAACCCCCCCTTGTTTGGGA‑ATATCCGCATAAC‑CCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa < 1:339910/148‑1 (MQ=255) cAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAcc > 2:913659/1‑151 (MQ=255) cAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAcc > 1:16643/1‑151 (MQ=255) gCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc > 1:846292/1‑151 (MQ=255) gCTCACCGACACCTTCGG‑CATTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc > 1:227779/1‑151 (MQ=255) ccccccgcccccTTTGG‑CTTTCATGGGGGGGGATATGGGGTCGTTATCCCCAAGGAAAATCCCCCCTTGTTTGGGG‑ATACCCGCATGAC‑CGGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc < 2:691583/141‑1 (MQ=255) ccctccacTTTTGC‑CCTTCGCCGGGGACGATGGGGGGTCGGTCTCAACCAGGAAAGCCCCCTCCTGTGTTGGC‑AACCCCCCACAAA‑CCGGCGCT‑TCCTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc < 1:256930/143‑1 (MQ=255) cccaccacCTTGGG‑CTTCCGTGGGGGGGGATAGGGGGTGGATA‑CAACCCAGAAAACCCCCCCCTTTTTTGGGAATAACC‑AATGAACCCGGCCCA‑GCAACCCCCGCCATATCGTGATTAACGAAGTAGccc < 1:103630/125‑1 (MQ=255) acCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc > 2:184796/1‑151 (MQ=255) ccTTTGG‑CCTTCATGGGGGAGGATATGGGGTCGGTACCATCCAGAAAAACCCCCTCTTTTTTGGGG‑ATACCCGCATAAC‑CCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc < 2:431416/150‑1 (MQ=255) tGG‑CCTTCATCGGGGGGGATATGGGGTCGTTATCATCCAGGAAAATCCCCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCCCC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc < 1:746552/149‑1 (MQ=255) tGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc > 2:137146/1‑149 (MQ=255) gtcctttctaaggggAGTATAGGGGGTCATTATCACCCAGAAAAACCCCCTCTGGTTTTGTA‑ATTCCCGGATAAC‑CCGCCTCT‑TCCTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc < 1:526316/138‑1 (MQ=255) gggAGGTTATGGGGCGGGACTCCCCAAGAAAAATCCCCTCCTTTTGTGGG‑ATACCCG‑ATGAACCCGCCCCC‑TCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat < 1:308953/151‑1 (MQ=255) gggAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTTTGGA‑ATACCCGCAGAAC‑CCGCCCCC‑CCATCCCCCGCCATATCGTGATTAACGAAGTGGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat < 2:257788/151‑1 (MQ=255) gaggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgc > 2:447433/1‑151 (MQ=255) atatGGGGTCGGTGTAATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgaccatgccca > 1:698686/1‑150 (MQ=255) tatGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat > 2:634492/1‑150 (MQ=255) tagggggTGGTTATCACCCAGAAAAACCCCCCCCTGTTTGGGG‑ATATCCGCAGGAC‑CCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat < 2:56001/147‑1 (MQ=255) gggggTGGGTACCACCCAGAAAAACCCCCTCCTTTTTGGGG‑AAACCCCAATAAC‑CCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt < 2:778120/150‑1 (MQ=255) ggTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg > 2:581415/1‑151 (MQ=255) atcatcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc > 2:813781/1‑150 (MQ=255) aGGAAAACCCCCTCCTGTTTTGGA‑ATACCCGCATAAC‑CCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAcc < 2:775002/150‑1 (MQ=255) gAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTgg > 2:335517/1‑151 (MQ=255) aaaaCCCCCTCCTGTTTTGGA‑ATATCCGCATGAA‑CCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg < 2:890994/149‑1 (MQ=255) accACCCCCTGTTGTGGA‑ACCTCGCCATGAG‑CCGGCCCC‑TAAGGCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCACCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgc < 2:710825/149‑1 (MQ=255) tCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg > 1:771282/1‑151 (MQ=255) tttgGGG‑ATACCCGCA‑GAA‑CCCGCCCCTTCAACCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGt < 1:848485/147‑1 (MQ=255) tgggtatatCCCCATGAC‑CCGCCCCC‑CCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGt < 1:452275/140‑1 (MQ=255) cATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCg > 2:876391/1‑151 (MQ=255) aa‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt > 1:299898/1‑151 (MQ=255) tggCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc > 1:190939/2‑150 (MQ=255) ccgcccccccATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc < 2:375634/138‑1 (MQ=255) ccc‑TCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagca < 1:700972/149‑1 (MQ=255) | ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAA‑CCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCA > NC_000913/298814‑299102 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |