| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,386,862 | 0 | G | A | 64.7% | 3.1 / 28.3 | 34 | intergenic (+31/‑72) | yfbL/yfbM | putative M28A family peptidase/DUF1877 family protein |
| Reads supporting (aligned to +/- strand): ref base G (3/9); new base A (22/0); total (25/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.19e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-06 | |||||||||||
GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑G‑GAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAC > NC_000913/2386725‑2387006 | gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAA‑CGACAA‑GCACatt > 1:237401/1‑151 (MQ=255) aaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAA‑CGACAT‑GCACattttattt > 2:311790/1‑151 (MQ=255) caaCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAA‑AGACATGGCAC‑TTTTATTTTTCCt > 1:370570/3‑149 (MQ=255) aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATc < 2:629223/151‑1 (MQ=255) tGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAA‑AGACAA‑GCAAaattttttttttcttctcttt > 1:407025/1‑140 (MQ=255) cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAA‑AGAAAA‑ACCCattttattttttccccccttttt > 1:690478/1‑138 (MQ=255) gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAA‑AGACAT‑ACACATTTTATTTTTTCTCCTCTGtttttaa > 1:401263/1‑148 (MQ=255) tACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGc < 1:636677/151‑1 (MQ=255) aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAG‑CGACa < 1:291068/114‑1 (MQ=255) ccGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAA‑AGACAT‑GCACAATTTATTTTTCCTCATCTGTTTTTAGATGCGa > 2:599858/1‑151 (MQ=255) acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAAc < 1:398325/148‑1 (MQ=255) gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAA‑CGACAT‑TCAAAATTTTTTTTTCCTCATCTTTTTTTAGATGCGAAAAAAGaaa > 2:590185/1‑151 (MQ=255) aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAt < 2:476475/151‑1 (MQ=255) tCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAA‑GGAAAT‑GCACAATTTTTTTTTTCTCATCTGTTTTAAGAGGCGGAACAAGAAAATCAACgga‑g > 2:221536/1‑151 (MQ=255) tCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAA‑AGACAA‑GCACATTTTTTTTTTCCTCCTCTGTTTTTAAATGCGAAAAATGAAAATCAAAgga‑g > 2:632663/1‑151 (MQ=255) gggaGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGGAAAAA‑CGACAA‑AAAAAATTTATTTTTTCTCATCTGTTTTTAGAAAAGAAAAAATAAAATAAACGGA‑G‑GGAAAAAAAATtggg > 1:448293/4‑150 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAA‑AGAAAA‑GGACATTTTATTTTTTCTCATCTTTTTTTAGAAGGGAAAAAAGAAAATTAACGGA‑G‑GGaaaaaaaaaaggg > 1:66777/1‑146 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAA‑AGAAAA‑GCACATTTTTTTTTTTCTCATCTGTTTTTAGAAGAGAAAAAAGAAAATTAAAGGA‑G‑GGAAAAAAAATGGGGATGa > 1:119950/1‑151 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAA‑CGACAA‑GCACATTTTTTTTTTCCTCCTCTGTTTTTAAATGCGGAAAAAGAAAATCAACGGAGG‑GAAAAAAAAATGGGatgg > 2:482381/1‑148 (MQ=255) tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAA‑AGACAT‑GCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAAAGAAAATTCACGGG‑G‑GGAAAATAAATGGGgatgat > 2:581311/1‑151 (MQ=255) ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAA‑CGACAT‑GCACAATTTATTTTTCCTCATCTTTTTTTAGATGCGAAAAATGAAAATTAACCGA‑G‑GAAAAATAAATGGGAATGaat > 2:83933/1‑149 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAA‑AGACAA‑GCACATTTTATTTTTTCTCATCTGTTTTTAGAAGCGAAAAAAGAAATTTAACGGA‑G‑GAAAAATAAATTGGGATGATTGGCTa > 1:573832/1‑151 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAA‑CGACAT‑GCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAAAATGAAATTCAAAGGA‑G‑GAAAAATAAATGGGAATGATTGGCTa > 2:418092/1‑151 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAA‑CGACAT‑GCACATTTTATTTTTTCTCATCTGTTTTTAGATGGGAAAAATTAAATTTAAAGGA‑G‑GAAAAATAAATGGGGATGATTGGCtt > 1:598687/1‑150 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGAGCGGAAAAAGCCAACAA‑TAAACTTTAATTTTCCTTCATCT‑TTTTGAAATTCGAAAAAATAAATTATACCGC‑GAGGAAAATAAATggatgtgtttgttg > 1:419186/1‑139 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAA‑CGACAT‑GCACAATTTTTTTTTCCCCCTCTGTTTTTAAATGCGAAAAATGAAAATCAAAGGA‑G‑GGAAAAAAAATGGGAATGATTTGGTATTTTTTAGaaa > 2:383229/1‑151 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAA‑AGACAA‑GAAAATTTTATTTTTCCTCCTCTGTTTTTAGATGCGGAAAAAGAAAATTAACGGAGG‑GAAAAAAAAATGGG‑ATTATTTGCTATTTTTTAGAAAATGAATCAGa > 1:428739/1‑151 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAA‑AGACAA‑ACAAAATTTATTTTTTCTCCTCTGTTTTTAGAAGGGAAAAAAGGAAATTAAAGGA‑G‑GAAAAATAAAAGGGGATGATTTGGTATTTTTCAGaaaattaatcaaa > 1:655488/1‑141 (MQ=255) ctcGAGGCGGGTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑G‑GAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTg < 2:690929/151‑1 (MQ=255) gCGGGTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑G‑GAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc < 2:34163/149‑1 (MQ=255) gTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑G‑GAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGa > 2:304846/1‑151 (MQ=255) aGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑G‑GAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa > 2:88364/1‑151 (MQ=255) aaaaaG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑G‑GAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc < 1:181397/150‑1 (MQ=255) aaaaaG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑G‑GAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc < 1:646602/150‑1 (MQ=255) | GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAG‑CGACAT‑GCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑G‑GAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAC > NC_000913/2386725‑2387006 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |