Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,997 (A)6→7 intergenic (‑198/‑125) 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,9911.A98.3% 231.7 / ‑2.6 59intergenic (‑192/‑131)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (29/29);  total (30/29)
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.45e-01

GACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAT  >  REL606/3911916‑3912061
                                                                            |                                                                      
gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAttttt                                                                         <  2:333006/76‑1 (MQ=255)
gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAttttt                                                                         >  1:315239/1‑76 (MQ=255)
gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAttttt                                                                         <  1:2087315/76‑1 (MQ=255)
 aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTa                                                                       <  1:1302911/76‑1 (MQ=255)
   aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaa                                                                     >  2:1908056/1‑76 (MQ=255)
   aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaa                                                                     >  1:365532/1‑76 (MQ=255)
     aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaa                                                                   <  1:1743754/76‑1 (MQ=255)
      ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaaa                                                                  <  2:1064986/76‑1 (MQ=255)
        aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaaaac                                                                >  2:66514/1‑74 (MQ=255)
             cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAAcaacaa                                                            <  2:1739631/76‑1 (MQ=255)
             cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAAcaacaa                                                            >  2:1523211/1‑76 (MQ=255)
              tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAAt                                                           <  1:228411/76‑1 (MQ=255)
               aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAAtt                                                          <  1:1710520/76‑1 (MQ=255)
                cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAAttt                                                         <  1:1480558/76‑1 (MQ=255)
                  aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTtat                                                       <  2:622733/76‑1 (MQ=255)
                  aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTtat                                                       >  2:135730/1‑76 (MQ=255)
                     aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATAtt                                                    >  2:1148646/1‑76 (MQ=255)
                     aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATAtt                                                    >  1:1193086/1‑76 (MQ=255)
                     aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATAtt                                                    >  1:685383/1‑76 (MQ=255)
                      aGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTg                                                   <  2:1280803/76‑1 (MQ=255)
                      aGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTg                                                   <  1:1146436/76‑1 (MQ=255)
                      aGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTg                                                   >  1:578701/1‑76 (MQ=255)
                       gggAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGa                                                  >  2:560961/1‑76 (MQ=255)
                       gggAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGa                                                  >  1:997192/1‑76 (MQ=255)
                        ggAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGaa                                                 <  1:152163/76‑1 (MQ=255)
                        ggAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGaa                                                 <  1:128696/76‑1 (MQ=255)
                         gAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGaaa                                                <  2:1563267/76‑1 (MQ=255)
                         gAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGaaa                                                >  2:1704356/1‑76 (MQ=255)
                         gAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGaaa                                                >  1:1509676/1‑76 (MQ=255)
                             aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAattat                                            <  1:596272/76‑1 (MQ=255)
                             aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAattat                                            >  2:1531200/1‑76 (MQ=255)
                             aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAattat                                            <  1:1903041/76‑1 (MQ=255)
                              acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAattatt                                           <  1:920796/76‑1 (MQ=255)
                               caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAAttatta                                          >  1:1837755/1‑76 (MQ=255)
                                aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaa                                         <  1:411020/76‑1 (MQ=255)
                                aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaa                                         >  1:1214293/1‑76 (MQ=255)
                                  tCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAAc                                       >  2:1614573/1‑76 (MQ=255)
                                   cATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACg                                      <  1:331276/76‑1 (MQ=255)
                                       cTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                                  >  2:707909/1‑76 (MQ=255)
                                           aaTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa                              <  1:732758/76‑1 (MQ=255)
                                           aaTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa                              <  1:2079423/76‑1 (MQ=255)
                                           aaTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa                              <  1:2077201/76‑1 (MQ=255)
                                            aTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                             <  2:1696593/76‑1 (MQ=255)
                                            aTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                             >  2:2000432/1‑76 (MQ=255)
                                             ttGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATTATaaaa                            >  2:2121544/1‑76 (MQ=255)
                                             ttGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                            <  2:1850672/76‑1 (MQ=255)
                                             ttGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                            <  1:373415/76‑1 (MQ=255)
                                              tGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa                           <  1:1933114/76‑1 (MQ=255)
                                                 aTTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCAaaaaaaatc                        <  2:1759761/76‑4 (MQ=255)
                                                    tttttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGcc                     >  1:427710/1‑76 (MQ=255)
                                                    tttttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGcc                     <  2:1349356/76‑1 (MQ=255)
                                                    tttttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGcc                     <  1:1416966/76‑1 (MQ=255)
                                                     ttttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                    >  2:119387/1‑76 (MQ=255)
                                                      tttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaa                   >  1:567765/1‑76 (MQ=255)
                                                       ttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa                  <  1:443702/76‑1 (MQ=255)
                                                       ttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa                  >  1:1800273/1‑76 (MQ=255)
                                                        tcactcacTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa                 <  2:250035/76‑1 (MQ=255)
                                                        tcactcacTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaa                >  2:86019/1‑76 (MQ=255)
                                                             cactATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAatat            >  2:1709125/1‑76 (MQ=255)
                                                               ctATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCt          >  1:1729068/1‑76 (MQ=255)
                                                                tattttatttTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt         >  1:1080372/1‑76 (MQ=255)
                                                                tattttatttTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt         >  1:991572/1‑76 (MQ=255)
                                                                 attttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg        >  1:1716929/1‑76 (MQ=255)
                                                                 attttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg        >  1:1312957/1‑76 (MQ=255)
                                                                   tttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTa      <  2:1640893/76‑1 (MQ=255)
                                                                   tttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTa      <  1:1670182/76‑1 (MQ=255)
                                                                    ttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAc     >  1:790330/1‑76 (MQ=255)
                                                                    ttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAc     >  1:352917/1‑76 (MQ=255)
                                                                    tatttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACt    <  1:1674387/73‑1 (MQ=255)
                                                                     atttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa   <  2:1808768/74‑1 (MQ=255)
                                                                     atttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa   >  1:503413/3‑76 (MQ=255)
                                                                     atttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa   >  1:373999/3‑76 (MQ=255)
                                                                     atttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa   >  1:1054877/3‑76 (MQ=255)
                                                                     atttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa   <  2:977482/74‑1 (MQ=255)
                                                                      tttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAt  >  2:1246920/2‑76 (MQ=255)
                                                                      tttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAt  >  1:659705/2‑76 (MQ=255)
                                                                            |                                                                      
GACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAT  >  REL606/3911916‑3912061

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