Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,035,587 | 0 | C | G | 26.5% | 49.8 / 10.6 | 35 | intergenic (+344/‑248) | yedS/hchA | pseudogene, outer membrane protein homology; putative outer membrane protein/glyoxalase III and Hsp31 molecular chaperone |
Reads supporting (aligned to +/- strand): ref base C (14/11); new base G (9/0); total (24/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.72e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCA > NC_000913/2035521‑2035654 | cTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAAtctct > 1:741668/1‑68 (MQ=255) cTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAAtctct > 1:571925/1‑68 (MQ=255) cTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAAtctct > 1:441732/1‑68 (MQ=255) ttttCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTcc < 1:637015/68‑1 (MQ=255) tACCAATCCGGTATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAc < 2:486692/68‑1 (MQ=255) cAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAcccc < 2:427633/68‑1 (MQ=255) cAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAcccc < 2:532485/68‑1 (MQ=255) tCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTa > 1:788326/1‑68 (MQ=255) cGGGATTTAAATATTCATCAGCACCACAATTTTCCACTGCTCTAATTCTCTCCCCGCCACCCCGTAcc < 1:780646/68‑1 (MQ=255) cGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTAcc < 1:413794/68‑1 (MQ=255) atatAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTg > 1:357893/1‑68 (MQ=255) aaTATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAAt > 1:39547/1‑68 (MQ=255) tatTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATgg < 2:741668/68‑1 (MQ=255) atTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTAcc < 1:314605/56‑1 (MQ=255) atTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGACACCCCGTAcc > 2:314605/1‑56 (MQ=255) tcatcaGCACCACAATATTCCACAGCACTAAATCTGTGCCCGCCACCCCGTACCTCTGATAATGGTCt > 2:738120/1‑68 (MQ=255) aGCACCACAATATTCCACAGCACTAAATCTGTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAAt > 2:474898/1‑68 (MQ=255) acaATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCAt < 2:5163/65‑1 (MQ=255) acaATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCAt > 1:5163/1‑65 (MQ=255) acaATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGa > 1:270594/1‑68 (MQ=255) aaTATTCCACAGCACTAAATCTCTCCCAGCCACCCCGTACCTCAGATAATGGTCAAAAATCATTGAAg > 1:604797/1‑68 (MQ=255) atTCCACAGCACTAAATCTCTGCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCa > 2:599194/1‑68 (MQ=255) tCCACAGCACTAAATCTGTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACt > 2:193949/1‑68 (MQ=255) ccACAGCACTAAATCTGTCCCAGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACtt > 2:687606/1‑68 (MQ=255) ccACAGCACTAAATCTGTCCACGCCACACCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCAcat > 2:405185/1‑66 (MQ=255) ccACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACtt > 1:407531/1‑68 (MQ=255) ccACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACtt > 1:406395/1‑68 (MQ=255) ccACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATAGAAGCCACtt > 2:1029758/1‑68 (MQ=255) acaGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGc > 2:1019666/1‑68 (MQ=255) acaGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGc > 2:457421/1‑68 (MQ=255) aGCACTAAATCTGTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcga > 1:834392/1‑68 (MQ=255) gCACTAACTCTGTCCCCGCCACCCCGTCCCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgac > 1:177459/1‑68 (MQ=255) aCTAAATCTGTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga > 2:902445/1‑68 (MQ=255) cTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCgacgac < 1:29585/68‑1 (MQ=255) aaTCTCTCCCCGCCACCCCGTACCTCTGATAATGTTCTAAAATCATTGAATCCACTTGCGACGACGtt < 2:604797/68‑1 (MQ=255) aTCTGTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTg > 1:86222/1‑68 (MQ=255) aTCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTg > 1:390402/1‑68 (MQ=255) aTCTATACCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTg > 1:696749/1‑68 (MQ=255) tctctgCCCGCCACCCCGTACCTCTTATAATTGTCTAAATTCTTTGTAGCCACTTTCGACGACGtttc < 2:596813/68‑3 (MQ=255) ctcCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCa < 1:457421/68‑1 (MQ=255) ctcCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCa < 2:792031/68‑1 (MQ=255) ctcCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCa < 1:193949/68‑1 (MQ=255) | CTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCA > NC_000913/2035521‑2035654 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |