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.96.7% 136.5 / ‑1.8 30intergenic (‑167/‑156)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (9/20);  total (9/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAA  >  REL606/3911891‑3912036
                                                                           |                                                                      
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                        <  1:729151/76‑1 (MQ=255)
 aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                       >  1:1000490/1‑76 (MQ=255)
      tgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttc                                                                  <  2:758143/76‑2 (MQ=255)
             aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTcactcact                                                          <  1:663109/76‑1 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTGCGNGGAAGGGAAAGACATTCAGACTGAAATTGAATTTTTTCACTCACtat                                                        <  1:122418/76‑1 (MQ=255)
                aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatt                                                       <  2:1060047/76‑1 (MQ=255)
                 aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtatt                                                       <  1:879770/76‑1 (MQ=255)
                      gCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttattt                                                 <  2:494579/76‑1 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTAttttt                                               <  2:836200/76‑1 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaa                                             <  1:864117/76‑1 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                            >  2:610292/1‑76 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                            <  1:240615/76‑1 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaa                                          <  2:29603/76‑1 (MQ=255)
                              aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTGTTaaaaaa                                         <  2:164394/76‑1 (MQ=255)
                               ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                        <  1:576058/76‑1 (MQ=255)
                                  tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaac                                     >  1:1185566/1‑76 (MQ=255)
                                    aaCTAGAAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                   <  2:1335950/76‑1 (MQ=255)
                                    aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                   <  2:1237887/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                            <  1:401522/76‑1 (MQ=255)
                                             gAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                          <  1:465843/76‑1 (MQ=255)
                                              aaGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                         >  1:1102267/1‑76 (MQ=255)
                                               aGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                        >  1:903414/1‑76 (MQ=255)
                                                gggAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaaa                       <  1:1103656/76‑1 (MQ=255)
                                                 ggAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat                      <  1:254151/76‑1 (MQ=255)
                                                  gAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                     <  1:416694/76‑1 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                    >  2:818497/1‑76 (MQ=255)
                                                    acaacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                   <  1:899607/76‑1 (MQ=255)
                                                         aTTCATACTGAAATTGAANTTTTTCACTCACTATTTTATGTTTAAAGAACAACAATTTATATTGAAATTATTAAAc              >  1:1339573/1‑76 (MQ=255)
                                                             aTACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat          >  2:749577/1‑76 (MQ=255)
                                                               aCTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca        <  1:1165797/76‑1 (MQ=255)
                                                                cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat       >  1:657707/1‑76 (MQ=255)
                                                                   aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa    >  2:997751/1‑76 (MQ=255)
                                                                     aTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa  <  1:383648/76‑1 (MQ=255)
                                                                           |                                                                      
GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAA  >  REL606/3911891‑3912036

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