Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,601,160 | 0 | G | A | 40.7% | 31.8 / 11.0 | 28 | F309F (TTC→TTT) | ftsX | putative ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base G (11/5); new base A (0/11); total (11/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.49e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAAT > NC_000913/3601017‑3601289 | acacGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑t < 2:270676/151‑1 (MQ=255) acGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑tgt < 1:48343/151‑1 (MQ=255) gCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGc < 1:165790/151‑1 (MQ=255) tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCg > 1:476500/1‑151 (MQ=255) ggCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAAccg < 2:76605/151‑1 (MQ=255) tGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGAc > 2:36358/1‑150 (MQ=255) gCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAAt > 1:145184/1‑151 (MQ=255) ttgtACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCAcc > 1:513618/1‑151 (MQ=255) caCCGTGGCAACCCACCCCGCCCCCCACCAAACCATAGAGAAGACCAGCACCAACAGCAGGCATTCACCAAATAAAAACCCTTTAATATAAACCTTCTTCCCAAAAC‑CC‑TGTCCCACTCCCGAAACCGCCGATGACAATCGCAGCACCa < 2:341175/148‑1 (MQ=255) gcggaaacccACGCCCCCCCCCAGCAAATCAGCGAGCAAACCGGCAGCAATAGAGGCCATCCAACAAAAAAAAACCAATTGATACAAAACTCCTTCCCAAAAC‑CC‑TGTTCCTCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAAtt < 2:579372/142‑1 (MQ=255) gTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAAttt > 1:270316/1‑151 (MQ=255) cccccctcccccccacccaATCACCAAGCATCCCAGAAGCAATAGCAGGTATCCATAGAATGATAACCTTTTAATAAAAACCTTC‑TTCCAAAAACCC‑TGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACa < 1:122201/138‑1 (MQ=255) cgcccccccccACCAAATCATCGACCAAACCAGCAAAAATAGAAGCCATCCACAGAAT‑AAAAACCATTTAAAAAAAACCTTCTTCCAAAAACCC‑GGTCCCCCTCCAACATCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAt < 2:350399/142‑1 (MQ=255) aGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGAcaacaa < 2:372136/151‑1 (MQ=255) caaacaatcGACCATCCCAACAGAAATAGCAGGCATTCATCAAATGAAAACCCATTAATATCAACCTCCTTCCCAAAAC‑CC‑TGTCCCCCTACCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaaca < 2:139356/142‑1 (MQ=255) cacagAGCATACCGGCAGAAATAGAATCCATTAAGCGATTAAAAAGCCGTGAATAACAACCTCCGTCGCAAAAC‑CA‑TGGTCCTGTCCCGAAACCGCCGATGACAATCGCAGCACCATAATTTCTGACAAAATTAATGacaaca < 2:464369/139‑1 (MQ=255) aCCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATccc > 2:353949/1‑150 (MQ=255) aaaaaacgcggccATTCACCGAATGAAAACCCATTGAAATCAAACTTCTTCCCAAAAC‑CC‑TGTCCCACTCCAAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcagcag < 2:45591/137‑1 (MQ=255) aaaaaTAGCAGCCTTTCATCGAATAAAAACCCATTAATAACAACCTTCGTCCCAAAAC‑CC‑TGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcagcag < 1:171421/146‑1 (MQ=255) cATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATAc > 1:36478/1‑150 (MQ=255) aTTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATAc > 2:509922/1‑149 (MQ=255) tCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACg > 2:260409/1‑151 (MQ=255) aaagaaaaGCCATTAAAAACAACCTCCGTTCCAAAAC‑CC‑TGTCCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGaa < 1:433092/140‑1 (MQ=255) tGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcgc > 1:111229/1‑151 (MQ=255) aaccatttaatctaaaaCTTCGTCCCAAACC‑CC‑TGCCCCACTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGa < 1:315770/137‑1 (MQ=255) aaGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAg > 1:418150/1‑149 (MQ=255) ttaataaCAAACTCCTTCCCAAAAC‑CC‑TGCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAAt < 1:43192/147‑1 (MQ=255) tctatacAAAATTCATTCACAAAA‑CCGTGCCCCACTCCAACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGa < 2:523650/142‑1 (MQ=255) | ACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAAT > NC_000913/3601017‑3601289 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |