Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,982,529 | 0 | C | A | 41.1% | 98.3 / 35.9 | 56 | intergenic (‑142/+25) | otsB/araH | trehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base C (3/30); new base A (23/0); total (26/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.21e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGT > NC_000913/1982388‑1982679 | ccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTAT‑GGGGGGGAAAAAACCTaaaat > 2:244662/1‑148 (MQ=255) gTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAAttgttg < 2:929961/151‑1 (MQ=255) ttCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTTT‑GGGGGGGAAAAAAGCCaaaatttt > 2:1122349/1‑144 (MQ=255) cctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGaa < 1:43402/151‑1 (MQ=255) ctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAg < 2:128647/148‑1 (MQ=255) cctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGaaa < 2:769177/149‑1 (MQ=255) cctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGaaa < 1:669528/149‑1 (MQ=255) ctctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGGGGGGAAAAAAGCCAAATTGTTGGAGaaaaaa > 2:438473/1‑148 (MQ=255) tctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTTT‑TGGGGGGAAAAAAACCAAAATGTTTTAGaaaaaaaaa > 1:96199/1‑147 (MQ=255) ttttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTCT‑TGGGGGGAAAAAAACCAAATTTTTGTCGAAAAAAACa > 2:71652/1‑148 (MQ=255) ttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTATGGGGGGGGAAAAAAACAAAATT‑TTGTAGaaaaaaaaaaaa > 2:731859/1‑143 (MQ=255) ttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTAT‑GGTGGGGGAAAAAGACAAAATGTTGTAGaaaaaaaaa > 1:1015613/1‑143 (MQ=255) tATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGAc < 2:377436/148‑1 (MQ=255) aTTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGt < 2:888790/150‑1 (MQ=255) tgtgTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTAT‑GGGGGGGAAAAAAACCAAAATTTTTCAGaaaaaaaaaaaacaaaagggg > 1:929359/1‑136 (MQ=255) gtTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGtt < 1:276578/150‑1 (MQ=255) tGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTc < 1:385962/150‑1 (MQ=255) gcgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGCTGGGTTT‑GGGGGGGAAAAAAGCCAAAATTTTGCAGAAAAAAACAACAGAAAggggg > 2:235548/1‑141 (MQ=255) cgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTTT‑GGGGGGGAAAAAAACCAAAATGTTGCCGAAAAAAAAATCAGACAGAGCGttt > 2:1198141/1‑147 (MQ=255) cAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttt < 1:490170/151‑1 (MQ=255) tttACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGCCAT‑GGTGGGGAAAAAAACTAAATTGTTTTAGAAAAAAAAATCAAACAGCGGGGTTCGCctttttttttt > 1:216124/1‑143 (MQ=255) ttACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTaa < 2:829151/150‑1 (MQ=255) ttACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTa < 1:738257/149‑1 (MQ=255) aCAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACgg < 2:877739/151‑1 (MQ=255) cAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGt < 1:858331/151‑1 (MQ=255) aGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTc < 2:821955/151‑1 (MQ=255) cATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGaa < 2:177162/151‑1 (MQ=255) aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGCTTT‑GGGGGGGAAAAAAGCCAAATTTTTTCAGAAAAAAAAAACAAACAAGGCGGTTTGCTTTTTTCTTGTAAAgg > 2:315331/1‑145 (MQ=255) aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGATGGGGTT‑GGTGGGGAAAAAAGCTCAAATGTTTTAGAAAAAAACAACAGACAAGGCGTTTTGCTTTTTTTTTTTAAAGGGCGAAGaac > 2:1026245/1‑148 (MQ=255) aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAc < 1:503063/151‑1 (MQ=255) taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTTT‑GGGGGGGGAAAAAACTAAAATTTTTTAAAAAAAAACAACAAACAGTGGGGGTTGCTTTTTTTTTGttaacg > 1:999488/1‑131 (MQ=255) aGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTAC‑GGTGGGGGAAAAAACCAAATTTTTGCAGCAAAAAAGATCAGAAAGTGCGTTTTTCTTTTTTGTTTTAACGGGCGGAGATTACCGccg > 2:750832/1‑149 (MQ=255) aGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCt < 1:482689/150‑1 (MQ=255) gACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAgg < 1:65431/150‑1 (MQ=255) cGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAtt < 2:529003/151‑1 (MQ=255) gggggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTTT‑GGTGGGGGAAAAAACCAAATTTTTTTAGAAAAAAAAAACAGACAGTGCGTTTTGCTTTTTTTTTTTTAAGGGCGAAAAAAACCGCCGCCAGCAGGagt > 2:249898/1‑148 (MQ=255) gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTAT‑GGGGGGGAAAAAAGATAAATTTTTGCAGAAAAAAAAAACAGAAAAGGCGGTTCGCTTTTTTTTTTTAACGGGCGAAGAACAACGCCGCCAGCAGGATTaaaac > 2:616966/1‑148 (MQ=255) gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGcc < 1:1184712/151‑1 (MQ=255) tCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTTT‑TGTGGGGAAAAAACCCAAAATGTTTCAGAAAAAAGAATAAAACAATGCGGTTTGCTTTTTTTTTGTTACGGTCGAAGATTACCGCTGCCCGCAGGAtttaag > 2:59631/1‑145 (MQ=255) cGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTCT‑TGTGGGGAAAAAAACTAAAAATTTTTAGAAAAAAACATCAGACAGTGGGGTTCGCTTTTTTTTTTTTAAGGGCGGAGAACCACGCTGCcaggaagaatca > 1:1121123/1‑137 (MQ=255) cGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 2:171198/151‑1 (MQ=255) gCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgcg < 1:285991/151‑1 (MQ=255) cATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGGGGGGAAAAAAACCAAAATTTTGGAGAAAAAAAAAACAAACAGTGGGGTTTGCTTTTTTTTTTTTACGGGCGAAGAACACCCCCGCCCGCAGGGTTAAAcc > 2:1013860/1‑146 (MQ=255) tttGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTAT‑GGGGGGGAAAAAAGATAAAATTTTGCAGAAAAAAAAATAAAACAGTGCGGTTTGCTTTTTTTTTTTTAAAGGCGAAGAACACCGCCGCCCGCAGGGATAAGACCCGAAcc > 2:122094/1‑151 (MQ=255) gCAACGAAGAACAGCCAACCGTTATGGCGGTTTGCTAT‑GGGGGGGAAAAAAGCTAAATTTTTTGAGAAAAAAAAATCAGACAATGCGGTTTGCTTTTTTTTTTTTAAGGGCGAAGATCACCCCCGCCACCAGGAtttaaac > 2:680866/1‑136 (MQ=255) gCAACGAAGAACAGCCAACCGTTATGGCGGTTGGTTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAAcc < 2:110916/148‑1 (MQ=255) gaagaaCAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGGAAAAAAAATAAATTTTTGCAGGAAAAAACATTACACAAGGGGGTTTGCCTTTTTCTTTTTAAGGTTCAACATCACCCCCGCCAGCCCGGTTAACTCCCGCACCACCTAc > 2:966101/1‑149 (MQ=255) gaaCAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgcgc < 2:694733/151‑1 (MQ=255) aCCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgaga < 1:1075070/151‑1 (MQ=255) ggTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGc > 1:39273/1‑150 (MQ=255) gTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAg < 2:658009/151‑1 (MQ=255) tGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAgg < 2:950366/150‑1 (MQ=255) ggCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGtt < 1:570958/151‑1 (MQ=255) aT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATg > 2:381089/1‑151 (MQ=255) ggtggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGc < 2:101161/151‑1 (MQ=255) ggAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGttt > 2:79920/1‑151 (MQ=255) cGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACgg > 2:392430/1‑150 (MQ=255) cGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGt < 2:901488/151‑1 (MQ=255) | CCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGT > NC_000913/1982388‑1982679 |
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: ‑ |