Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,082,685 G→A 92.5% A494V (GCT→GTT)  yegI ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,082,6850G→A92.5% 3.9 0.0 0.0531481yegIhypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/29):  ref base (2/2):  total (22/31)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAGCATT  >  REL606/2082652‑2082717
                                 |                                
cGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACa                                >  1:4384262/1‑36 (MQ=255)
cGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACa                                <  1:7465202/36‑1 (MQ=255)
 gTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACAc                               <  1:4944403/36‑1 (MQ=255)
 gTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACAc                               <  1:6336366/36‑1 (MQ=255)
  ttCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACg                              >  1:8370386/1‑36 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             <  1:1834560/36‑1 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             >  1:3070238/1‑36 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             <  1:7818117/36‑1 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             <  1:6483190/36‑1 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             <  1:2529283/36‑1 (MQ=255)
    ctcTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtt                            <  1:7834386/34‑1 (MQ=38)
    cccTTGCGGTGTCGTGAAGTGCTTCCAGAACACGtt                            <  1:1168245/36‑1 (MQ=38)
    cccTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtt                            <  1:2086479/36‑1 (MQ=255)
      cTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtttt                          >  1:7539215/1‑36 (MQ=255)
      cTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtttt                          <  1:1832838/36‑1 (MQ=255)
       ttGCGGTGTCGTGAAGTGCTGCCAGAGCACGttttt                         >  1:4194478/1‑36 (MQ=255)
       ttGCGGTGTCGTGAAGTGCTGCCAGAACACGttttt                         <  1:498475/36‑1 (MQ=255)
       ttGCGGTGTCGTGAAGTGCTGCCAGAACACGttttt                         >  1:4292120/1‑36 (MQ=255)
        tGCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttc                        <  1:7689237/36‑1 (MQ=255)
        tGCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttc                        >  1:1889307/1‑36 (MQ=255)
         gCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttct                       >  1:609238/1‑36 (MQ=255)
         gCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttct                       >  1:1114816/1‑36 (MQ=255)
          cGGTGTCGTGAAGTGCTGCCAGAACACGTTTttctt                      <  1:1334268/36‑1 (MQ=255)
          cGGTGTCGTGAAGTGCTGCCAGAACACGTTTttctt                      <  1:830310/36‑1 (MQ=255)
           ggTGTCGTGAAGTGCTGCCAGAACACGTTTttcttc                     >  1:5376249/1‑36 (MQ=255)
            gtgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcc                    >  1:2388726/1‑36 (MQ=255)
            gtgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcc                    >  1:7924034/1‑36 (MQ=255)
             tgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcct                   <  1:6383614/36‑1 (MQ=255)
             tgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcct                   <  1:2608116/36‑1 (MQ=255)
             tgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcct                   <  1:7408024/36‑1 (MQ=255)
               tCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctcc                 <  1:2941814/36‑1 (MQ=255)
               tCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctcc                 <  1:3403547/36‑1 (MQ=255)
               tCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctcc                 <  1:3957952/36‑1 (MQ=255)
                cctGACGTGCTGCCAGCACACGTTTTTCTTCCTCCg                <  1:239801/34‑1 (MQ=18)
                 gTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCgt               >  1:5035266/1‑35 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               <  1:4561601/36‑1 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               >  1:2364092/1‑36 (MQ=255)
                     aGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTAAc           >  1:6638649/1‑36 (MQ=255)
                     aGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTAAc           <  1:6168883/36‑1 (MQ=255)
                       tgctgcCAGAACACGTTTTTCTTCCTCCGGTAAccc         <  1:7419659/36‑1 (MQ=255)
                        gctgcCAGAGCACGTTTTTCTTCCTCCGGTAACCCg        <  1:7204564/36‑1 (MQ=255)
                        gctgcCAGAGCACGTTTTTCTTCCTCCGGTAACCCg        <  1:4634174/36‑1 (MQ=255)
                        gctgcCAGAACACGTTTTTCTTCCTCCGGTAACCCg        >  1:8030366/1‑36 (MQ=255)
                        gctgcCAGAACACGTTTTTCTTCCTCCGGTAACCCg        <  1:7636836/36‑1 (MQ=255)
                         ctgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGa       <  1:3517409/36‑1 (MQ=255)
                         ctgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGa       >  1:385531/1‑36 (MQ=255)
                         ctgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGa       >  1:1557729/1‑36 (MQ=255)
                          tgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGAg      <  1:6622612/36‑1 (MQ=255)
                            ccAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCa    <  1:6065593/36‑1 (MQ=255)
                            ccAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCa    >  1:5311961/1‑36 (MQ=255)
                             cAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCAt   >  1:232393/1‑36 (MQ=255)
                             cAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCAt   >  1:2740156/1‑36 (MQ=255)
                             cAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCAt   >  1:8223256/1‑36 (MQ=255)
                              agaaCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt  <  1:5353674/36‑1 (MQ=255)
                              agaaCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt  >  1:848738/1‑36 (MQ=255)
                                 |                                
CGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAGCATT  >  REL606/2082652‑2082717

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG