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.T97.0% 670.2 / 9.6 167intergenic (‑166/‑157)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base T (46/116);  total (47/120)
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

Only 100 of 190 total aligned reads displayed.

GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAA  >  REL606/3911891‑3912035
                                                                           |                                                                       
gAATGATGCTCGCAAAAATGCAGCGGCCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                        <  2:6508141/76‑1 (MQ=255)
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                        <  1:6132185/76‑1 (MQ=255)
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                        <  1:3562760/76‑1 (MQ=255)
  attatgCTCGCAACAAGTCAGCGGACAAAGGTTTAACTACGAGGAAGGGAACAACATTCATACTGAaatttaattt                                                                      <  2:4676495/73‑1 (MQ=255)
  atgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                                                      <  1:7831440/76‑1 (MQ=255)
  atgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                                                      <  2:7114965/76‑1 (MQ=255)
   tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                     <  2:7239841/76‑1 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                  <  2:35415/76‑1 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                  >  2:3736576/1‑76 (MQ=255)
      tgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttt                                                                 <  2:5191962/76‑1 (MQ=255)
       gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttttt                                                                >  1:5036643/1‑75 (MQ=255)
       gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttttt                                                                >  1:4007966/1‑75 (MQ=255)
       gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttttt                                                                >  2:5639675/1‑75 (MQ=255)
        cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttc                                                               <  1:6563562/76‑3 (MQ=255)
         tCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttca                                                              <  2:6291156/76‑4 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                            <  1:2835533/76‑1 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                            >  2:3311660/1‑76 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                            <  2:2956290/76‑1 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactca                                                           <  2:5911496/76‑1 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactca                                                           <  1:7153287/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                        <  2:2999251/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                        <  2:1194214/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                        <  1:5457192/76‑1 (MQ=255)
                 aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtat                                                       <  2:5727791/76‑1 (MQ=255)
                 aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtat                                                       <  1:8003093/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatt                                                      <  1:5053518/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatt                                                      >  2:1791357/1‑76 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                     <  1:1640529/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                     <  1:4042694/76‑1 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttt                                                    <  2:5017044/76‑1 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttt                                                    <  1:6793724/76‑1 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttt                                                    <  1:115315/76‑1 (MQ=255)
                      gCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttat                                                  <  1:5954097/76‑1 (MQ=255)
                      gCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttat                                                  <  1:1314155/76‑1 (MQ=255)
                       cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttatt                                                 <  2:2829605/76‑1 (MQ=255)
                         gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTattttatttt                                               <  2:158691/76‑1 (MQ=255)
                         gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTattttatttt                                               <  1:981820/76‑1 (MQ=255)
                           aaaaaGATGAACTACGAAGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTa                                             <  1:3092175/75‑1 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttatt                                                 >  2:2608532/1‑70 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaa                                           >  2:2240421/1‑76 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaa                                           <  1:6899516/76‑1 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaa                                           <  2:7366003/76‑1 (MQ=255)
                              aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaa                                          >  2:6677235/1‑76 (MQ=255)
                                 aTGAACTACGAGGAAGGTAACAACATTCATACGGTAATTGCATTTTGTTTCTCACACTATTTTATTTTTAAAAaac                                       <  2:4505872/76‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaac                                       <  2:5543736/76‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaac                                       <  1:6219127/76‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaac                                       <  1:6748400/76‑1 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaaca                                      >  1:4287826/1‑76 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                     <  2:1936724/76‑1 (MQ=255)
                                    aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaac                                    <  2:4703996/76‑1 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTGTTTTATTTTTAAAAAACAACAAt                                 <  2:7828717/76‑1 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                 >  2:6762065/1‑76 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                >  2:6794037/1‑76 (MQ=255)
                                         cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                               <  1:3366197/76‑1 (MQ=255)
                                         cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                               <  1:5932194/76‑1 (MQ=255)
                                          gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                              <  2:2814478/76‑1 (MQ=255)
                                          gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                              <  2:2709791/76‑1 (MQ=255)
                                          gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                              <  2:8086655/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                             <  1:3694034/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                             <  1:2967974/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                             >  2:5765491/1‑76 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                             <  2:6160680/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                           <  2:6997746/76‑1 (MQ=255)
                                              aaGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATAtt                          <  1:5119281/76‑1 (MQ=255)
                                              aaGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATAtt                          <  1:1772615/76‑1 (MQ=255)
                                               aGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                         >  2:6912784/1‑76 (MQ=255)
                                               aGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                         >  1:3901148/1‑76 (MQ=255)
                                                gggAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                        >  2:3444763/1‑76 (MQ=255)
                                                gggAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                        <  2:4468656/76‑1 (MQ=255)
                                                 ggAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                       <  1:1932519/76‑1 (MQ=255)
                                                 ggAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                       <  2:6551037/76‑1 (MQ=255)
                                                  gAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAGCAATTTATATTGaaa                      <  1:765783/76‑1 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat                     <  1:7804680/76‑1 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                    <  1:6121020/76‑1 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAata                    <  2:5383138/76‑2 (MQ=255)
                                                     caacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                   >  1:5872263/1‑76 (MQ=255)
                                                     caacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                   >  1:6656638/1‑76 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                  <  1:2899370/76‑1 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                  <  2:6123785/76‑1 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                  >  1:2995684/1‑76 (MQ=255)
                                                       acaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                 <  1:6321212/76‑1 (MQ=255)
                                                       acaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                 >  1:6406447/1‑76 (MQ=255)
                                                          ttCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa              <  1:5538480/76‑1 (MQ=255)
                                                          ttCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa              <  2:1466279/76‑1 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc             >  2:3674683/1‑76 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc             >  1:1097775/1‑76 (MQ=255)
                                                            cATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg            <  1:5614404/76‑1 (MQ=255)
                                                             aTACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc           <  1:233618/76‑1 (MQ=255)
                                                                cTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc        >  1:115852/1‑76 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      >  1:2590722/1‑76 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      <  1:4653544/76‑1 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      <  2:4753609/76‑1 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      <  1:8025017/76‑1 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      <  1:7759491/76‑1 (MQ=255)
                                                                   aaaTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa     <  2:4228362/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa  >  2:107664/1‑76 (MQ=255)
                                                                      ttGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa  >  2:3024407/1‑76 (MQ=255)
                                                                      ttGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa  <  1:4440558/76‑1 (MQ=255)
                                                                      ttGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa  <  2:390123/76‑1 (MQ=255)
                                                                      ttGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa  <  1:6006931/76‑1 (MQ=255)
                                                                           |                                                                       
GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAA  >  REL606/3911891‑3912035

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