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.89.4% 179.0 / 19.7 47intergenic (‑167/‑156)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base . (13/29);  total (13/34)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAT  >  REL606/3911891‑3912037
                                                                           |                                                                       
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                         <  1:179163/76‑1 (MQ=255)
 aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                        >  1:1319516/1‑76 (MQ=255)
   tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                      >  2:712850/1‑76 (MQ=255)
   tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                      >  2:712832/1‑76 (MQ=255)
   tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                      <  2:1432783/76‑1 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                    >  1:1072004/1‑76 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                    <  1:80882/76‑1 (MQ=255)
      tgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttc                                                                   <  2:74251/76‑2 (MQ=255)
        cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttcac                                                                 >  1:261881/1‑73 (MQ=255)
            cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactcac                                                            <  1:1063074/76‑1 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTcactcact                                                           >  2:931373/1‑76 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTcactcact                                                           <  1:1762029/76‑1 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTcactcact                                                           <  1:1278279/76‑1 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTcactcact                                                           <  1:1261223/76‑1 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtat                                                         <  2:1028851/76‑1 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtat                                                         <  1:759758/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtat                                                         <  2:1734148/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                      <  2:273998/76‑1 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                      >  1:226822/1‑76 (MQ=255)
                  tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                      <  1:1329031/76‑1 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                    >  2:68856/1‑76 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTAttttt                                                <  2:1259332/76‑1 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTAttttt                                                >  2:21827/1‑76 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTAttttt                                                <  2:1056376/76‑1 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                             >  1:273731/1‑76 (MQ=255)
                               ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                         <  1:1126178/76‑1 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                     <  1:1299646/76‑1 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                     <  1:271242/76‑1 (MQ=255)
                                    aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                     <  1:1754885/76‑1 (MQ=255)
                                    aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                    <  1:563339/76‑1 (MQ=255)
                                     aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                   <  2:1429860/76‑1 (MQ=255)
                                      caaCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  <  2:1394935/74‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  <  1:1667924/76‑1 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                 <  1:1743421/76‑1 (MQ=255)
                                         cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                               <  2:522216/76‑1 (MQ=255)
                                         cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                               >  2:1129723/1‑76 (MQ=255)
                                         cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                               >  1:489959/1‑76 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                              <  2:632568/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                           >  2:572485/1‑76 (MQ=255)
                                               aGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                         >  1:492310/1‑76 (MQ=255)
                                                  gAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                      >  1:383863/1‑76 (MQ=255)
                                                     caacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                   >  2:1017221/1‑76 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                  >  2:1059869/1‑76 (MQ=255)
                                                       acaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa                 <  1:402657/76‑1 (MQ=255)
                                                         aTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc               <  2:434042/76‑1 (MQ=255)
                                                         aTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc               <  1:705407/76‑1 (MQ=255)
                                                          ttCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                           >  2:546037/1‑63 (MQ=255)
                                                          ttCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                           <  1:546037/63‑1 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc             <  1:243227/76‑1 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc             <  2:381721/76‑1 (MQ=255)
                                                               aCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc          <  2:122092/76‑1 (MQ=255)
                                                                   aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa     <  1:1123167/76‑1 (MQ=255)
                                                                    aaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa    <  1:495771/76‑1 (MQ=255)
                                                                    aaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa    <  1:174509/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa   <  2:399428/76‑1 (MQ=255)
                                                                       tGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt  <  1:198328/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 < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: