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.G98.0% 145.7 / ‑5.5 45pseudogene (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/1);  major base G (21/22);  minor base C (0/1);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

CGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACT  >  AP009048/547762‑547894
                                                                      |                                                                
cgtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAg                                                                 <  1:4555226/71‑1 (MQ=255)
     gtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggga                                                           <  1:3100342/71‑3 (MQ=255)
          cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                       <  1:2399479/71‑1 (MQ=255)
             aTGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                    <  1:4445299/71‑1 (MQ=255)
              tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                      >  1:2618850/1‑68 (MQ=255)
               gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                        >  1:5327539/1‑65 (MQ=255)
               gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                   >  1:3586173/1‑70 (MQ=255)
               gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                   >  1:1859535/1‑70 (MQ=255)
                tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                                                  <  1:3364501/70‑1 (MQ=255)
                tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                                                  <  1:2779176/70‑1 (MQ=255)
                  tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggga                                                           <  1:3752176/58‑3 (MQ=255)
                     tCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTcc                                             <  1:291322/70‑1 (MQ=255)
                       gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATAATGGAGCAGGGGGAAGTg                                                        >  1:2983135/1‑57 (MQ=38)
                             tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                    >  1:4161053/1‑71 (MQ=255)
                              gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                    <  1:403424/70‑1 (MQ=255)
                              gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                   >  1:1248224/1‑71 (MQ=255)
                                   aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                       <  1:1629986/46‑1 (MQ=38)
                                       aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          >  1:4826154/1‑71 (MQ=255)
                                        tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          <  1:4210573/70‑1 (MQ=255)
                                        tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                         <  1:3531652/71‑1 (MQ=255)
                                           aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg                      >  1:4618993/1‑71 (MQ=255)
                                           aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGGGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                       <  1:652005/70‑1 (MQ=255)
                                             cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                    <  1:4379555/71‑1 (MQ=255)
                                              gcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc                   >  1:4355594/1‑71 (MQ=255)
                                              gcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc                   >  1:5182206/1‑71 (MQ=255)
                                              gcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc                   >  1:4970448/1‑71 (MQ=255)
                                               cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTcc                                             <  1:3281504/44‑1 (MQ=38)
                                               cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                       <  1:609871/50‑1 (MQ=39)
                                               cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                       <  1:4576482/50‑1 (MQ=39)
                                               cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                  <  1:4248552/71‑1 (MQ=255)
                                                aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                 <  1:598339/71‑1 (MQ=255)
                                                   cgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg              <  1:218933/71‑1 (MQ=255)
                                                      tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc           >  1:1231384/1‑71 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                  <  1:1399929/62‑1 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         >  1:941793/1‑71 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         >  1:1040753/1‑71 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg          <  1:3594506/70‑1 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg          <  1:1416989/70‑1 (MQ=255)
                                                         cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         >  1:5007988/1‑70 (MQ=255)
                                                          aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                 >  1:3745350/1‑45 (MQ=255)
                                                           aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt      <  1:1129969/71‑1 (MQ=255)
                                                            tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAtt     <  1:3262338/71‑1 (MQ=255)
                                                            tACTGGAGCACGGGGTAGTGAACTTCATTCCTTGTCACCATGTGAAt                             <  1:1165760/47‑1 (MQ=255)
                                                              cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc   >  1:4765799/1‑71 (MQ=255)
                                                              cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc   >  1:5104381/1‑71 (MQ=255)
                                                               tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt  >  1:1608872/1‑71 (MQ=255)
                                                                ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt  >  1:4673566/1‑70 (MQ=255)
                                                                ggAGCAGGGTGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt  <  1:4246625/70‑1 (MQ=255)
                                                                      |                                                                
CGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACT  >  AP009048/547762‑547894

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: