Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,226,149 G→C 6.8% G145A (GGC→GCC)  malK → fused maltose transport subunit, ATP‑binding component of ABC superfamily/regulatory protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,226,1490G→C 6.8% 3.2 0.0 ‑0.244434malKfused maltose transport subunit, ATP‑binding component of ABC superfamily/regulatory protein
Reads supporting (aligned to +/- strand):  new base (2/1):  ref base (20/21):  total (22/22)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.70e-01

CTCTCCGGTGGTCAGCGTCAGCGTGTGGCGATTGGCCGTACGCTGGTGGCCGAGCCAAGCGTATTTTTGC  >  REL606/4226115‑4226184
                                  |                                   
ctctCCGGTGGTCAGCGTCAGCGTGTGGCGATTGGc                                    <  1:3095428/36‑1 (MQ=255)
ctctCCGGTGGTCAGCGTCAGCGAGTGGCGATTGGc                                    >  1:66462/1‑36 (MQ=255)
 tctcGGGTGGTCAGCGTGAGCGTGTGGCGATTGGcc                                   <  1:1796862/36‑1 (MQ=21)
 tctcCGGTGGTCAGCGTCAGCGTGTGGCGATTGGcc                                   <  1:6951326/36‑1 (MQ=255)
    ccGGTGGTCAGCGTCAGCGTGTGGCGATTGGCCGTa                                <  1:123734/36‑1 (MQ=255)
     cGGTGGTCAGCGTCAGCGTGTGTCGATTGGCCGTAc                               <  1:3211204/36‑1 (MQ=255)
      ggtggtCAGCGTCAGCGTGTGGCGATTGGCCGTACg                              <  1:7193019/36‑1 (MQ=255)
      ggtggtCAGCGTCAGCGTGTGGCGATTGGCCGTACg                              >  1:5032789/1‑36 (MQ=255)
      ggtggtCAGCGTCAGCGTGTGGCGATTGGCCGTACg                              >  1:3743314/1‑36 (MQ=255)
      ggtggtCAGCGTCAGCGTGTGGCGATTGGCCGTACg                              >  1:2934319/1‑36 (MQ=255)
        tggtCAGCGTCAGCGTGTGGCGATTGGCCGTACGCt                            <  1:4925107/36‑1 (MQ=255)
        tggtCAGCGTCAGCGTGTGGCGATTGGCCGTACGCt                            >  1:6734359/1‑36 (MQ=255)
         ggtCAGCGTCAGCGTGTGGCGATTGGCCGTACGCtg                           <  1:4740630/36‑1 (MQ=255)
         ggtCAGCGTCAGCGTGTGGCGATTGGCCGTACGCtg                           >  1:3290519/1‑36 (MQ=255)
          gtcagcgtcagcGTGTGGCGATTGGCCGTACGCtgg                          >  1:7251430/1‑36 (MQ=255)
          gggagcggcagcgtGTGGCGATTGCCCGTGCGCtgg                          <  1:2045627/33‑1 (MQ=11)
              gcgtcagcgtGTGGCGATTGGCCGTACGCTGGTGGc                      >  1:2445099/1‑36 (MQ=255)
               ggtcagcgtGTGGCGATTGGCCTTACGCTGGTGGcc                     <  1:793289/35‑1 (MQ=255)
               cgtcagcgtGTGGCGATTGGCCGTACGCTGGTGGcc                     <  1:5444508/36‑1 (MQ=255)
               cgccagcgtGTGGCGATTGGCCGTACGCTGGTGGcc                     <  1:863526/33‑1 (MQ=255)
                gtcagcgtGTGGCGATTGGCCGTACGCTGGTGGCCg                    >  1:4467497/1‑36 (MQ=255)
                gtcagcgtGTGGCGATTGGCCGTACGCTGGTGGCCg                    >  1:5609819/1‑36 (MQ=255)
                 tcagcgtGTGGCGATTTGCCGTACGCTGGTGGCCGa                   >  1:2889450/1‑36 (MQ=255)
                   agcgtGTGGCGATTGGCCGTACGCTGGTGGCCGAGc                 >  1:6911240/1‑36 (MQ=255)
                   agcgtGTGGCGATTGGCCGTACGCTGGTGGCCGAGc                 >  1:4370113/1‑36 (MQ=255)
                   agcgtGTGGCGATTGCCCGCGCGCTGGCGGGGGAGc                 >  1:3204512/1‑36 (MQ=14)
                   agcgtGTGGCGATTGCCCGCGCGCTGGCGGGGGAGc                 >  1:5324786/1‑36 (MQ=14)
                    gcgtGTTGCGATTGGCCGTACGCTGGTGGCCGAGcc                <  1:188816/36‑1 (MQ=255)
                    gcgtGTGGCGATTGGCCGTACGCTGGTGGCCGAGcc                <  1:385031/36‑1 (MQ=255)
                     cgtGTGGCGATTGGCCGTACGCTGGTGGCCGAGCCa               >  1:3107908/1‑36 (MQ=255)
                      gtgtGGCGATTGGCCGTACGCTGGTGGCCGAGCCaa              <  1:3611878/36‑1 (MQ=255)
                      gtgtGGCGATTGGCCGTACGCTGGTGGCCGAGCCaa              <  1:525163/36‑1 (MQ=255)
                       tgtgGCGATTGGCCGTACGCTGGTGGCCGAGCCAAg             <  1:4657598/36‑1 (MQ=255)
                       tgtgGCGATTGGCCGTACGCTGGTGGCCGAGCCAAg             <  1:2899731/36‑1 (MQ=255)
                        gtgGCGATTGGCCGTACGCTGGTGGCCGAGCCAAGc            >  1:4951802/1‑36 (MQ=255)
                        gtgGCGATTGGCCGTACGCTGGTGGCCGAGCCAAGc            <  1:3641336/36‑1 (MQ=255)
                        gtgGCGATTGGCCGTACGCTGGTGGCCGAGCCAAGc            <  1:3020009/36‑1 (MQ=255)
                         tgGCGATTGGCCGTACGCTGGTGGCCGAGCCAAGCg           >  1:3622845/1‑36 (MQ=255)
                          ggCGATTGGCCGTACGCTGGTGGCCGAGCCAAGCGt          >  1:1646689/1‑36 (MQ=255)
                          ggCGATTGGCCGTACGCTGGTGGCCGAGCCAAGCGt          >  1:2941886/1‑36 (MQ=255)
                          ggCGATTGGCCGTACGCTGGTGGCCGAGCCAAGCGt          <  1:839260/36‑1 (MQ=255)
                             gATTGGCCGTACGCTGGTGGCCGAGCCAAGCGTAtt       >  1:1857072/1‑36 (MQ=255)
                             gATTGGCCGTACGCTGGTGGCCGAGCCAAGCGTAtt       >  1:356711/1‑36 (MQ=255)
                              atgggCCGTACGCTGGTGGCCGAGCCAAGCGTAttt      <  1:2351009/33‑1 (MQ=255)
                               ttGGCCGTACGCTGGTGGCCGAGCCAAGCGTAtttt     <  1:897828/36‑1 (MQ=255)
                                  gCCGTACGCTGGTGGCCGAGCCAAGCGTATTTTTGc  >  1:4800731/1‑36 (MQ=255)
                                  gCCGTACGCTGGTGGCCGAGCCAAGCGTATTTTTGc  <  1:1579517/36‑1 (MQ=255)
                                  |                                   
CTCTCCGGTGGTCAGCGTCAGCGTGTGGCGATTGGCCGTACGCTGGTGGCCGAGCCAAGCGTATTTTTGC  >  REL606/4226115‑4226184

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG