Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,911,966 Δ1 bp 26.1% intergenic (‑167/‑156) yifB ← / → ilvL predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,911,9660T→.26.1% 19.6 0.0 ‑0.123intergenic (‑167/‑156)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (2/15):  total (2/21)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.66e-01

ACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACA  >  REL606/3911931‑3911999
                                   |                                 
aCGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                   <  1:3091158/36‑1 (MQ=255)
  gAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                 >  1:1833122/1‑36 (MQ=255)
      aaGGGAACAACATTCATACTGAAATTGAAttttttt                             <  1:3171777/36‑1 (MQ=255)
       aGGTAACAGCATTCATACTGAAATTGAAttttttcg                            <  1:5929105/36‑3 (MQ=25)
       aGGGAACAACATTCATACTGAAATTGAATTTTTTtc                            <  1:2156440/36‑1 (MQ=255)
           aacaacATTCATACTGAAATTGAATTTTTTtcactc                        <  1:6891450/36‑1 (MQ=255)
            acaacaTTCATACTGAAATTGAATTTTTTtcactca                       <  1:4956491/36‑1 (MQ=255)
            acaacaTTCATACTGAAATTGAA.TTTTTtcactcac                      <  1:1299665/36‑1 (MQ=37)
            aaaacaTTCATACTGAAATTGAA.TTTTTtcactcac                      <  1:7579796/34‑1 (MQ=37)
             caacaTTCATACTGAAATTGAA.TTTTTTcactcact                     <  1:1757060/36‑1 (MQ=37)
             aaacaTTCATACTGAAATTGAATTTTTTtcactcac                      <  1:1778208/35‑1 (MQ=255)
              aacaTTCATACTGAAATTGAATTTTTTTcactcact                     <  1:3261338/36‑1 (MQ=255)
                caTTCATACTGAAATTGAATTTTTTTCACTCACtat                   <  1:7478651/36‑1 (MQ=255)
                caTTCATACTGAAATTGAATTTTTTTCACTCACtat                   <  1:3253909/36‑1 (MQ=255)
                 aTTCATACTGAAATTGAATTTTTTTCACTCACtatt                  <  1:7642510/36‑1 (MQ=255)
                  ttCATACTGAAATTGAA.TTTTTTCACTCACtatttt                <  1:3748270/36‑1 (MQ=37)
                   tCATACTGAAATTGAATTTTTTTCACTCACtatttt                >  1:6661237/1‑36 (MQ=255)
                   tCATACTGAAATTGAATTTTTTTCACTCACtatttt                <  1:1281439/36‑1 (MQ=255)
                    cATACTGAAATTGAATTTTTTTCACTCACtatttta               <  1:7146703/36‑1 (MQ=255)
                    cATACGGAAATTGAA.TTTTTTCACTCACtattttat              <  1:2942732/36‑1 (MQ=25)
                     aTACTGAAATTGAATTTTTTTCACTCACtattttat              <  1:7399524/36‑1 (MQ=255)
                        cTGAAATTGAATTTTTTTCACTCACTattttatttt           <  1:1427854/36‑1 (MQ=255)
                           aaaTTGAATTTTTTTCACTCACTATTTTATTTTTaa        <  1:6330867/36‑1 (MQ=255)
                            aaTTGAATTTTTTTCACTCACTATTTTATTTTTaaa       <  1:6261067/36‑1 (MQ=255)
                              ttGAATTTTTTTCACTCACTATTTTATTTTTaaaaa     <  1:4026996/36‑1 (MQ=255)
                              ttGAA.TTTTTTCACTCACTATTTTATTTTTaaaaaa    <  1:4614293/36‑1 (MQ=37)
                                 aaTTTTTTTCACTCACTATTTTATTTTTAAAAaaca  >  1:5396752/1‑36 (MQ=255)
                                   |                                 
ACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACA  >  REL606/3911931‑3911999

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 35 ≤  ATCG  < 37 ≤  ATCG