Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 547,831:1 +G 100% pseudogene (251/261 nt)
pseudogene (251/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
*AP009048547,8311.G100.0% 134.8 / NA 39pseudogene (251/261 nt)
pseudogene (251/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 . (0/0);  new base G (19/20);  total (19/20)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC  >  AP009048/547761‑547896
                                                                       |                                                                 
gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:3641502/71‑1 (MQ=255)
gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:3869189/71‑1 (MQ=255)
           cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                         <  1:1030143/71‑1 (MQ=255)
            gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                        >  1:906758/1‑71 (MQ=255)
                  cTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                                  >  1:2773098/1‑71 (MQ=255)
                   tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                                >  1:391748/1‑55 (MQ=255)
                     ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                                  <  1:2245117/68‑1 (MQ=255)
                       cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCt                                              <  1:679014/70‑1 (MQ=255)
                        gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                             <  1:989723/70‑1 (MQ=255)
                         aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                           <  1:4036552/71‑1 (MQ=255)
                                ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                         <  1:2223687/66‑1 (MQ=255)
                                 gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                         >  1:4338434/1‑65 (MQ=255)
                                   gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  >  1:4371828/1‑70 (MQ=255)
                                   gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTACATTCCTTGTCACCATGTGa                                 <  1:2059092/71‑1 (MQ=255)
                                     aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaac                               >  1:538644/1‑70 (MQ=255)
                                      aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                              >  1:4461634/1‑71 (MQ=255)
                                      aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                              >  1:4547390/1‑71 (MQ=255)
                                           aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:4461392/62‑1 (MQ=255)
                                           aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                          <  1:118926/70‑1 (MQ=255)
                                             gcgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCgg                       >  1:4300686/1‑71 (MQ=255)
                                              cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                      >  1:396642/1‑71 (MQ=255)
                                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                    <  1:3039295/71‑1 (MQ=255)
                                                 aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                   <  1:2361891/71‑1 (MQ=255)
                                                  ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                  <  1:736124/71‑1 (MQ=255)
                                                  ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                  >  1:1310618/1‑71 (MQ=255)
                                                   gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:1632222/54‑1 (MQ=39)
                                                   gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                                    >  1:4352982/1‑52 (MQ=39)
                                                    cgcTCGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg                <  1:1885665/71‑1 (MQ=255)
                                                     gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                                    >  1:2882025/1‑50 (MQ=39)
                                                     gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg               <  1:3413505/71‑1 (MQ=255)
                                                       tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc             >  1:2306816/1‑71 (MQ=255)
                                                        ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                    <  1:1058583/63‑1 (MQ=255)
                                                         gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                      >  1:3214029/1‑44 (MQ=38)
                                                         gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg           >  1:25638/1‑71 (MQ=255)
                                                            aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt        <  1:996340/71‑1 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                  >  1:2903920/1‑60 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg              <  1:897351/64‑1 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                      >  1:3817926/1‑38 (MQ=38)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc     <  1:2516778/71‑1 (MQ=255)
                                                                tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg           >  1:2245726/1‑64 (MQ=255)
                                                                  gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                   >  1:538904/1‑54 (MQ=255)
                                                                  gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc  <  1:3847598/71‑1 (MQ=255)
                                                                       |                                                                 
GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC  >  AP009048/547761‑547896

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: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: