| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,184,467 | 1 | . | A | 65.4% | 47.9 / 26.3 | 27 | intergenic (‑3/‑371) | ribB/yqiC | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/BMFP family putative fusogenic protein |
| Reads supporting (aligned to +/- strand): ref base . (5/4); new base A (7/10); total (13/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.83e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCATATCACCTTCGTTTTCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑T‑A‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGC > NC_000913/3184335‑3184602 | tCATATCACCTTCGTTTTCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCGCGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGAGGCATGG‑T‑a‑aacgaagcc > 2:300613/1‑136 (MQ=255) ttCGTTTTCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAACAGAGGCAAGTAGCGTCTGAGTCATGG‑TAA‑AA‑AA‑A‑ACCTCACGAAAAt > 2:61811/1‑143 (MQ=255) gTTTGCACGGTCTTCATCATCAAGCACCGTTGCACCGCGTCCATCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTAt < 2:326163/143‑1 (MQ=255) ttCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTACCGCACCTTCGACAGGCGTACCAAAAGAGGAAATTAGCGTCTGACCCATGG‑A‑A‑AA‑AACA‑ACAGCATAAAAATTATGGt > 1:130555/1‑143 (MQ=255) tcGGTCTACAGCATCAACCACCATTACACCGCGTCCGTCAACAAGAGCAGCAAGGGCATTTTCAACACGATCGAAAGGCGTAACAAAAGAAGAAAGTAGCATCTGATACATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTAc < 1:333169/142‑1 (MQ=255) tcatcaAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑T‑ACAA‑AA‑A‑ACCTCACTACAATTCTGGTTACCAGAATCCgg > 2:367832/1‑143 (MQ=255) cagcaAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCgggg < 2:79444/140‑1 (MQ=255) tcaAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATGCATGG‑T‑A‑AACAA‑A‑ACCTCACTAAAATTATGGTTACCAGAACTAGGGCa > 1:450342/1‑143 (MQ=255) gCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑T‑A‑AACAA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGCCt > 1:459652/1‑143 (MQ=255) aacATTACACCGCGTACTTCACGCAGCGCAGCCAGTGCATGTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCAAAG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTa < 1:383144/141‑1 (MQ=255) aCCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAACTATGGTTACCAGAAACAGGGCAGTCTTa > 1:455013/1‑143 (MQ=255) aCCAATACACCGCGTCCCTCACGCAGCGCAGCCAGTGCATTTTCAACACTTTCGAAAGGCGTACCAAAAGAGTAAAGTAGCGTCTGATTCATAG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTa < 2:376228/143‑1 (MQ=255) cATTACACCGCGTCCTTCACGCAGCGCAGCCAGGGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAgg < 1:212812/143‑1 (MQ=255) ccGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCCTAGGACTggcc > 2:243590/1‑142 (MQ=255) cgcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAACGTAGCGTCTGATTCATGG‑T‑A‑AAGAA‑A‑ACCGCACTAAACTTATGGTTACCAGTATCAGGGCAGTCTTAGGAGTggcgg > 2:439476/1‑143 (MQ=255) ttCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCCGACGCGCTGGTGTGTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCAt > 2:118483/1‑143 (MQ=255) acacGTTCTAAAGGCGTACCAAAAGAGGAAAGTAGCGTATGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTATCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGccc < 1:334273/143‑1 (MQ=255) aaaGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TGA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGCCAGACTc > 2:489131/1‑143 (MQ=255) aaaaGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACtctctc < 2:492831/143‑1 (MQ=255) aaaaGAGGAAAGTAGCGACTGATACATGGAA‑A‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACtctctc < 1:282039/143‑1 (MQ=255) aaaGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTcc > 1:441365/1‑143 (MQ=255) gAAAGTAGCGTCTGATTCATGG‑T‑A‑AA‑AA‑ACACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCg < 2:475688/143‑1 (MQ=255) aaGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGa < 2:487623/143‑1 (MQ=255) aggaGCGACTGAAACATGG‑TAA‑AA‑AA‑A‑ACCTCACTAACATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGAc < 1:329567/140‑1 (MQ=255) gtggagTCGGATGCGTGG‑T‑A‑AAGAA‑A‑ACCTCACGACAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGAct < 1:156100/138‑1 (MQ=255) ctAGCGTCTGATTCAGGG‑C‑A‑AA‑AA‑AGACCTCACTAAAATGATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGAct < 2:293424/142‑1 (MQ=255) gCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCGCATCCGGACTCTa > 1:3735/1‑143 (MQ=255) tCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATACCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTAATATATCACATCAGGACACTAACCggcggc > 2:42291/1‑143 (MQ=255) | TCATATCACCTTCGTTTTCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑T‑A‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGC > NC_000913/3184335‑3184602 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |