Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,972 (T)7→6 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,9660T.98.8% 810.1 / 2.7 162intergenic (‑167/‑156)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (62/98);  total (63/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 174 total aligned reads displayed.

AATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGG  >  REL606/3911892‑3912040
                                                                          |                                                                          
aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                           >  2:1632836/1‑76 (MQ=255)
  tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                         >  2:2267962/1‑76 (MQ=255)
  tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                         >  1:3660140/1‑76 (MQ=255)
  tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                         >  1:5313547/1‑76 (MQ=255)
     tgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttc                                                                      >  2:3889633/1‑75 (MQ=255)
     tgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttc                                                                      <  1:2165969/76‑2 (MQ=255)
      gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttca                                                                     <  2:5836621/76‑3 (MQ=255)
       cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttcac                                                                    <  1:3084139/76‑4 (MQ=255)
         cGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactc                                                                 >  1:2630082/1‑76 (MQ=255)
          gCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactca                                                                <  2:340623/76‑1 (MQ=255)
          gCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactca                                                                >  1:1983412/1‑76 (MQ=255)
           cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactcac                                                               <  2:631249/76‑1 (MQ=255)
             aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACta                                                             <  1:5518670/76‑1 (MQ=255)
              aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtat                                                            <  1:3826099/76‑1 (MQ=255)
              aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtat                                                            <  1:2316641/76‑1 (MQ=255)
               aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatt                                                           <  1:3870660/76‑1 (MQ=255)
               aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatt                                                           <  1:3720152/76‑1 (MQ=255)
               aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatt                                                           <  1:5880961/76‑1 (MQ=255)
                aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattt                                                          <  2:3172894/76‑1 (MQ=255)
                aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattt                                                          <  1:4808964/76‑1 (MQ=255)
                 tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                         <  1:2130145/76‑1 (MQ=255)
                 tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                         <  2:292939/76‑1 (MQ=255)
                  gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttta                                                        <  1:3357646/76‑1 (MQ=255)
                  gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttta                                                        >  1:484082/1‑76 (MQ=255)
                   cAGCGGACAAAGGATGAACTACGAGGAAGGGAGCAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                       <  1:250350/76‑1 (MQ=255)
                   cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                       <  2:687506/76‑1 (MQ=255)
                   cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                       <  1:1553841/76‑1 (MQ=255)
                     gCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttattt                                                     <  1:6277586/76‑1 (MQ=255)
                      cGGACAAAGGATGAACTCCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTattttatttt                                                    >  2:1073101/1‑76 (MQ=255)
                      cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTattttatttt                                                    <  2:3040040/76‑1 (MQ=255)
                      cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTattttatttt                                                    >  1:2236520/1‑76 (MQ=255)
                        gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTa                                                  <  1:2991480/76‑1 (MQ=255)
                        gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTa                                                  >  1:4072280/1‑76 (MQ=255)
                         aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaa                                                 <  2:2531164/76‑1 (MQ=255)
                         aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaa                                                 <  2:1731760/76‑1 (MQ=255)
                         aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaa                                                 <  2:2055446/76‑1 (MQ=255)
                          cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                                <  1:5952588/76‑1 (MQ=255)
                           aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaa                                               >  2:4892593/1‑76 (MQ=255)
                            aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaa                                              >  2:2068786/1‑76 (MQ=255)
                            aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaa                                              <  2:412192/76‑1 (MQ=255)
                            aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaa                                              <  1:1893181/76‑1 (MQ=255)
                            aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaa                                              <  2:700838/76‑1 (MQ=255)
                             aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaaa                                             >  2:873201/1‑76 (MQ=255)
                             aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaaa                                             <  1:2888302/76‑1 (MQ=255)
                             aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaaa                                             >  1:341153/1‑76 (MQ=255)
                              ggATGAGCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                            <  2:6152396/76‑1 (MQ=255)
                              ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                            <  2:4973874/76‑1 (MQ=255)
                              ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                            <  1:2888756/76‑1 (MQ=255)
                               gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaaca                                           <  1:1342173/76‑1 (MQ=255)
                                aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaccaa                                          >  2:1443003/1‑76 (MQ=255)
                                aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                          >  1:3672654/1‑76 (MQ=255)
                                  gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                        <  2:842694/76‑1 (MQ=255)
                                   aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                       >  1:5606938/1‑76 (MQ=255)
                                    aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                      <  1:1809122/76‑1 (MQ=255)
                                      tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                    <  1:1059989/76‑1 (MQ=255)
                                      tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                    <  1:1737748/76‑1 (MQ=255)
                                      tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                    >  1:192553/1‑76 (MQ=255)
                                       aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                    >  1:5057256/1‑76 (MQ=255)
                                       aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                                   >  2:6158367/1‑76 (MQ=255)
                                       aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                                   >  2:5133639/1‑76 (MQ=255)
                                        cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                                  <  1:556135/76‑1 (MQ=255)
                                        cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                                  <  2:5170704/76‑1 (MQ=255)
                                          aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                                >  1:6139574/1‑76 (MQ=255)
                                          aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                                >  1:2390366/1‑76 (MQ=255)
                                          aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                                >  1:3417210/1‑76 (MQ=255)
                                            gAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                              <  1:4923966/76‑1 (MQ=255)
                                             aaGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                             <  1:1264644/76‑1 (MQ=255)
                                             aaGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                             <  1:2696928/76‑1 (MQ=255)
                                                  aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                        >  2:1920519/1‑76 (MQ=255)
                                                    caacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                      >  2:4741573/1‑76 (MQ=255)
                                                    caacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                      <  2:2400949/76‑1 (MQ=255)
                                                    caacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                      <  2:3937895/76‑1 (MQ=255)
                                                    caacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                      <  2:3256522/76‑1 (MQ=255)
                                                     aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                     >  2:3748191/1‑76 (MQ=255)
                                                     aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                     <  1:3350813/76‑1 (MQ=255)
                                                     aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                     <  1:877477/76‑1 (MQ=255)
                                                     aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                     >  2:4391670/1‑76 (MQ=255)
                                                     aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                     >  2:4572879/1‑76 (MQ=255)
                                                      acaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa                    <  2:4303302/76‑1 (MQ=255)
                                                       caTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa                   >  2:3499673/1‑76 (MQ=255)
                                                        aTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc                  <  2:2525620/76‑1 (MQ=255)
                                                        aTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc                  >  2:5664130/1‑76 (MQ=255)
                                                        aTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc                  >  2:1268750/1‑76 (MQ=255)
                                                         ttCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg                 >  1:935762/1‑76 (MQ=255)
                                                         ttCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg                 >  1:5613168/1‑76 (MQ=255)
                                                            aTACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat              <  1:5262217/76‑1 (MQ=255)
                                                               cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat           <  2:4671601/76‑1 (MQ=255)
                                                               cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat           <  1:2679102/76‑1 (MQ=255)
                                                               cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat           <  1:4563643/76‑1 (MQ=255)
                                                                tGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa          >  2:307277/1‑76 (MQ=255)
                                                                tGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa          <  2:6179559/76‑1 (MQ=255)
                                                                tGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa          >  1:933101/1‑76 (MQ=255)
                                                                 gAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa         >  1:4423387/1‑76 (MQ=255)
                                                                 gAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa         >  2:301004/1‑76 (MQ=255)
                                                                  aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa        <  1:345410/76‑1 (MQ=255)
                                                                  aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa        <  1:4396151/76‑1 (MQ=255)
                                                                  aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa        <  1:4730334/76‑1 (MQ=255)
                                                                   aaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa       <  2:4786098/76‑1 (MQ=255)
                                                                     ttGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt     >  1:1086724/1‑76 (MQ=255)
                                                                     ttGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt     >  2:4433421/1‑76 (MQ=255)
                                                                          |                                                                          
AATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGG  >  REL606/3911892‑3912040

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: