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.7% | 34.2 / 17.1 | 35 | E304G (GAA→GGA) | ftsX | putative ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base T (14/5); new base C (0/16); total (14/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.65e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGG > NC_000913/3601026‑3601326 | tattatTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACtt > 2:902560/1‑151 (MQ=255) attattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTc < 1:728753/151‑1 (MQ=255) tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTcc < 2:846672/150‑1 (MQ=255) tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCg > 1:740050/1‑151 (MQ=255) aGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAAcc < 2:380003/151‑1 (MQ=255) aaGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATg > 2:471451/1‑151 (MQ=255) cGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATc > 1:630118/1‑151 (MQ=255) aGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGAc > 2:638933/1‑150 (MQ=255) cccgcccccccccAGTCAATCATCGACAACCCCAGCAGAAATAGAGGCCATCCAGAAAATAATACGCCTTTGAT‑CCAAACCTCCCCCCCAAAA‑CCTT‑GTCCCTCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaa < 1:346108/139‑1 (MQ=255) aCGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAt > 2:218915/1‑151 (MQ=255) ctccccccaacccAATCACCGACCATCCCAACAGAAATAGCAGCCATTCACCAAATAAAAAGCCATTGATA‑AAAACCTCCTTCCCAAAACCC‑T‑GTCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTa < 1:870487/139‑1 (MQ=255) tGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAAt > 2:145179/1‑151 (MQ=255) cccccccACCAAATCATCGAACAAACCAGCAGAAATAGCAGCCATTCATAGAATGATAAGCCTTTGAAAAAAATCTTCG‑TCCAAAAA‑CCGT‑GCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAAt < 2:823562/146‑1 (MQ=255) ccccccacccaATCATCGAGCATACCAGCACCATTAGCAGGCTTTCATCAAATAATAAGCATTTAAAAACAACCTCCTTCCCAAAAA‑CC‑T‑GTCCCCCTCCCGAATCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGa < 2:841340/148‑1 (MQ=255) ccccccACCAAATAATCGAGCATACCAGAAGAAATAGAAGCCATTAAAAGAATAAAAAGCAATAAATACAAATCTTCGTCCCAAAAC‑CC‑T‑GTCCCTCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGa < 2:583323/148‑1 (MQ=255) cAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgc > 2:601141/1‑151 (MQ=255) accGACCATACCAGAAAAAATAGCAGCCATTAATAAAATAAAAAGCATTTAAAAAAAAACTTCGTCCCAAAAC‑CC‑T‑GTCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaaca < 1:708244/140‑1 (MQ=255) accGACCATACCACCAAAAAAAGCAGGCTTTCATCAAATAATAACCTTTTAATATCAATCTTCTTTCCAAAAC‑CC‑T‑GTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCa < 2:857133/148‑1 (MQ=255) aGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaaa > 1:461096/1‑150 (MQ=255) aTTAATCAAAGGATATCCCATTGAAATAAATCTTCTCCCCGAAAC‑CC‑G‑GCCCCCCTTCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATa < 1:594565/148‑1 (MQ=255) aTCAAATGATAAGCCTTTGATAAAAACCTTCGTCCAAAAAC‑CC‑T‑GTCCCCCTCCCGAATCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCAt < 2:687736/143‑1 (MQ=255) gAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACgg > 1:798/1‑150 (MQ=255) aataatAAGCTATTGATATAAAACTTCTTCCCAAAAA‑CC‑T‑GTCCCACTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACg < 1:114521/148‑1 (MQ=255) aaTGAAAAGCCATTAATACAAAACTCCGAACCGAAAC‑CC‑T‑GACCCACCTCCAAACCCGCCGATGACAATCGCGGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACggg < 1:930096/150‑1 (MQ=255) tGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcgc > 2:608238/1‑151 (MQ=255) aTAAGCAATGGATACAAAACTCCCTCCAAAAAC‑CC‑T‑GTCCCACCCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgc < 1:751523/147‑1 (MQ=255) agcatttGATAAAAACCTTC‑TTCCCAAAA‑CC‑TGGTCCCCCCCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt < 2:738100/146‑1 (MQ=255) ataaCAACCTTCGTCCAAAAAC‑CC‑T‑GCCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCAt < 1:780654/151‑1 (MQ=255) aaaCTTCTTCCCAAAAA‑CC‑T‑GTCCCACCCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGt < 1:543374/150‑1 (MQ=255) ccttcttccaaaaac‑CC‑T‑GTCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGc < 2:588246/140‑1 (MQ=255) gAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt > 2:560677/1‑151 (MQ=255) aaa‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt > 2:254182/1‑151 (MQ=255) acccc‑T‑GTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt < 1:572768/148‑1 (MQ=255) aaccc‑T‑GTCCCCCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt < 2:707698/148‑1 (MQ=255) t‑gtgCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGt > 2:169009/1‑150 (MQ=255) tgCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACg > 1:646847/1‑151 (MQ=255) tCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATgg > 2:336501/1‑151 (MQ=255) | TATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGG > NC_000913/3601026‑3601326 |
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: ‑ |