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.5% 569.1 / 0.7 132intergenic (‑166/‑157)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (30/100);  total (30/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 160 total aligned reads displayed.

GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAA  >  REL606/3911891‑3912035
                                                                           |                                                                      
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                       <  1:3254467/76‑1 (MQ=255)
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                       >  1:4842098/1‑76 (MQ=255)
 aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATCCTGAAATTGAAtt                                                                      >  2:4145052/1‑76 (MQ=255)
 aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                      >  1:4840386/1‑76 (MQ=255)
  atgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                                                     <  1:4080152/76‑1 (MQ=255)
  atgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                                                     <  2:4012745/76‑1 (MQ=255)
   tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                    <  1:92934/76‑1 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                  <  1:356164/76‑1 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                  >  2:3821267/1‑76 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                  >  1:7393329/1‑76 (MQ=255)
       gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTTAAtttttttt                                                                <  2:7211224/76‑2 (MQ=255)
       gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttttt                                                                >  2:6666120/1‑75 (MQ=255)
       gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttttt                                                                <  1:4144335/76‑2 (MQ=255)
        cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttc                                                               <  1:6758459/76‑3 (MQ=255)
        cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttc                                                               <  2:5320391/76‑3 (MQ=255)
         tCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttca                                                              <  1:7021603/76‑4 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                            <  1:967113/76‑1 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                            <  1:6087993/76‑1 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                            >  1:4607944/1‑76 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                            <  1:4555087/76‑1 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactca                                                           >  2:2817707/1‑76 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactca                                                           <  2:3089032/76‑1 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactca                                                           >  1:4276260/1‑76 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactca                                                           >  1:1053659/1‑76 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactca                                                           <  1:1515275/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                        <  2:7652978/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                        <  1:3275348/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                        <  1:7784167/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                        >  2:6177194/1‑76 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACta                                                        <  1:1866120/76‑1 (MQ=255)
                 aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtatt                                                      <  2:667385/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtattt                                                     <  1:6589937/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                     <  2:4337393/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                     <  2:6002096/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                     <  2:3755311/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                     <  2:2754570/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                     <  2:2110674/76‑1 (MQ=255)
                   gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattt                                                     >  2:1371598/1‑76 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttt                                                    <  2:5917706/76‑1 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttt                                                    >  2:1298249/1‑76 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttt                                                    <  1:4801004/76‑1 (MQ=255)
                     aGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttta                                                   <  1:381867/76‑1 (MQ=255)
                     aGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttta                                                   <  1:4863554/76‑1 (MQ=255)
                      gCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttat                                                  <  1:4036763/76‑1 (MQ=255)
                       cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttatt                                                 <  1:3143781/76‑1 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttattt                                                <  2:4182611/76‑1 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttattt                                                <  1:2798753/76‑1 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTAttttt                                              <  2:7486275/76‑1 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTAttttt                                              <  2:2662125/76‑1 (MQ=255)
                           cAAAGGATAAACTACGAGGAAGGGACCAACATTCATACTGAAATTGATTTTTTTTCACTCACTATTTTATTTTTaa                                            <  2:5558217/76‑1 (MQ=255)
                            aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaa                                            <  2:195498/76‑1 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaa                                           <  2:486875/76‑1 (MQ=255)
                              aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaa                                          <  1:74034/76‑1 (MQ=255)
                              aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaa                                          <  1:3800117/76‑1 (MQ=255)
                              aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaa                                          <  1:5104126/76‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaac                                       >  1:2702822/1‑76 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaaca                                      <  1:7792812/76‑1 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaaca                                      >  2:1734846/1‑76 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaaca                                      <  2:1463077/76‑1 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaaca                                      <  2:6117484/76‑1 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                     <  2:4543866/76‑1 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                     <  1:1031169/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                   <  1:5695853/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                   <  2:1625393/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                   <  1:5479263/76‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                  <  2:1089915/76‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                  <  1:2421604/76‑1 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                 >  2:3507236/1‑76 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                 <  2:5902428/76‑1 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                 <  1:6988482/76‑1 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                 >  1:5979166/1‑76 (MQ=255)
                                         cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                               <  2:2588270/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                             <  1:759953/76‑1 (MQ=255)
                                            ggAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                            <  1:2887356/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                           <  1:2974050/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                           <  2:255528/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                           <  1:2615345/76‑1 (MQ=255)
                                                gggAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                        <  1:5274224/76‑1 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat                     >  1:3154159/1‑76 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                    <  2:3947025/76‑1 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                    <  2:6586652/76‑1 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                    <  1:1735770/76‑1 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                  >  2:3177123/1‑76 (MQ=255)
                                                        caTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                <  1:2209938/76‑1 (MQ=255)
                                                        caTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                <  1:6421951/76‑1 (MQ=255)
                                                        caTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                <  1:7500968/76‑1 (MQ=255)
                                                        caTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                <  1:5535388/76‑1 (MQ=255)
                                                         aTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa               >  1:7259309/1‑76 (MQ=255)
                                                         aTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa               >  1:1897944/1‑76 (MQ=255)
                                                            cATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg            >  1:5966924/1‑76 (MQ=255)
                                                            cATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg            <  1:7568204/76‑1 (MQ=255)
                                                              tACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca          >  2:6354731/1‑76 (MQ=255)
                                                              tACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca          <  1:2549703/76‑1 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      <  1:5191392/76‑1 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      >  2:563838/1‑76 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      <  1:4894215/76‑1 (MQ=255)
                                                                  gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat      <  1:3400042/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa   >  1:6029785/1‑76 (MQ=255)
                                                                     aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa   <  1:2567511/76‑1 (MQ=255)
                                                                      ttGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa  <  2:6928012/76‑1 (MQ=255)
                                                                           |                                                                      
GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAA  >  REL606/3911891‑3912035

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