Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,991 (T)5→4 intergenic (‑192/‑131) 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,9870T.92.1% 590.6 / 49.6 141intergenic (‑188/‑135)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (4/7);  new base . (55/74);  total (60/81)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 159 total aligned reads displayed.

GCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAC  >  REL606/3911913‑3912058
                                                                           |                                                                       
  ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTattttatttt                                                                      <  2:1919241/76‑1 (MQ=255)
  ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTattttatttt                                                                      <  2:2043838/76‑1 (MQ=255)
    aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTattttattttaa                                                                    >  1:45372/1‑74 (MQ=255)
    aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTa                                                                    >  1:1554215/1‑76 (MQ=255)
     cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaa                                                                   <  1:785639/76‑1 (MQ=255)
     cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaa                                                                   <  2:489922/76‑1 (MQ=255)
     cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaa                                                                   >  2:4384669/1‑76 (MQ=255)
      aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaa                                                                  <  1:4372996/76‑1 (MQ=255)
      aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaa                                                                  <  2:1606090/76‑1 (MQ=255)
       aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAaaaa                                                                 <  1:778258/76‑1 (MQ=255)
         ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaaaac                                                               >  2:832949/1‑75 (MQ=255)
          gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaaca                                                             >  1:893067/1‑76 (MQ=255)
           aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaa                                                            >  1:3600658/1‑76 (MQ=255)
           aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaa                                                            <  2:4415330/76‑1 (MQ=255)
            tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaac                                                           <  1:3720377/76‑1 (MQ=255)
            tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaac                                                           <  2:975898/76‑1 (MQ=255)
            tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaac                                                           <  2:1095947/76‑1 (MQ=255)
            tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaac                                                           <  1:4224367/76‑1 (MQ=255)
             gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAacaaca                                                          <  2:2505458/76‑1 (MQ=255)
             gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAacaaca                                                          <  1:4681748/76‑1 (MQ=255)
             gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAacaaca                                                          <  1:2587245/76‑1 (MQ=255)
              aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAAcaacaa                                                         <  2:2443296/76‑1 (MQ=255)
              aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAAcaacaa                                                         >  2:4162846/1‑76 (MQ=255)
              aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAAcaacaa                                                         <  2:1006895/76‑1 (MQ=255)
               aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAAt                                                        <  1:2536743/76‑1 (MQ=255)
               aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAAt                                                        <  2:3599427/76‑1 (MQ=255)
               aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAAt                                                        <  1:1009925/76‑1 (MQ=255)
                cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAAtt                                                       >  1:4945285/1‑76 (MQ=255)
                cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAAtt                                                       <  1:4582933/76‑1 (MQ=255)
                 tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAAttt                                                      <  2:1576295/76‑1 (MQ=255)
                    gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTtata                                                   <  1:3044186/76‑1 (MQ=255)
                    gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTtata                                                   <  2:2190721/76‑1 (MQ=255)
                    gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTtata                                                   >  2:252561/1‑76 (MQ=255)
                     aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTatat                                                  >  2:4822228/1‑76 (MQ=255)
                      ggAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATAtt                                                 <  1:1587717/76‑1 (MQ=255)
                        aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGa                                               <  1:1033898/76‑1 (MQ=255)
                         aGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGaa                                              <  2:1990145/76‑1 (MQ=255)
                         aGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGaa                                              >  2:1834501/1‑76 (MQ=255)
                             aacaacATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAatta                                          >  2:4007473/1‑76 (MQ=255)
                                aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAAttatta                                       <  2:4472456/76‑1 (MQ=255)
                                aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAAttatta                                       >  1:3671680/1‑76 (MQ=255)
                                aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAAttatta                                       <  1:2937191/76‑1 (MQ=255)
                                aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAAttatta                                       >  1:561530/1‑76 (MQ=255)
                                  caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                     <  2:1775003/76‑1 (MQ=255)
                                  caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                     >  2:456896/1‑76 (MQ=255)
                                  caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                     <  2:2711779/76‑1 (MQ=255)
                                  caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                     <  1:4169787/76‑1 (MQ=255)
                                   aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAAc                                    >  2:2809802/1‑76 (MQ=255)
                                    ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACg                                   <  1:4484331/76‑1 (MQ=255)
                                     tCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc                                  <  1:61345/76‑1 (MQ=255)
                                       aTACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat                                <  1:1475813/76‑1 (MQ=255)
                                        tACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                               >  1:2060092/1‑76 (MQ=255)
                                        tACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                               >  1:2498991/1‑76 (MQ=255)
                                        tACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                               <  1:2787457/76‑1 (MQ=255)
                                         aCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca                              >  2:4278415/1‑76 (MQ=255)
                                         aCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca                              <  1:1691410/76‑1 (MQ=255)
                                          cTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                               <  1:1508911/76‑1 (MQ=255)
                                           tGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                            <  2:309910/76‑1 (MQ=255)
                                           tGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                            >  1:2684083/1‑76 (MQ=255)
                                           tGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                            <  2:4732293/76‑1 (MQ=255)
                                             aaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                          <  1:3111458/76‑1 (MQ=255)
                                             aaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                          <  2:2835829/76‑1 (MQ=255)
                                             aaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                          <  1:3302525/76‑1 (MQ=255)
                                              aaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                         <  2:1697685/76‑1 (MQ=255)
                                              aaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                         <  2:2603004/76‑1 (MQ=255)
                                               aTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa                        >  1:4680755/1‑76 (MQ=255)
                                               aTTGAATTTGTTTCGCTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa                        <  1:1239307/76‑1 (MQ=255)
                                                ttGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt                       >  2:1205952/1‑76 (MQ=255)
                                                ttGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt                       <  1:4498666/76‑1 (MQ=255)
                                                ttGAATTTTTTTCACTCACTATTTTGTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt                       >  1:5004415/1‑76 (MQ=255)
                                                 tGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATc                      >  2:920909/1‑76 (MQ=255)
                                                 tGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATc                      >  1:891672/1‑76 (MQ=255)
                                                  gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                     >  1:3404391/1‑76 (MQ=255)
                                                  gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                     >  1:4898908/1‑76 (MQ=255)
                                                  gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                     >  1:1622107/1‑76 (MQ=255)
                                                  gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                     <  2:3378945/76‑1 (MQ=255)
                                                  gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                     >  2:4510438/1‑76 (MQ=255)
                                                  gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                     >  1:185942/1‑76 (MQ=255)
                                                  gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                     >  1:3223622/1‑76 (MQ=255)
                                                   aaTTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCgg                    >  2:3303005/1‑76 (MQ=255)
                                                   aaTTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCgg                    >  2:3886056/1‑76 (MQ=255)
                                                   aaTTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCgg                    <  2:605818/76‑1 (MQ=255)
                                                    aTTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGc                   <  2:522990/76‑1 (MQ=255)
                                                       ttttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                 <  2:2722452/76‑1 (MQ=255)
                                                       ttttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                 >  1:3929752/1‑76 (MQ=255)
                                                       ttttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                 >  1:1349111/1‑76 (MQ=255)
                                                         ttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa               <  1:961249/76‑1 (MQ=255)
                                                          tttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa              >  2:3535052/1‑76 (MQ=255)
                                                          tttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa              <  2:621883/76‑1 (MQ=255)
                                                           ttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaa             <  1:669997/76‑1 (MQ=255)
                                                               ctcactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAatat         <  1:4105293/76‑1 (MQ=255)
                                                                 cactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCt       >  1:2559957/1‑76 (MQ=255)
                                                                   ctATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg     >  1:3376093/1‑76 (MQ=255)
                                                                    tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt    >  1:2238022/1‑76 (MQ=255)
                                                                    tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt    >  1:1152025/1‑76 (MQ=255)
                                                                     attttattttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt    <  1:3247538/76‑1 (MQ=255)
                                                                     attttattttAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTa   >  1:1254620/1‑76 (MQ=255)
                                                                     attttattttAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTa   >  1:1072899/1‑76 (MQ=255)
                                                                      ttttattttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTa   <  2:1021382/76‑1 (MQ=255)
                                                                      ttttattttaAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAc  >  1:1965544/1‑76 (MQ=255)
                                                                           |                                                                       
GCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAC  >  REL606/3911913‑3912058

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