Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,972 (T)7→8 intergenic (‑173/‑150) yifB ← / → ilvL predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,911,9651.T98.3% 489.1 / 0.6 118intergenic (‑166/‑157)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (39/77);  total (40/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 125 total aligned reads displayed.

GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA  >  REL606/3911891‑3912034
                                                                           |                                                                     
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACGTTCATACTGAAATTGAAt                                                                      <  2:4833658/76‑1 (MQ=255)
 aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                     >  1:4679948/1‑76 (MQ=255)
   tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                   >  1:4547348/1‑76 (MQ=255)
        cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttc                                                              >  2:3490008/1‑74 (MQ=255)
        cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttc                                                              <  1:5120635/76‑3 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                           <  1:3710897/76‑1 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                           <  2:4205440/76‑1 (MQ=255)
              aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactcac                                                         >  2:2084070/1‑76 (MQ=255)
              aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactcac                                                         <  1:128520/76‑1 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTcactcact                                                        <  1:1983025/76‑1 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTcactcact                                                        <  2:1557258/76‑1 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTcactcact                                                        <  1:953368/76‑1 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTcactcact                                                        <  1:2424089/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                       <  2:3336199/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                       <  2:4072239/76‑1 (MQ=255)
                 aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtat                                                      <  1:4876661/76‑1 (MQ=255)
                 aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtat                                                      <  2:520544/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatt                                                     >  1:29173/1‑76 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatt                                                     <  2:893861/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatt                                                     <  2:1580096/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatt                                                     >  1:1677109/1‑76 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                    <  1:2596028/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                    <  2:1492459/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                    <  2:4946880/76‑1 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttt                                                   <  2:1175736/76‑1 (MQ=255)
                      gCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttat                                                 <  1:2038331/76‑1 (MQ=255)
                       cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttatt                                                <  1:1759561/76‑1 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttattt                                               <  1:4602528/76‑1 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttattt                                               <  2:182944/76‑1 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttattt                                               >  1:4686343/1‑76 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttattt                                               <  2:3546379/76‑1 (MQ=255)
                         gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATGCTGAAATTGAATTTTTTTTCACTCACTattttatttt                                              <  1:2657724/76‑1 (MQ=255)
                         gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTattttatttt                                              <  1:3049554/76‑1 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTAttttt                                             <  2:2449106/76‑1 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTAttttt                                             >  2:1515499/1‑76 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTAttttt                                             >  2:885303/1‑76 (MQ=255)
                            aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaa                                           <  1:4724171/76‑1 (MQ=255)
                            aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaa                                           >  2:3651299/1‑76 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaa                                          <  1:2067676/76‑1 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaa                                          <  2:4597542/76‑1 (MQ=255)
                              aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaa                                         <  1:2232205/76‑1 (MQ=255)
                               ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaaa                                        <  1:384491/76‑1 (MQ=255)
                               ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaaa                                        >  1:3910933/1‑76 (MQ=255)
                                gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaaaa                                       <  2:3507217/76‑1 (MQ=255)
                                gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaaaa                                       >  2:2423770/1‑76 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                       <  2:4167774/59‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                       >  1:4167774/1‑59 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaaca                                     <  1:3931481/76‑1 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                    <  1:2334769/76‑1 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                    >  1:2794292/1‑76 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                    >  1:685278/1‑76 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                    <  1:782050/76‑1 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                    >  2:130131/1‑76 (MQ=255)
                                    aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaac                                   <  1:1833977/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                  <  2:3815602/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                  <  1:1993417/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                  <  1:1182270/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                  <  2:4185216/76‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                >  1:4755927/1‑76 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                >  2:1186567/1‑76 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                >  2:2470273/1‑76 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                >  1:700114/1‑76 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                               <  2:3157875/76‑1 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                               <  2:1764333/76‑1 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                               >  1:1988812/1‑76 (MQ=255)
                                          gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                             >  1:3548001/1‑76 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  1:1874190/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  2:3132188/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  2:3834065/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  1:3973356/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  1:2794657/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                          <  1:352991/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                          <  1:2720036/76‑1 (MQ=255)
                                              aaGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATAtt                         <  1:1126864/76‑1 (MQ=255)
                                               aGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        >  1:1194877/1‑76 (MQ=255)
                                               aGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        <  1:3259885/76‑1 (MQ=255)
                                                 ggAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                      >  2:1712726/1‑76 (MQ=255)
                                                 ggAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                      >  2:1212984/1‑76 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                   >  2:1403458/1‑76 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                   <  2:1227225/76‑1 (MQ=255)
                                                     caacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                  <  1:303883/76‑1 (MQ=255)
                                                       acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta               <  1:3299186/76‑1 (MQ=255)
                                                        caTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta               <  2:1043220/76‑1 (MQ=255)
                                                         aTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa              <  2:3795714/76‑1 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc            >  1:2170843/1‑76 (MQ=255)
                                                             aTACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc          >  1:3304542/1‑76 (MQ=255)
                                                              tACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca         >  2:3141903/1‑76 (MQ=255)
                                                              tACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca         <  1:3427195/76‑1 (MQ=255)
                                                               aCTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        >  1:372568/1‑76 (MQ=255)
                                                               aCTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        >  2:4950272/1‑76 (MQ=255)
                                                               aCTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        >  1:415543/1‑76 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat     <  1:3199738/76‑1 (MQ=255)
                                                                   aaaTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa    >  1:3495052/1‑76 (MQ=255)
                                                                    aaTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa   <  2:674398/76‑1 (MQ=255)
                                                                    aaTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa   <  2:2959792/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  <  2:2344672/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  <  1:4071796/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  >  1:323721/1‑76 (MQ=255)
                                                                     aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  <  1:690130/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  <  1:2822241/76‑1 (MQ=255)
                                                                           |                                                                     
GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA  >  REL606/3911891‑3912034

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