| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,184,469 | 1 | . | A | 65.4% | 48.0 / 26.2 | 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/3184337‑3184604 | tCATATCACCTTCGTTTTCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCGCGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGAGGCATGG‑T‑a‑aacgaagcc > 1:300613/1‑136 (MQ=255) ttCGTTTTCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAACAGAGGCAAGTAGCGTCTGAGTCATGG‑TAA‑AA‑AA‑A‑ACCTCACGAAAAt > 1:61811/1‑143 (MQ=255) gTTTGCACGGTCTTCATCATCAAGCACCGTTGCACCGCGTCCATCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTAt < 1:326163/143‑1 (MQ=255) ttCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTACCGCACCTTCGACAGGCGTACCAAAAGAGGAAATTAGCGTCTGACCCATGG‑A‑A‑AA‑AACA‑ACAGCATAAAAATTATGGt > 2:130555/1‑143 (MQ=255) tcGGTCTACAGCATCAACCACCATTACACCGCGTCCGTCAACAAGAGCAGCAAGGGCATTTTCAACACGATCGAAAGGCGTAACAAAAGAAGAAAGTAGCATCTGATACATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTAc < 2:333169/142‑1 (MQ=255) tcatcaAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑T‑ACAA‑AA‑A‑ACCTCACTACAATTCTGGTTACCAGAATCCgg > 1:367832/1‑143 (MQ=255) cagcaAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCgggg < 1:79444/140‑1 (MQ=255) tcaAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATGCATGG‑T‑A‑AACAA‑A‑ACCTCACTAAAATTATGGTTACCAGAACTAGGGCa > 2:450342/1‑143 (MQ=255) gCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑T‑A‑AACAA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGCCt > 2:459652/1‑143 (MQ=255) aacATTACACCGCGTACTTCACGCAGCGCAGCCAGTGCATGTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCAAAG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTa < 2:383144/141‑1 (MQ=255) aCCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAACTATGGTTACCAGAAACAGGGCAGTCTTa > 2:455013/1‑143 (MQ=255) aCCAATACACCGCGTCCCTCACGCAGCGCAGCCAGTGCATTTTCAACACTTTCGAAAGGCGTACCAAAAGAGTAAAGTAGCGTCTGATTCATAG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTa < 1:376228/143‑1 (MQ=255) cATTACACCGCGTCCTTCACGCAGCGCAGCCAGGGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAgg < 2:212812/143‑1 (MQ=255) ccGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCCTAGGACTggcc > 1:243590/1‑142 (MQ=255) cgcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAACGTAGCGTCTGATTCATGG‑T‑A‑AAGAA‑A‑ACCGCACTAAACTTATGGTTACCAGTATCAGGGCAGTCTTAGGAGTggcgg > 1:439476/1‑143 (MQ=255) ttCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCCGACGCGCTGGTGTGTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCAt > 1:118483/1‑143 (MQ=255) acacGTTCTAAAGGCGTACCAAAAGAGGAAAGTAGCGTATGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTATCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGccc < 2:334273/143‑1 (MQ=255) aaaGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑TGA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGCCAGACTc > 1:489131/1‑143 (MQ=255) aaaaGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACtctctc < 1:492831/143‑1 (MQ=255) aaaaGAGGAAAGTAGCGACTGATACATGGAA‑A‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACtctctc < 2:282039/143‑1 (MQ=255) aaaGAGGAAAGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTcc > 2:441365/1‑143 (MQ=255) gAAAGTAGCGTCTGATTCATGG‑T‑A‑AA‑AA‑ACACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCg < 1:475688/143‑1 (MQ=255) aaGTAGCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGa < 1:487623/143‑1 (MQ=255) aggaGCGACTGAAACATGG‑TAA‑AA‑AA‑A‑ACCTCACTAACATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGAc < 2:329567/140‑1 (MQ=255) gtggagTCGGATGCGTGG‑T‑A‑AAGAA‑A‑ACCTCACGACAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGAct < 2:156100/138‑1 (MQ=255) ctAGCGTCTGATTCAGGG‑C‑A‑AA‑AA‑AGACCTCACTAAAATGATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGAct < 1:293424/142‑1 (MQ=255) gCGTCTGATTCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCGCATCCGGACTCTa > 2:3735/1‑143 (MQ=255) tCATGG‑TAA‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATACCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTAATATATCACATCAGGACACTAACCggcggc > 1:42291/1‑143 (MQ=255) | TCATATCACCTTCGTTTTCACGGTCTTCATCATCAAGCACCATTACACCGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGG‑T‑A‑AA‑AA‑A‑ACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGC > NC_000913/3184337‑3184604 |
| 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: ‑ |