| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | W3110S.gb | 907,713 | G→A | 100% | T331I (ACC→ATC) | ybjT ← | conserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 907,713 | 0 | G | A | 75.0% | 32.1 / 7.7 | 20 | T331I (ACC→ATC) | ybjT | conserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain | 
| Reads supporting (aligned to +/- strand): ref base G (5/0); new base A (8/7); total (13/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.14e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| ATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAA  >  W3110S.gb/907662‑907764                                                   |                                                   aTTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATa                                                   <  1:1513329/54‑1 (MQ=255)        cACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAAc                          >  1:969825/1‑71 (MQ=255)        cACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAAc                          >  1:2455512/1‑71 (MQ=255)        cACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAAc                          >  1:2396306/1‑71 (MQ=255)           tACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATa                                                   <  1:1473092/43‑1 (MQ=255)           tACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAAt                             <  1:1435956/65‑1 (MQ=255)           tACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAAt                             <  1:37523/65‑1 (MQ=255)            aCCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTa                      >  1:533854/1‑71 (MQ=255)            aCCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTa                      >  1:1250561/1‑71 (MQ=255)              cTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACt                    >  1:2454961/1‑71 (MQ=255)              cTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACt                    >  1:374457/1‑71 (MQ=255)              cTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACt                    >  1:1377931/1‑71 (MQ=255)                gCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTcc                  >  1:1651015/1‑71 (MQ=255)                gCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTcc                  >  1:1448944/1‑71 (MQ=255)                      aaGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACACCGCCTg                                        >  1:1383299/1‑43 (MQ=255)                             aaGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCggg     >  1:760631/1‑71 (MQ=255)                               ggCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCa   <  1:1758223/71‑1 (MQ=255)                               ggCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCa   <  1:814250/71‑1 (MQ=255)                                gCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCaa  <  1:480307/71‑1 (MQ=255)                                            tttAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGt               >  1:250316/1‑46 (MQ=255)                                                   |                                                   ATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAA  >  W3110S.gb/907662‑907764 | 
| Alignment Legend | 
|---|
| Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG | 
| Unaligned base: atcgMasked matching base:atcgAlignment gap:‑Deleted base:‑ |