Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,964 2 bp→TT intergenic (‑165/‑157) 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,9640AT95.2% 156.2 / ‑0.8 42intergenic (‑165/‑158)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (16/24);  total (17/25)
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.79e-01
*REL6063,911,9650AT95.2% 157.7 / ‑0.7 42intergenic (‑166/‑157)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (16/24);  total (17/25)
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.18e-01

TGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA  >  REL606/3911894‑3912034
                                                                      ||                                                                     
tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGtttttt                                                                   <  2:1112164/76‑1 (MQ=255)
       cGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTtctcc                                                            >  1:209137/1‑73 (MQ=255)
           aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTcactcact                                                        <  1:796638/76‑1 (MQ=255)
            aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACta                                                       >  2:402873/1‑76 (MQ=255)
              aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACtatt                                                     >  1:1234795/1‑76 (MQ=255)
                 cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACtatttta                                                  >  1:413494/1‑76 (MQ=255)
                     ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTattttatttt                                              >  1:823162/1‑76 (MQ=255)
                      gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTAttttt                                             <  1:1447913/76‑1 (MQ=255)
                         aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTaaa                                          <  1:150552/76‑1 (MQ=255)
                            ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTaaaaaa                                       >  2:1203783/1‑76 (MQ=255)
                            ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTaaaaaa                                       >  1:1049317/1‑76 (MQ=255)
                             gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAaac                                      <  2:1490244/76‑1 (MQ=255)
                              aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAaaca                                     <  2:103899/76‑1 (MQ=255)
                               tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                    <  2:1247413/76‑1 (MQ=255)
                               tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                    <  1:665505/76‑1 (MQ=255)
                                gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAaacaac                                   <  1:1547395/76‑1 (MQ=255)
                                  aCTACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                 <  1:754284/76‑1 (MQ=255)
                                    tACGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                               <  1:306521/76‑1 (MQ=255)
                                     aCGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                              <  1:926820/76‑1 (MQ=255)
                                     aCGAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                              <  1:56330/76‑1 (MQ=255)
                                       gAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  1:730472/76‑1 (MQ=255)
                                       gAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  2:345745/76‑1 (MQ=255)
                                       gAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  2:678853/76‑1 (MQ=255)
                                       gAGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  1:1111426/76‑1 (MQ=255)
                                        aGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                           <  1:1047284/76‑1 (MQ=255)
                                        aGGAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                           <  2:1060550/76‑1 (MQ=255)
                                         tgAAGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                          >  2:1140690/2‑76 (MQ=255)
                                           aaGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        >  1:885167/1‑76 (MQ=255)
                                           aaGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        <  2:503207/76‑1 (MQ=255)
                                           aaGGGAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        >  1:829107/1‑76 (MQ=255)
                                             gggAACAACATTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                      >  2:1612752/1‑76 (MQ=255)
                                                 acaacaTTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                  >  1:1335138/1‑76 (MQ=255)
                                                   aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                >  2:366857/1‑76 (MQ=255)
                                                     caTTCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa              <  1:237968/76‑1 (MQ=255)
                                                     caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa              <  2:553993/76‑1 (MQ=255)
                                                       ttCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc            <  2:1399482/76‑1 (MQ=255)
                                                       ttCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc            >  2:706634/1‑76 (MQ=255)
                                                        tCATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg           <  1:956394/76‑1 (MQ=255)
                                                         cATACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc          <  1:1154509/76‑1 (MQ=255)
                                                           tACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        >  1:964614/1‑76 (MQ=255)
                                                           tACTGAAATTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        >  1:1383434/1‑76 (MQ=255)
                                                                aaaTTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa   >  2:209415/1‑76 (MQ=255)
                                                                 aaTTGTTTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  <  2:1324657/76‑1 (MQ=255)
                                                                      ||                                                                     
TGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA  >  REL606/3911894‑3912034

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