Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 1,443,724 | 0 | T | . | 20.0% | 53.7 / 10.0 | 20 | intergenic (‑297/+56) | USA300HOU_1341/brnQ2 | IS1272 transposase/LIVCS family branched chain amino acid:cation symporter |
Reads supporting (aligned to +/- strand): ref base T (6/10); new base . (3/1); total (9/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.85e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGACCGAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTAATATCGTTGCAATAA > CP000730/1443627‑1443853 | aGACCGAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCtt < 1:193437/141‑1 (MQ=34) gACCGAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCttt < 2:71528/141‑1 (MQ=35) agaCGGCACACTATGAAAGCGAAGCCATTGAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCg < 2:71210/141‑1 (MQ=255) gCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTAtt < 1:202812/141‑1 (MQ=255) aGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATAtt > 2:265828/1‑141 (MQ=255) cATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGa < 2:364669/141‑1 (MQ=255) ttCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATc > 2:263876/1‑141 (MQ=255) cGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGtttt > 1:122305/1‑141 (MQ=255) gAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTa < 2:287240/141‑1 (MQ=255) aaGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTAc < 1:106375/141‑1 (MQ=255) ataAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGcc < 1:245548/141‑1 (MQ=255) aaTAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGAc < 1:123196/141‑1 (MQ=255) gaACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATa > 1:29848/1‑141 (MQ=255) cagTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTaa > 1:156723/1‑141 (MQ=255) tAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTAatat < 1:103016/141‑1 (MQ=255) aattttcaattagAAAA‑TATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCt < 2:192135/118‑1 (MQ=38) aattttcaattagAAAA‑TATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCt > 1:192135/13‑130 (MQ=38) tttcaattagAAAA‑TATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTAATATCGTTGCaat > 2:312537/10‑141 (MQ=255) tcaattagAAAA‑TATCTTACTGCTGTTTTTTTTTTTTTTTTTTTTACAACCAGATTCGTATTATTCTTTATTCCGGTAGATGGTTAATTGGATCTTGTTTTACAATGATTGAGAGTATATAGCATAATATCGTAGCaatta > 2:88178/8‑139 (MQ=255) cTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGTTAGCCTAATATCGTTGCaataa > 1:97165/1‑141 (MQ=255) | AGACCGAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTAATATCGTTGCAATAA > CP000730/1443627‑1443853 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |