Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,807,994 C→T 9.1% intergenic (‑96/+54) luxS ← / ← gshA S‑ribosylhomocysteinase/gamma‑glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,807,9940CT 9.1% 194.7 / 8.7 66intergenic (‑96/+54)luxS/gshAS‑ribosylhomocysteinase/gamma‑glutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (21/39);  new base T (2/4);  total (23/43)
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.97e-01

AAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTT  >  minE/1807924‑1808060
                                                                      |                                                                  
aaaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATcc                                                                    <  1:405129/71‑1 (MQ=255)
   aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCt                                                                  >  1:1148587/1‑70 (MQ=255)
   aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCt                                                                  >  1:2863482/1‑70 (MQ=255)
   aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCt                                                                  >  1:1705919/1‑70 (MQ=255)
           aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCag                                                          <  1:900332/70‑1 (MQ=255)
           aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCag                                                          <  1:723498/70‑1 (MQ=255)
           aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCa                                                           >  1:2841036/1‑69 (MQ=255)
           aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCa                                                           >  1:2828300/1‑69 (MQ=255)
            aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGaa                                                               <  1:2814013/64‑1 (MQ=255)
            aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGaa                                                               <  1:1803226/64‑1 (MQ=255)
            aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGaa                                                               <  1:2537464/64‑1 (MQ=255)
            aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGaa                                                               <  1:2345543/64‑1 (MQ=255)
             ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCTCTGAATTCAggga                                                       <  1:2570641/71‑1 (MQ=255)
                   cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCTCTGAATTCAggga                                                       <  1:530720/65‑1 (MQ=255)
                   cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCTCTGAATTCAggga                                                       <  1:1372473/65‑1 (MQ=255)
                   cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCTCTGAATTCAggga                                                       <  1:2651309/65‑1 (MQ=255)
                   cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAt                                                 >  1:631393/1‑71 (MQ=255)
                   cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAt                                                 >  1:3590223/1‑71 (MQ=255)
                   cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAt                                                 >  1:1156614/1‑71 (MQ=255)
                      cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtt                                                <  1:1389405/69‑1 (MQ=255)
                        cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTgtt                                                               <  1:3129601/52‑3 (MQ=255)
                                tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    <  1:1933337/71‑1 (MQ=255)
                                tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    >  1:2813021/1‑71 (MQ=255)
                                 cTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    >  1:2625708/1‑70 (MQ=255)
                                  tGAGTATATGAAAGACGCGCATTTGTTATCATCATCTCTGAAtt                                                             >  1:2743370/1‑44 (MQ=255)
                                  tGAGTATATGAAAGACGCGCATTTGTTATCATCATCTCTGAAtt                                                             >  1:407692/1‑44 (MQ=255)
                                  tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    >  1:2392413/1‑69 (MQ=255)
                                  tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    >  1:3057796/1‑69 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                   <  1:3660041/68‑1 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                   <  1:3271884/68‑1 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                   <  1:1900570/68‑1 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                  <  1:2918261/69‑1 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                 <  1:3419298/70‑1 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                 <  1:3432084/70‑1 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                 <  1:3456758/70‑1 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                <  1:3484880/71‑1 (MQ=255)
                                    aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                <  1:2550997/71‑1 (MQ=255)
                                     gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                <  1:1755226/70‑1 (MQ=255)
                                         atGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                   <  1:3090550/63‑1 (MQ=255)
                                         atGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                   <  1:1618614/63‑1 (MQ=255)
                                          tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                   <  1:698450/62‑1 (MQ=255)
                                           gAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttt                         >  1:45190/1‑71 (MQ=255)
                                           gAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttt                         >  1:1897471/1‑71 (MQ=255)
                                            aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                        >  1:2206986/1‑71 (MQ=255)
                                            aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                        >  1:737189/1‑71 (MQ=255)
                                             aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                       <  1:2995433/71‑1 (MQ=255)
                                             aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                        >  1:434125/1‑70 (MQ=255)
                                               gACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                     <  1:1572661/71‑1 (MQ=255)
                                                  gcgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                   >  1:3728925/1‑54 (MQ=255)
                                                   cgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                 <  1:1024299/55‑1 (MQ=255)
                                                    gcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                 <  1:2916145/54‑1 (MQ=255)
                                                       tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctttt                  <  1:1932340/66‑1 (MQ=255)
                                                       tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctttt                  <  1:2242654/66‑1 (MQ=255)
                                                       tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctttt                  <  1:2459238/66‑1 (MQ=255)
                                                       tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctttt                  <  1:2602025/66‑1 (MQ=255)
                                                          gTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    <  1:2603497/45‑1 (MQ=255)
                                                          gTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    <  1:1627021/45‑1 (MQ=255)
                                                          gTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    <  1:2330466/45‑1 (MQ=255)
                                                          gTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    <  1:3452931/45‑1 (MQ=255)
                                                           ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    <  1:2607980/44‑1 (MQ=255)
                                                           ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                    <  1:3681731/44‑1 (MQ=255)
                                                            tATCATCATCCCTGAATTCAGAGATGAAATTTTGGc                                           >  1:822831/1‑36 (MQ=255)
                                                            tATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                  >  1:3122521/1‑45 (MQ=255)
                                                            tATCATCATCCCTGAATTCAGAGATGAAATTTGGCCACTCACGa                                  >  1:1877195/1‑44 (MQ=255)
                                                              tcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg      <  1:3521033/71‑1 (MQ=255)
                                                                atcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtgt     <  1:2665247/70‑1 (MQ=255)
                                                                   atcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtttt  <  1:3487767/70‑1 (MQ=255)
                                                                      |                                                                  
AAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTT  >  minE/1807924‑1808060

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