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,6940AG97.9% 156.2 / ‑4.1 47pseudogene (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 (1/0);  new base G (21/25);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGG  >  AP009048/547630‑547759
                                                                |                                                                 
gtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                              <  1:5295621/70‑1 (MQ=255)
 tcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                            <  1:2558096/71‑1 (MQ=255)
 tcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                            <  1:21143/71‑1 (MQ=255)
   gcgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                                           <  1:2566496/70‑1 (MQ=255)
   gcgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                                           <  1:2518529/70‑1 (MQ=255)
   gcgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                                           <  1:1965163/70‑1 (MQ=255)
       tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGc                                                      >  1:650211/1‑71 (MQ=255)
          tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTc                                                   <  1:1141049/71‑1 (MQ=255)
                ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                             <  1:3177983/71‑1 (MQ=255)
                   gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgc                                          >  1:3310127/1‑71 (MQ=255)
                      ggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:3973091/1‑53 (MQ=255)
                       gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                      <  1:5193464/71‑1 (MQ=255)
                        aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                      >  1:2782924/1‑70 (MQ=255)
                             cAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgc                                 <  1:5061814/70‑1 (MQ=255)
                               aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                              <  1:852879/39‑1 (MQ=255)
                                 cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                            >  1:1720299/1‑71 (MQ=255)
                                   gCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:4867370/1‑40 (MQ=255)
                                     ttcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                         >  1:1891320/1‑70 (MQ=255)
                                      tcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                              >  1:1751406/1‑64 (MQ=255)
                                       cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCt                                                     <  1:2224935/40‑1 (MQ=255)
                                        ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGc                                                      <  1:4368080/38‑1 (MQ=255)
                                          cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag                      <  1:4511973/68‑1 (MQ=255)
                                          cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                    <  1:1310209/70‑1 (MQ=255)
                                          cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                    <  1:4562903/70‑1 (MQ=255)
                                          cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                    <  1:21369/70‑1 (MQ=255)
                                          cACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAgg                        >  1:1011501/1‑66 (MQ=255)
                                            cTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgc                                          <  1:747492/46‑1 (MQ=255)
                                                tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGa             >  1:3630765/1‑71 (MQ=255)
                                                 cAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGag            >  1:1367672/1‑70 (MQ=255)
                                                  aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                         <  1:4020426/57‑1 (MQ=255)
                                                   gCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                      <  1:4199646/43‑1 (MQ=255)
                                                    cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgaga                     >  1:3802117/1‑59 (MQ=255)
                                                    cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa         >  1:3616885/1‑71 (MQ=255)
                                                    cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa         >  1:3283599/1‑71 (MQ=255)
                                                      aaCTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                      <  1:240727/40‑1 (MQ=255)
                                                      aaCTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                      <  1:1876800/40‑1 (MQ=255)
                                                       aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgc                                 >  1:4794668/1‑44 (MQ=255)
                                                       aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATg       >  1:5351643/1‑70 (MQ=255)
                                                       aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa      <  1:4236773/71‑1 (MQ=255)
                                                           aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg  >  1:3948795/1‑71 (MQ=255)
                                                           aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg  <  1:3621835/71‑1 (MQ=255)
                                                           aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg  >  1:128883/1‑71 (MQ=255)
                                                           aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg  <  1:931419/71‑1 (MQ=255)
                                                           aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg   >  1:615834/1‑70 (MQ=255)
                                                            aCGAGGGCAAACCACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGCCCCATGa      >  1:2443469/1‑66 (MQ=255)
                                                            aCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa      >  1:4871109/1‑66 (MQ=255)
                                                             cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGctct                              >  1:2424362/1‑41 (MQ=37)
                                                                |                                                                 
GTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGG  >  AP009048/547630‑547759

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: