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% 124.2 / NA 36pseudogene (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 (20/16);  total (20/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG  >  W3110S.gb/547628‑547761
                                                                  |                                                                   
tcgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGc                                                                  >  1:544034/1‑70 (MQ=255)
  gtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                <  1:970674/70‑1 (MQ=255)
   tcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                              <  1:882048/71‑1 (MQ=255)
   tcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                              <  1:1675215/71‑1 (MQ=255)
     gcgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                              >  1:1700390/1‑69 (MQ=255)
     gcgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                              >  1:1004111/1‑69 (MQ=255)
     gcgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                              >  1:1554540/1‑69 (MQ=255)
     gcgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                              >  1:1408236/1‑69 (MQ=255)
                         gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                        <  1:1322332/71‑1 (MQ=255)
                         gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                        <  1:1324489/71‑1 (MQ=255)
                         gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                        <  1:1533848/71‑1 (MQ=255)
                         gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgc                                         >  1:62652/1‑70 (MQ=255)
                                          ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                         <  1:1674191/69‑1 (MQ=255)
                                           tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      >  1:1408559/1‑71 (MQ=255)
                                           tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      >  1:1068731/1‑71 (MQ=255)
                                           tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      >  1:817435/1‑71 (MQ=255)
                                                 aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                 >  1:551985/1‑70 (MQ=255)
                                                   cAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      <  1:1512011/63‑1 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:1341698/1‑64 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:1079598/1‑64 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:120052/1‑64 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:1610887/1‑64 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:454837/1‑64 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:214348/1‑64 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc             <  1:1203262/71‑1 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc             <  1:878806/71‑1 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAc              <  1:123206/70‑1 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAc              <  1:807188/70‑1 (MQ=255)
                                                     gCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAc              <  1:697312/69‑1 (MQ=255)
                                                         aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa        <  1:238149/71‑1 (MQ=255)
                                                         aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa        <  1:1344278/71‑1 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                      >  1:239012/1‑37 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                      >  1:271454/1‑37 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                      >  1:821972/1‑37 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg    <  1:1385950/71‑1 (MQ=255)
                                                               cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgggg  >  1:1306744/1‑71 (MQ=255)
                                                                  |                                                                   
TCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG  >  W3110S.gb/547628‑547761

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