Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 547,694 A→G 100% pseudogene (114/261 nt)
pseudogene (114/1259 nt)
ylbE →
ylbE →
ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment
ECK0512:JW0508+JW0507:b4572; hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb547,6940AG100.0% 105.7 / NA 32pseudogene (114/261 nt)
pseudogene (114/1259 nt)
ylbE
ylbE
ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment
ECK0512:JW0508+JW0507:b4572; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/10);  total (22/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG  >  W3110S.gb/547627‑547761
                                                                   |                                                                   
ttCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAggg                                                                   >  1:2596453/1‑70 (MQ=255)
               cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca                                                            >  1:563611/1‑62 (MQ=255)
               cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca                                                            >  1:1425701/1‑62 (MQ=255)
                  acTGGCTGGATGCGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                <  1:300756/55‑1 (MQ=255)
                   ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                <  1:2439146/54‑1 (MQ=255)
                                    cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                              >  1:2312496/1‑71 (MQ=255)
                                          cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACCCCGGGCCGCCAATGCGCTGGCAg                           >  1:485269/1‑68 (MQ=255)
                                           ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                           >  1:46562/1‑67 (MQ=255)
                                             cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgcc                                           >  1:1822527/1‑49 (MQ=255)
                                             cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgcc                                           >  1:2358025/1‑49 (MQ=255)
                                             cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgcc                                           >  1:24855/1‑49 (MQ=255)
                                                 tATCAGCGAACTAAACGAGGGCAAACCACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                 >  1:552065/1‑71 (MQ=255)
                                                 tATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                 >  1:2647691/1‑71 (MQ=255)
                                                 tATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                 >  1:1366637/1‑71 (MQ=255)
                                                 tATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                 >  1:2064701/1‑71 (MQ=255)
                                                 tATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                 >  1:1098865/1‑71 (MQ=255)
                                                  aTCAGCGAACTAAACGAGAGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:1121441/1‑67 (MQ=255)
                                                   tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcg                                  <  1:2196468/52‑1 (MQ=255)
                                                   tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcg                                  <  1:1226785/52‑1 (MQ=255)
                                                     aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc             >  1:2575271/1‑71 (MQ=255)
                                                     aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc             >  1:1534581/1‑71 (MQ=255)
                                                     aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc             >  1:1218662/1‑71 (MQ=255)
                                                          aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa        <  1:1559684/71‑1 (MQ=255)
                                                            tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      <  1:2155879/55‑1 (MQ=255)
                                                            tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      <  1:2023395/55‑1 (MQ=255)
                                                            tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      <  1:1929578/55‑1 (MQ=255)
                                                            tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      <  1:1886642/55‑1 (MQ=255)
                                                            tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      <  1:777163/55‑1 (MQ=255)
                                                             aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCa                                       >  1:774988/1‑37 (MQ=255)
                                                             aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGa               >  1:1998034/1‑61 (MQ=255)
                                                              aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAc                   >  1:1048323/1‑56 (MQ=255)
                                                              aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAc                   >  1:2581907/1‑56 (MQ=255)
                                                                cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgggg  >  1:916909/1‑71 (MQ=255)
                                                                   |                                                                   
TTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG  >  W3110S.gb/547627‑547761

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: