Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 29,486 | 0 | C | A | 37.8% | 80.6 / 21.9 | 45 | intergenic (+291/‑165) | dapB/carA | dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase |
Reads supporting (aligned to +/- strand): ref base C (24/4); new base A (0/17); total (24/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.43e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT > NC_000913/29336‑29635 | tttttGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCcac > 2:124369/1‑151 (MQ=255) ttttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAaataat > 2:661046/1‑150 (MQ=255) tttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAataata < 1:1093/150‑1 (MQ=255) aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATaaa > 1:577989/1‑150 (MQ=255) gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATaaaaa > 2:18715/1‑150 (MQ=255) aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATccc > 2:433988/1‑151 (MQ=255) aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATccc > 1:613654/1‑151 (MQ=255) tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCg < 1:124332/150‑1 (MQ=255) ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg > 2:793037/1‑150 (MQ=255) cctcTTTTTTCTCCTTTTTTTTATTTCATCCAATCTTCTTGCTGCGCAAGCTTTTCCAAAAACGG‑TTT‑AAATTATTTTTTTGTCCCTAAATCCCTTAAAAACATGCA‑TGACCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAc < 2:139244/146‑1 (MQ=255) ccATTTTTTTCTTTTTTTTTTTATTTCATGCAATCTTCTTGCTGCGAAAGC‑TTTTCCAAAAACG‑GTTAAAATAATTTTTTTGTCTTTTAATCCCTTTAAAACATGCA‑TGAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACt < 1:282949/150‑1 (MQ=255) ttttttttGCTCTTTTTTTTTTAATGCATTTTTTTTGCTCCAAATGCTTTTTCAAAAAAAG‑GTA‑AGATAATTTTTTTGTCGTTTATTCCCGAAAAAAGTTGCT‑TGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACtttt < 1:169831/149‑1 (MQ=255) ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg > 2:52126/1‑151 (MQ=255) gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc > 1:279304/1‑151 (MQ=255) tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc > 1:754187/1‑150 (MQ=255) ttttATTTCATGCATTTTTTTTGCTGCGAAAGCTTTTCCAAAAACAG‑GTT‑AGATATTTTTTTTGCCGTTTATTCCCTAAAAAACATGAA‑TGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc < 1:784151/150‑1 (MQ=255) tCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc > 2:780930/1‑151 (MQ=255) aTGCATTCTTCTTGCGGCG‑AAGGCTTTTACAAAAAAGTGTT‑AAAAAATTTTTTTGACGTTTAGGCCCTAAAAAACGTGTT‑GGAGCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc < 1:660987/149‑1 (MQ=255) tGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc > 2:231453/1‑151 (MQ=255) ttcttGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc > 2:132716/1‑151 (MQ=255) gctgcGCAAGCGTTTTCCAGTACAG‑GTT‑AGATGATCTTTTTGTGGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa > 1:96022/1‑149 (MQ=255) tgcGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt > 1:359277/1‑151 (MQ=255) aagcttttcaaaaaaCAG‑TTT‑AAATATTTTTTTTGTCGTTAATTCCCTGAAAAACATTCA‑TGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc < 1:158349/137‑1 (MQ=255) aaGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg > 1:750671/1‑151 (MQ=255) aaGCGTTTCCAAAAACGG‑GTT‑AAATGATCTTTTTGTCGTTTAATGCCTGAAAAACATGCA‑TGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt < 2:565673/149‑1 (MQ=255) tttCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc > 1:94899/1‑151 (MQ=255) ttCCAGAACAG‑GTT‑AGATGATCTTTTTGGCGCTTAAAGCCTGTAACACCTGCA‑TAGGCCACCAAATACTCTAAAAAATCCCGCCATTAAGTTAACTGTTAGCGCCCATATCTCCCGAATGCCGCCGTTTGCAACAAATTCGTCGGTAAGc > 1:529128/1‑150 (MQ=255) ccAAAACAG‑GTA‑AGTTGATTTTTTGTTCTTTTATGGCCTTAAAAAATTGCA‑GGAACAAAAAAATAAAAAAAAACACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa < 1:354375/149‑1 (MQ=255) aaattatttttttGTCGCTTAATGCTTAAAAAACAGGAA‑TGAGCCAAAAAATAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa < 1:106500/141‑1 (MQ=255) atgtttttttttGCGCTAAATGCCTTAAAAACAGGTA‑GGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa < 2:485951/141‑1 (MQ=255) taatttttttGGCGAATATTGCCTTAAAAACAAGCA‑GGACCAAAAAAAAAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa < 1:708301/145‑1 (MQ=255) tttttttGTCATTTATTCCTTGAAAAACATGAA‑TGAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt < 2:790300/148‑1 (MQ=255) tCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGTAGATTTGCATTGATTTACGtc > 2:423739/1‑151 (MQ=255) ttttttGTCGTTTAATCCTTGAAAAACTTGCA‑TGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc < 2:614755/149‑1 (MQ=255) tttttGTCGTTTAATGCCTGAAAAACATGCA‑TGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca < 2:776181/150‑1 (MQ=255) tttGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc > 2:455373/1‑150 (MQ=255) tGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt > 1:106089/1‑151 (MQ=255) ttAAAACTTTCATTGACCCAAAAAATAATATAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg < 1:45468/149‑1 (MQ=255) aaacaatgct‑tgAGCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa < 1:11599/146‑1 (MQ=255) aCAGGCA‑GGACCAACAAAATAATAAAAAAAAACCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa < 1:540004/150‑1 (MQ=255) catgca‑tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa > 2:435888/1‑150 (MQ=255) catgca‑tgAACCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa < 2:105881/150‑1 (MQ=255) ca‑tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg > 1:680790/1‑150 (MQ=255) a‑tgACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtga < 2:297622/150‑1 (MQ=255) ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGTTTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga > 1:511518/1‑149 (MQ=255) ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat > 2:301258/1‑151 (MQ=255) acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat > 2:681467/1‑151 (MQ=255) | TTTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT > NC_000913/29336‑29635 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |