Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 52.9% | 27.2 / 51.7 | 34 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (14/2); new base C (0/18); total (14/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.62e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-08 |
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT > NC_000913/357730‑358022 | aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt < 2:233908/151‑1 (MQ=255) tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc > 1:547058/1‑151 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa < 1:582490/151‑1 (MQ=255) tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa > 1:816645/1‑150 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 1:368459/151‑1 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag > 2:615794/1‑150 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 2:679397/1‑151 (MQ=255) ccGCAGCCCTGTGTCGGCGGCTACCTTGATTCATTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAAATGCTTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc < 1:700426/147‑1 (MQ=255) cGCAGCCCTGTGTCAGCGGCTCCCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTATCCACCCCCATGTGAAAAAATGCTTTCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg < 2:97996/149‑1 (MQ=255) gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg > 1:155485/1‑149 (MQ=255) ggCTACCGTAATTCATTCCCGCCAACACCCGCGCATTTATCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTa < 2:816736/150‑1 (MQ=255) tcccGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCGGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAg < 1:252466/147‑1 (MQ=255) ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGA‑G‑a‑aaaaaaatttagaac > 1:684199/1‑137 (MQ=25) ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGG‑A‑AAAAACTTTAAAccc > 1:296535/1‑151 (MQ=25) ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑G‑A‑AAAAACATTAAAccc > 1:921209/1‑150 (MQ=25) cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑G‑AGAAAAGTTTTAAGc > 1:182016/1‑148 (MQ=25) cGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑A‑AAAAAATTTTAACCCCGCATGTTTTGTgg > 2:709007/1‑150 (MQ=25) acaccCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg < 2:710308/151‑1 (MQ=255) acccGCGCATTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 2:707408/149‑1 (MQ=255) gcttttATCCAACCCCATGTGAAAAAATGCCTCCGGGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct < 2:634723/145‑1 (MQ=255) aTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑A‑AAAAAATTTTAAACCGCCCTGTTGTTTGGGGAAACCTGCT‑GGGAACaa > 2:257698/1‑150 (MQ=21) ccACCGCCTTGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACa < 1:590049/147‑1 (MQ=255) ccAACGCCATGTCAAAAAATCCCTTCGCGGCGGCTGTCTGCGAGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAtgt < 1:248131/151‑1 (MQ=255) cGCCATGTGAAAAAATGCTTGCGCGGCGGTTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGa < 1:562627/150‑1 (MQ=255) ccttgtcaaaaaaTCCTTCCGCGGCGGTTGTCTGCCACCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGaaa < 2:622272/142‑1 (MQ=255) cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAA‑AAAAATTTTAAGCCC‑TCATGTTTTGTGGGGATAGCCGCT‑GGGGAAAATGTGGGAAgg > 2:887453/1‑151 (MQ=25) cATGTCCAAAAATTCTTTCGCGGCGGCTGTCCGCCAGCTTTTTTCCCGCCGCACCAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAg < 2:463042/150‑1 (MQ=255) tgaaaaaaTGCCTTCGGGGCGCCTGTCACCCAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGa < 1:708221/146‑1 (MQ=255) tgaaaaaaGGCCTTCGGGGGGGCTTTCTGCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGa < 2:719489/146‑1 (MQ=255) aaaaaTGCCTTCGCGGCGCCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGt < 2:860759/149‑1 (MQ=255) aaaaaTCCCTTCGCGGCGGCTGTCGGCAACCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTgc < 2:604989/151‑1 (MQ=255) cgggggTTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 1:338929/147‑1 (MQ=255) gCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGCGA‑G‑G‑AAAAAAATTAAAAACCTCCTGTTGTGTGGGCATCGCTGCT‑GGTGAAAATAAGGGAAGGGGAGTGCGGCGGAGGAGCCCCACAAc > 1:767619/1‑145 (MQ=255) cgatcttttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg < 1:876150/141‑1 (MQ=255) cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGG‑A‑AAAAATTTTAACACCACCTAGTT‑TGTGGGAATAGCCGCT‑GGGGAAAATGGGGGGAGGGGAGGGGGGGGGATCCACTCCAGCAGCGCGCTAAt > 1:204919/1‑149 (MQ=255) tATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc < 2:661488/144‑1 (MQ=255) ttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGACGG‑G‑A‑AAAAAATTTTAGCCCCCCATGTTTTGTGGGGAACGCTGCTGGGAAAAAAGGGGGAAAGGGAGGGGGGGGGAA‑CAACTCCAGAAAAGCGCCaaat > 1:701817/1‑144 (MQ=255) | ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑A‑AATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT > NC_000913/357730‑358022 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |