| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,180,506 | 0 | A | C | 28.6% | 24.8 / 12.0 | 21 | E29A (GAG→GCG) | ygiC | ATP‑Grasp family ATPase |
| Reads supporting (aligned to +/- strand): ref base A (7/8); new base C (3/3); total (10/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAC > NC_000913/3180369‑3180641 | tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGGGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatggt < 6:138965/149‑6 (MQ=255) cGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCACCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgtgt > 8:146397/1‑135 (MQ=255) cGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCACCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgtgt < 7:146397/135‑1 (MQ=255) tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTACATGGtgtg < 8:147844/133‑1 (MQ=39) tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTACATGGtgtg > 7:147844/1‑133 (MQ=39) ctctTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGtt > 5:4645/1‑149 (MQ=255) ctctTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGtt > 5:2146/1‑149 (MQ=255) ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGCGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAc > 3:140930/1‑149 (MQ=255) ttCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGc < 7:60641/149‑1 (MQ=255) ccGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg < 6:2146/148‑1 (MQ=255) ccGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg < 6:4645/148‑1 (MQ=255) tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGa < 5:33360/127‑1 (MQ=255) tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGa > 6:33360/1‑127 (MQ=255) tACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGAAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGc > 5:32196/1‑149 (MQ=255) gcCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtg < 4:128592/149‑1 (MQ=255) gaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg > 1:25449/1‑148 (MQ=255) aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTAAAAGtggtgg > 4:125342/1‑148 (MQ=255) ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAg < 6:151946/94‑1 (MQ=255) tCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTAAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg < 3:125342/148‑1 (MQ=255) cacaCCA‑GTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga < 3:116844/148‑1 (MQ=255) accagTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc > 1:169163/5‑149 (MQ=255) | TGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAC > NC_000913/3180369‑3180641 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |