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.97.0% 628.2 / 12.1 135intergenic (‑167/‑156)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base . (38/93);  total (39/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 149 total aligned reads displayed.

GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAT  >  REL606/3911891‑3912037
                                                                           |                                                                       
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAAGAAGGGAACAACATTCATACTGAAATTGAAt                                                                         <  2:4161396/76‑1 (MQ=255)
 aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                        <  1:5919075/76‑1 (MQ=255)
  atgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                                                       <  1:1602314/76‑1 (MQ=255)
    gatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttt                                                                     >  1:4928639/1‑76 (MQ=255)
    gatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttt                                                                     <  1:5780707/76‑1 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                    <  1:5579631/76‑1 (MQ=255)
       gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttca                                                                  <  1:4695225/76‑3 (MQ=255)
          cGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactc                                                              >  2:3437792/1‑76 (MQ=255)
           gCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactca                                                             >  1:741291/1‑76 (MQ=255)
           gCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactca                                                             <  2:1905657/76‑1 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactcac                                                            >  2:4500527/1‑76 (MQ=255)
              aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACta                                                          <  1:5419512/76‑1 (MQ=255)
              aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACta                                                          <  1:5014138/76‑1 (MQ=255)
              aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACta                                                          <  1:662392/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatt                                                        <  2:5530771/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatt                                                        <  1:2137507/76‑1 (MQ=255)
                 aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattt                                                       <  1:1461145/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                      <  1:2457107/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                      <  2:1042325/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                      <  1:3127801/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                      >  1:1242883/1‑76 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                    <  2:1145879/76‑1 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                    >  1:5619400/1‑76 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                    <  1:4088573/76‑1 (MQ=255)
                      gCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttattt                                                  <  1:2320372/76‑1 (MQ=255)
                         gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTa                                               <  1:1737457/76‑1 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaa                                              <  1:514372/76‑1 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaa                                              <  2:5720966/76‑1 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                             >  1:4608228/1‑76 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                             >  2:5104344/1‑76 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                             >  1:4396300/1‑76 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                             <  2:3908997/76‑1 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                             <  1:5796575/76‑1 (MQ=255)
                              aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaaa                                          <  2:1519476/76‑1 (MQ=255)
                                gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaaca                                        <  2:4753954/76‑1 (MQ=255)
                                gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaaca                                        <  1:739268/76‑1 (MQ=255)
                                gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaaca                                        <  1:2478071/76‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                       <  2:2146284/76‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                       <  2:400921/76‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                       <  2:5874415/76‑1 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaac                                      <  1:5071624/76‑1 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaac                                      >  2:4217670/1‑76 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaac                                      <  1:1370266/76‑1 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                     >  1:5399795/1‑76 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                   <  2:1755649/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                   <  1:3807879/76‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  >  2:258932/1‑76 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  <  2:3100169/76‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  <  2:1158046/76‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  <  1:1875751/76‑1 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                 <  2:37545/76‑1 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                 <  1:4285361/76‑1 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                 >  2:3858790/1‑76 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                                <  2:521553/76‑1 (MQ=255)
                                         cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                               <  1:5310776/76‑1 (MQ=255)
                                          gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                              <  2:1917176/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                             <  1:3923588/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                           <  1:255963/76‑1 (MQ=255)
                                              aaGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                          >  1:2530857/1‑76 (MQ=255)
                                                gggAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaaa                        <  1:5407357/76‑1 (MQ=255)
                                                 ggAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat                       <  1:17706/76‑1 (MQ=255)
                                                  gAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                      <  1:325170/76‑1 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                     <  2:3532383/76‑1 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                     >  1:6129768/1‑76 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                     <  1:4403298/76‑1 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                    <  1:1191469/76‑1 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                    >  1:2630384/1‑76 (MQ=255)
                                                     caacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                   <  1:4638760/76‑1 (MQ=255)
                                                     caacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                   >  2:3660816/1‑76 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                  <  2:4547413/76‑1 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                  >  1:1049252/1‑76 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                  >  2:1497797/1‑76 (MQ=255)
                                                       acaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa                 <  1:5245762/76‑1 (MQ=255)
                                                        caTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa                <  2:742841/76‑1 (MQ=255)
                                                        caTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa                >  2:3121212/1‑76 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc             <  1:3773502/76‑1 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc             >  1:1040987/1‑76 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc             <  2:1036367/76‑1 (MQ=255)
                                                            cATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca            >  1:160985/1‑76 (MQ=255)
                                                            cATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca            >  1:6308981/1‑76 (MQ=255)
                                                            cATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca            <  2:1779154/76‑1 (MQ=255)
                                                             aTACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat           >  2:2590122/1‑76 (MQ=255)
                                                              tACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc          <  2:3337616/76‑1 (MQ=255)
                                                              tACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc          <  2:2367539/76‑1 (MQ=255)
                                                               aCTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca         <  1:4857106/76‑1 (MQ=255)
                                                               aCTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca         >  2:5642051/1‑76 (MQ=255)
                                                               aCTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca         <  1:1757354/76‑1 (MQ=255)
                                                               aCTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca         <  1:3218547/76‑1 (MQ=255)
                                                               aCTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca         >  2:2977561/1‑76 (MQ=255)
                                                                cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat        >  1:3452603/1‑76 (MQ=255)
                                                                cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat        <  1:4691790/76‑1 (MQ=255)
                                                                cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat        >  1:5831951/1‑76 (MQ=255)
                                                                 tGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa       >  2:1742522/1‑76 (MQ=255)
                                                                  gAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa      <  1:4368153/76‑1 (MQ=255)
                                                                   aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa     <  1:3687728/76‑1 (MQ=255)
                                                                   aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa     <  1:1131775/76‑1 (MQ=255)
                                                                    aaTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa     <  2:1683510/76‑1 (MQ=255)
                                                                    aaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa    <  2:4913975/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa   <  2:810383/76‑1 (MQ=255)
                                                                      ttGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt  <  2:2596385/76‑1 (MQ=255)
                                                                           |                                                                       
GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAT  >  REL606/3911891‑3912037

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