Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,963,239 G→T 8.9% A295E (GCA→GAA)  araF ← L‑arabinose transporter subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,963,2390G→T 8.9% 3.1 ‑3.5 0.061884araFL‑arabinose transporter subunit
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (19/36):  total (24/37)
Fisher's exact test strand distribution p-value = 3.07e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAGCATTTCGCTGGATTTATA  >  REL606/1963204‑1963273
                                   |                                  
tAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTg                                    <  1:6248941/36‑1 (MQ=255)
  aCTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCt                                  <  1:7723518/36‑1 (MQ=255)
  aCTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCt                                  <  1:3321096/36‑1 (MQ=255)
  aCTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCt                                  <  1:4539482/36‑1 (MQ=255)
   cTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTa                                 <  1:340255/36‑1 (MQ=255)
   cTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTa                                 <  1:1713334/36‑1 (MQ=255)
    ttCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAc                                <  1:7436784/36‑1 (MQ=255)
     tCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAcc                               <  1:4475181/36‑1 (MQ=255)
      cGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc                              <  1:6387124/36‑1 (MQ=255)
      cGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc                              <  1:6721206/36‑1 (MQ=255)
       ggTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCa                             <  1:6569783/36‑1 (MQ=255)
        gTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAg                            <  1:6020900/36‑1 (MQ=255)
         tAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGt                           <  1:516817/36‑1 (MQ=255)
         tAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGt                           <  1:7405653/36‑1 (MQ=255)
          aaaTTTTGGCGGTTCTACGTCTTTTTTTACCCCGtt                          >  1:1739797/1‑36 (MQ=19)
          aaaTTTTGGCGGTTCAACGTCTTTTGCTCCCCAGtt                          <  1:2847265/36‑1 (MQ=255)
          aaaTTTTGGCGGTTCAACGTCTTTTGCTACCCAGtt                          <  1:7016702/36‑1 (MQ=255)
          aaaTTTTGGCGGTTCAACGTCTTTTGCTACCCAGtt                          <  1:1837458/36‑1 (MQ=255)
          aaaTTTTGGCGGTTCAACGTCTTTTGCTACCCAGtt                          <  1:1548698/36‑1 (MQ=255)
           aaTTTTGGCGGTTCTACGTCTTTTTCTTCCCCGttt                         >  1:3547734/1‑35 (MQ=16)
           aaTTTTGGCGGTTCAACGTCTTTTTCTACCCCGttt                         >  1:7486543/1‑35 (MQ=25)
           aaTTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTg                         <  1:3677341/36‑1 (MQ=255)
            aTTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGt                        <  1:3001442/36‑1 (MQ=255)
              tttGGCGGTTCAACGTCTTTTGCTTCCCAGTTGTaa                      <  1:4245118/36‑1 (MQ=255)
              tttGGCGGTTCAACGTCTTTTGCTACCCAGTTGTaa                      <  1:4823323/36‑1 (MQ=255)
               ttGGCGGTTCAACGTCTTTTTCTACCCCGTTGTaaa                     >  1:5175756/1‑36 (MQ=37)
               ttGGCGGTTCAACGTCTTTTGCTACCCAGTTGTaaa                     <  1:3430953/36‑1 (MQ=255)
                tGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAg                    <  1:3524710/36‑1 (MQ=255)
                tGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAg                    <  1:1815584/36‑1 (MQ=255)
                 ggCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAGc                   <  1:5634214/36‑1 (MQ=255)
                 ggCGGTTCAACGGCTTTTTGTACCCAGTTGTGAAGc                   >  1:1449748/1‑36 (MQ=19)
                  gCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAGCa                  >  1:1925122/1‑36 (MQ=255)
                   cGGTTCAACGTCTTTTGCTACCCAGTTGTAAAGCAt                 >  1:5351157/1‑36 (MQ=255)
                   cGGTTCAACGTCTTTTGCTACCCAGTTGTAAAGCAt                 >  1:2430941/1‑36 (MQ=255)
                   cGGTTCAACGTCTTTTGCTACCCAGTTGTAAAGCAt                 <  1:1041367/36‑1 (MQ=255)
                    ggTTCAACGTCTTTTGCTACCCAGTTGTAAAGCAtt                >  1:4643922/1‑36 (MQ=255)
                    ggTTCAACGTCTTTTGCTACCCAGTTGTAAAGCAtt                >  1:6466858/1‑36 (MQ=255)
                     gTTCAACGTCTTTTGCTACCCAGTTGTAAAGCAttt               <  1:6609993/36‑1 (MQ=255)
                     gTTCAACGTCTTTTGCTACCCAGTTGTAAAGCAttt               >  1:4774969/1‑36 (MQ=255)
                        cAACGTCTTTTGCTACCCAGTTGTAAAGCATTTCGc            <  1:1699500/36‑1 (MQ=255)
                        cAACGTCTTTTGCTACCCAGTTGTAAAGCATTTCGc            <  1:7276294/36‑1 (MQ=255)
                        cAACGTCTTTTGCTACCCAGTTGTAAAGCATTTCGc            <  1:7819777/36‑1 (MQ=255)
                         aaCGTCTTTTGCTACCCAGTTGTAAAGCATTTCGCt           <  1:1442206/36‑1 (MQ=255)
                          aCGTCTTTTGCTACCCAGTTGTAAAGCATTTCGCTg          >  1:4374177/1‑36 (MQ=255)
                           cGTCTTTTGCTACCCAGTTGTAAAGCATTTCGCTgg         >  1:2203924/1‑36 (MQ=255)
                           cGTCTTTTGCTACCCAGTTGTAAAGCATTTCGCTgg         >  1:5676242/1‑36 (MQ=255)
                           cGTCTTTTGCTACCCAGTTGTAAAGCATTTCGCTgg         <  1:6296183/36‑1 (MQ=255)
                           cGTCTTTTGCTACCCAGTTGTAAAGCATTTCGCTgg         >  1:7320338/1‑36 (MQ=255)
                            gTCTTTTTCTCCCCATTTGTAAAGCATTTCGCTGGa        <  1:4918455/36‑1 (MQ=25)
                            gTCTTTTGCTACCCAGTTGTAAAGCATTTCGCTGGa        <  1:856376/36‑1 (MQ=255)
                            gTCTTTTGCTACCCAGTTGTAAAGCATTTCGATGGa        >  1:3108463/1‑36 (MQ=255)
                             tCTTTTGCTACCCAGTTGTAAAGGATTTCGCTGGAt       >  1:3367871/1‑36 (MQ=255)
                             tCTTTTGCTACCCAGTTGTAAAGCATTTCGCTGGAt       <  1:3074039/36‑1 (MQ=255)
                             tCTTTTGCTACCCAGTTGTAAAGCATTTCGCTGGAt       <  1:7888848/36‑1 (MQ=255)
                             tCTTTTGCTACCCAGTTGTAAAGCATTTCGCTGGAt       >  1:4163352/1‑36 (MQ=255)
                               ttttGCTACCTAGTTGTAAAGCATTTCGCTGGAttt     >  1:5089720/1‑36 (MQ=255)
                               ttttGCTACCCAGTTGTAAAGCCTTTCGCTGGAttt     <  1:856/36‑1 (MQ=255)
                               ttttGCTACCCAGTTGTAAAGCATTTCGCTGGAttt     >  1:1395945/1‑36 (MQ=255)
                               ttttGCTACCCAGTTGTAAAGCATTTCGCTGGAttt     >  1:2260847/1‑36 (MQ=255)
                               ttttGCTACCCAGTTGTAAAGCATTTCGCTGGAttt     >  1:6896798/1‑36 (MQ=255)
                                 ttGCTACCCAGTTGGAAAGCATTTCGCTGGATTtat   >  1:5579714/1‑36 (MQ=255)
                                  tGCTACCCAGTTGTAAAGCATTTCGCTGGATTtata  >  1:5340204/1‑36 (MQ=255)
                                   |                                  
TAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAGCATTTCGCTGGATTTATA  >  REL606/1963204‑1963273

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG