| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_007779 | 2,648,595 | 0 | A | C | 47.2% | 39.3 / 19.8 | 37 | G783G (GGT→GGG) | yfhM | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (16/3); new base C (0/17); total (16/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.35e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCG > NC_007779/2648497‑2648691 | ttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTAcc < 1:751340/101‑1 (MQ=255) aTCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACccg < 3:918141/101‑1 (MQ=255) cccGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCCCCACAGCGGCCCTTCACTGGACTCGACCATCGCATACCCTTTACCCGCCGt < 1:1278096/100‑1 (MQ=255) ggAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCCCCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTCCCCGCCGTTgg < 1:330859/101‑1 (MQ=255) gAATCGTCAGATCCAGCCCTTGAGCCCGAACTTCAATCTCTTGCCCCCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTCCCCGCCGTTGgc < 1:163144/101‑1 (MQ=255) tCAGACCCAGCCCTTGAGCCCAAACATCAATCTCTTGCCCCCACAGCGGCCCTTCACGGGACTCGACCATCGCATACCCTTTCCCCGCCGTTGGCgcggcg < 1:948087/101‑1 (MQ=255) agcccccccccTGGACCCCAAACATCATTCTCTTCCCCCCCCGGCGCCCCTTCCCTGGCCTCGCCCACCGCATACCCTTTCCCCGCCGTTGGCGCGGCGat < 1:14060/93‑1 (MQ=255) ccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACtt > 3:90916/1‑101 (MQ=255) ccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACtt > 1:1530778/1‑101 (MQ=255) gCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgt > 1:430171/1‑101 (MQ=255) ccGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtgt > 3:1411767/1‑101 (MQ=255) cGAACATCAATCTCTTGCCCCCACGGCGGCCCTTCCCTGGACTCGCCCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTc < 1:1548689/101‑1 (MQ=255) gAACATCAATCTCTTGCCCCCACGGCGGCCCTTCACTGGCCTCGACCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACGGAAGGgtgtgg < 4:12410/101‑3 (MQ=255) aCATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGcc > 3:677056/1‑101 (MQ=255) tCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAgg > 1:1285346/1‑101 (MQ=255) ctctTCCCCCCCCGGCGCCCCTTCACGGGCCTCGCCCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGAt < 1:962216/101‑1 (MQ=255) ctTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATaa > 1:852649/1‑101 (MQ=255) cccaccccaccGTCCCTACACGGGACTCCCCCACCGCTTACCCTTGTCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTg < 3:1271430/100‑1 (MQ=255) cccccccggcgcCCCTTCACGGGACTCGACCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTgg < 1:319341/93‑1 (MQ=255) ccaccaCAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTgg > 3:1582189/1‑101 (MQ=255) ttCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCa > 1:397676/1‑101 (MQ=255) ttCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCa > 1:182085/1‑101 (MQ=255) tccccTGGCCCCGCCCACCGAATACCCTTCCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCa < 1:947283/97‑1 (MQ=255) ccctgcccTCGCCCACCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGc < 1:711111/94‑1 (MQ=255) cACTGGACTCGACCACCGCATACCCTTTCCCCGCCGTTGGCGCGGTGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGc < 3:1245674/101‑1 (MQ=255) tcggcaCTCGCCCTTCCCATATCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCg < 1:117476/96‑1 (MQ=255) aCTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCg > 1:728442/1‑101 (MQ=255) aCTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCg > 1:1257610/1‑101 (MQ=255) aCTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCg > 1:1391754/1‑101 (MQ=255) gggccTCGCCCACCCCATACCCTTCCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTg < 1:293828/97‑1 (MQ=255) gccgcgcccACCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGac < 3:1200357/94‑1 (MQ=255) gcccATCGCATACCCTTCCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggt < 1:718879/99‑1 (MQ=255) cGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtc > 1:1262565/1‑101 (MQ=255) aTAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCa > 1:803024/1‑101 (MQ=255) ttttCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTgcg < 1:1109526/100‑1 (MQ=255) cTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTgcg < 1:694183/101‑1 (MQ=255) cTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTgcg > 3:434143/1‑101 (MQ=255) | TTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCG > NC_007779/2648497‑2648691 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |