Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 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
*AP009048547,6940AG100.0% 170.0 / NA 50pseudogene (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 (29/21);  total (29/21)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG  >  AP009048/547628‑547761
                                                                  |                                                                   
tcgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGc                                                                  >  1:1273018/1‑70 (MQ=255)
 cgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCa                                                                 >  1:4341758/1‑70 (MQ=255)
              cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                   >  1:3465881/1‑71 (MQ=255)
                 acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAAATAAACGAGGGCAAAACACTGCTTCACGcc                                                >  1:3648563/1‑71 (MQ=255)
                  ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                               <  1:2683914/71‑1 (MQ=255)
                   tggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                           >  1:2051942/1‑58 (MQ=255)
                    ggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCggg                                             >  1:4319677/1‑71 (MQ=255)
                     gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTTCTTCACGCCggg                                             <  1:2897804/70‑1 (MQ=255)
                     gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgc                                            >  1:3959355/1‑71 (MQ=255)
                        ggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgc                                         >  1:2271620/1‑71 (MQ=255)
                         gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgc                                         >  1:3023/1‑70 (MQ=255)
                           tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                      >  1:2241334/1‑71 (MQ=255)
                               cAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgc                                   >  1:1486595/1‑70 (MQ=255)
                                aaCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcgc                                 <  1:2346551/71‑1 (MQ=255)
                                 aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCggg                                             <  1:4034763/58‑1 (MQ=255)
                                   cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                           >  1:912239/1‑42 (MQ=255)
                                    ggCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                <  1:2873691/36‑1 (MQ=255)
                                     gCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                        <  1:2069274/59‑1 (MQ=255)
                                        tcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                                <  1:11334/64‑1 (MQ=255)
                                        tcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                           >  1:4240982/1‑69 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc                                                <  1:5960/47‑1 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc                                                <  1:2871772/47‑1 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                           >  1:1348914/1‑68 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                         <  1:99707/70‑1 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                         <  1:2032351/70‑1 (MQ=255)
                                          ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGGCGGgccgcc                                        <  1:4265322/54‑1 (MQ=255)
                                           tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                   <  1:1469294/42‑1 (MQ=255)
                                             aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTGACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:258823/1‑71 (MQ=255)
                                             aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                     <  1:3329798/70‑1 (MQ=255)
                                             aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:4353189/1‑71 (MQ=255)
                                             aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    >  1:4248612/1‑71 (MQ=255)
                                               ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATg                                    <  1:788831/53‑1 (MQ=255)
                                               ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAcc                  >  1:3186532/1‑71 (MQ=255)
                                                tATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa                                                    >  1:4526579/1‑36 (MQ=255)
                                                  tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                                >  1:2167155/1‑54 (MQ=255)
                                                   cAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAc              >  1:753889/1‑71 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc             >  1:1641172/1‑71 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAc              <  1:2947354/70‑1 (MQ=255)
                                                      cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa           >  1:1170161/1‑71 (MQ=255)
                                                       gAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                 >  1:202727/1‑64 (MQ=255)
                                                       gAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAccc            <  1:3023862/69‑1 (MQ=255)
                                                         aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa        <  1:382973/71‑1 (MQ=255)
                                                            aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg     >  1:1807189/1‑71 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                      <  1:270088/37‑1 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGa               >  1:464778/1‑60 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg     <  1:291319/70‑1 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg     <  1:1157943/70‑1 (MQ=255)
                                                               cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaa       >  1:424979/1‑66 (MQ=255)
                                                               cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgggg  >  1:3247826/1‑71 (MQ=255)
                                                               cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgggg  >  1:2203911/1‑71 (MQ=255)
                                                                  |                                                                   
TCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG  >  AP009048/547628‑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: