Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,601,176 | 0 | T | C | 45.5% | 53.6 / 19.6 | 45 | E304G (GAA→GGA) | ftsX | putative ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base T (21/3); new base C (0/20); total (21/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.02e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCT > NC_000913/3601029‑3601314 | tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTcc < 1:192993/150‑1 (MQ=255) tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCg > 1:530004/1‑151 (MQ=255) tCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCa > 1:129310/1‑150 (MQ=255) cAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCaa < 2:155292/150‑1 (MQ=255) aGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAAcc < 2:34274/151‑1 (MQ=255) gCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAAt > 1:112736/1‑151 (MQ=255) tAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCgc > 1:698862/1‑151 (MQ=255) aaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCgc > 2:148353/1‑149 (MQ=255) aCCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGa > 1:89447/1‑150 (MQ=255) cgggggCAACCCCCGCCCCCCCCCACCAAATCATCGAGCATCCCAGCAGAAAGAGCAGGCTTT‑CATCAAATGATAAGCCATTAATAACAACCTTC‑GTCCCAAAAC‑CCTGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAAt < 1:34611/147‑1 (MQ=255) aaccccccctcccccccACCAAACCATCGAGCATCCCAGCAGAAAAAGCAGGCATT‑CATCAAATAAAAACCCATTAATATCAACCTTC‑TTCCCAAAAC‑CCTGGCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAAtt < 1:750660/131‑1 (MQ=255) cccccccccccccccacccaATCATCGCGCACACCAGCAAAAATAGC‑GGCATTCCACAAAATAATAACCCATTGATA‑CAAACCTT‑GTTCCCAAAACCCTGTGCCCCTCCGCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaa < 1:688183/139‑1 (MQ=255) ccccccccccACCCCGCCAATCACCAACCATCCCACAACAAATAGCCGCCATT‑ATTTAATTAATAACCAATTGATAACAAGCCTC‑GTCCCAAAAC‑CCTGTGCCATTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaaa < 2:200560/143‑1 (MQ=255) cccccAGCAAAACATCGAGCACCCAAGCAGAAATAGCAGGCTTT‑CATAAAATAATAAGCTTTTGATATAAACCTTC‑GTCCCAAAAC‑CCTGTCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGac < 1:465296/149‑1 (MQ=255) cccccACCAAATCACCGAACATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCAATTAATATCAAACTTC‑GTCCAAAAAC‑CCTGCCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTaa < 1:572441/145‑1 (MQ=255) catcGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGcc > 1:276350/1‑150 (MQ=255) atcGAGCATACCAGCAGAAAAAGCGGCCATT‑AAGCAAATAATAAGCTTTTAATAAAAACCTTC‑TTCCCAAAAC‑CCTGCCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgcg < 1:459070/147‑1 (MQ=255) agcatccaagaaaaaaTAGCAGGCATT‑CATCAAATGATAACCCATTAATATCAATCTTC‑GTCCAAAAAA‑CCTGTCCCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAg < 2:759362/133‑1 (MQ=255) tACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATccc > 1:737468/1‑151 (MQ=255) caagcagaAATAGCGGATATT‑CATCGAATGAAAAGCTTTTGA‑AACAAACCTCTGTCCCGAACC‑CCGTTTCCCCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaaaa < 1:259037/143‑1 (MQ=255) agcagcAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcagc > 1:813187/1‑150 (MQ=255) gctttcagaaaaTGATAAGCATTTGAAATCAACCTTC‑TTTCAAAAAC‑CCGGGCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATa < 1:378513/140‑1 (MQ=255) cGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACgg > 1:280835/1‑151 (MQ=255) gAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACggg > 2:649568/1‑151 (MQ=255) aataatAACCCATTAATAACAATCTTC‑TTCCCAAAAC‑CCTGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACg < 1:243548/148‑1 (MQ=255) aataaaacgCCATTAATATAAATCTTC‑‑TTCCAAAAACCCTGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgc < 1:340049/143‑1 (MQ=255) aaTGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgc > 1:216965/1‑151 (MQ=255) aaTGAAAACCCTTTAATATCAACCTTC‑‑TTCCAAAAACCCTGTCCCCCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACg < 1:622739/148‑1 (MQ=255) aTGATAAGCCATTAATATCAACCTCC‑GTCCCGAAAC‑CCTGCCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcg < 2:339818/151‑1 (MQ=255) aggcaaTTAATACCAAACTCC‑GTTCCAAAAA‑CCTGCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGa < 2:839234/146‑1 (MQ=255) aaGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAgg > 2:298051/1‑150 (MQ=255) atatCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCAt > 1:669276/1‑151 (MQ=255) aacctccgtccaaaaac‑CCTGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTg < 2:153242/136‑1 (MQ=255) aaCTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCg > 2:221341/1‑151 (MQ=255) aaCTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCg > 2:826263/1‑151 (MQ=255) aCTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGc > 2:795821/1‑151 (MQ=255) gtcccaaaac‑CCTGTCCCCCTCCCGCAGCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCaa < 1:837286/145‑1 (MQ=255) ttCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAAt > 2:560677/1‑151 (MQ=255) ccaaaac‑CCTGTCCCACTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAAt < 1:365896/146‑1 (MQ=255) cGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATc > 1:85774/1‑149 (MQ=255) gAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt > 2:136168/1‑151 (MQ=255) aaac‑CCTGTTCCACCACCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCaa < 2:213284/145‑1 (MQ=255) aacccTGTGCCCCTCCCGCAACCGCCGATGACAATCGTAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt < 1:614970/148‑1 (MQ=255) aa‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTc > 1:180973/1‑151 (MQ=255) ccTGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCt < 1:763872/150‑1 (MQ=255) | TATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATT‑CATCGAATGATAAGCCATTGATATCAAACTTC‑GTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCT > NC_000913/3601029‑3601314 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |