Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 1,530,173 | 0 | A | C | 33.3% | 20.1 / 15.7 | 18 | *399S (TAA→TCA) noncoding (62/244 nt) | ftsZ PP_mr19 | Cell division protein FtsZ Xoo8 ncRNA |
Reads supporting (aligned to +/- strand): ref base A (4/8); new base C (3/3); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCC > AE015451/1529963‑1530375 | agaCCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGt < 1:226124/235‑1 (MQ=255) aGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTAt > 1:167806/1‑237 (MQ=255) aCTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCTGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATAAGGGGTATAAAGGt > 2:403829/1‑226 (MQ=255) aaTCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTa > 1:238649/1‑173 (MQ=255) aaTCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTa < 2:238648/173‑1 (MQ=255) tGCGCCGGTGCGCCAGGGGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCGCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCa < 1:163398/234‑1 (MQ=255) cgcCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGt < 1:139697/184‑1 (MQ=255) ccTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTCCCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATACAGGTGATTGGTATTCATCACAGGTTGGTGCTGtntt > 2:53393/1‑181 (MQ=255) cTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTTCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGc < 2:167805/213‑1 (MQ=255) ggAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTatcatc < 1:93577/185‑1 (MQ=255) cATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCACTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGc < 1:206897/236‑1 (MQ=255) ccAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt > 2:65425/1‑110 (MQ=255) ccAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt < 1:65426/110‑1 (MQ=255) ccGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATAAAAGGTTg < 1:403831/121‑1 (MQ=255) cGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGt < 1:48611/62‑1 (MQ=255) cGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGt > 2:48611/1‑62 (MQ=255) cGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGc > 1:325271/1‑237 (MQ=255) atatCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGcc > 2:215194/1‑236 (MQ=255) cATTCCTGCGTCGTCAGTCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGc < 1:365083/226‑1 (MQ=255) cgtcgtCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCa < 1:24409/186‑1 (MQ=255) | ATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCC > AE015451/1529963‑1530375 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |