Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,893,550 +G 4.8% intergenic (+5/‑51) kup → / → insJ‑5 potassium transporter/IS150 hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,893,5501.→G 4.8% 7.4 0.0 ‑0.162intergenic (+5/‑51)kup/insJ‑5potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/2):  ref base (26/33):  total (27/35)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.14e-01

TATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..A..GACTGTACTGCACCCATTTTGTTGGACGATGAAAT  >  REL606/3893515‑3893585
                                      |                                    
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..a                                       <  1:8296114/36‑1 (MQ=255)
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..a                                       >  1:5228505/1‑36 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCC..A..g                                    >  1:8070720/1‑36 (MQ=255)
  tCGAACTGGGTACTCAGGTCGAAATCTAACGCC..A..Ga                                   >  1:2632029/1‑36 (MQ=255)
   tgtttaactgtccacctttttggggtcagtaC..A..GAc                                  >  1:5633225‑M1/32‑36 (MQ=255)
   tgtttaactgtccaactttttggggtcagtaC..A..GAc                                  >  1:268887‑M1/32‑36 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCC..A..GAc                                  <  1:4513083/36‑1 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCC..A..GAc                                  <  1:515615/36‑1 (MQ=255)
    gtttaactgtccaactttttggggtcagtaC..A..GACt                                 <  1:5109585‑M1/6‑1 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCC..A..GACt                                 >  1:4248728/1‑36 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCC..A..GACt                                 >  1:6077975/1‑36 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACGCC..A..GACTGt                               <  1:7813515/36‑1 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACGCC..A..GACTGt                               <  1:2776400/36‑1 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACGCC..A..GACTGt                               >  1:2176603/1‑36 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACACC..A..GACTGt                               <  1:2524519/36‑1 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCC..A..GACTGTa                              >  1:4729795/1‑36 (MQ=255)
        tGGGTACTCAGGTCGAAATCTAACGCC..A..GACTGTAc                             <  1:5377606/36‑1 (MQ=255)
         gggTACTCAGGTCGAAATCTAACGCC..A..GACTGTACt                            <  1:991601/36‑1 (MQ=255)
          ggTACTCAGGTCGAAATCTAACGCC..A..GACTGTACTg                           <  1:2925741/36‑1 (MQ=255)
           gTACTCAGGTCGAAATCTAACGCC..A..GACTGTACTGc                          >  1:1858315/1‑36 (MQ=255)
            tACTCAGGTCGAAATCTAACGCC..AG.GACTGTACTGc                          <  1:5999114/36‑1 (MQ=37)
            tACTCAGGTCGAAATCTAACGCC..A..GACTGTACTGCa                         <  1:7968813/36‑1 (MQ=255)
            tACTCAGGTAGAAATCTATCGCA..A..GACTGTACTGCa                         >  1:1452899/1‑36 (MQ=25)
             aCTCAGGTCGAAATCTAACGCC..AGGGACTGTACTGc                          <  1:1355945/36‑1 (MQ=25)
             aCTCAGGTCGAAATCTAACGCC..A..GACTGTACTGCAc                        <  1:6850951/36‑1 (MQ=255)
             aCTCAGGTCGAAATCTAACGCC..A..GACTGTACTGCAc                        <  1:7519249/36‑1 (MQ=255)
               caactttttggggtcagtaC..A..GACTGTACTGCAccc                      >  1:7595546‑M1/20‑36 (MQ=255)
                cAGGTCGAAATCTAACGCC..A..GACTGTACTGCACCCa                     <  1:4292360/36‑1 (MQ=255)
                 aGGTCGAAATCTAACGCC..A..GACTGTACTGCACCCAt                    >  1:7994305/1‑36 (MQ=255)
                 aGGCCGAAATCTAACGCC..A..GACTGTACTGCACCCAt                    <  1:3467933/36‑1 (MQ=255)
                   gTCGAAATCTAACGCC..A..GACTGTACTGCACCCAttt                  >  1:1596533/1‑36 (MQ=255)
                   gTCGAAATCAAACGCC..A..GACTGTACTGCACCCAttt                  <  1:2349551/36‑1 (MQ=255)
                    tCGAAATCTAACGCC..A..GACTGTACTGCCCCCAtttt                 <  1:13831/36‑1 (MQ=255)
                     cGAAATCTAACGCC..A..GACTGTACTGCACCCATTttg                <  1:5766651/36‑1 (MQ=255)
                      gAAATCTAACGCC..A..GACTGTACTGCACCCATTttgt               <  1:5788687/36‑1 (MQ=255)
                       aaaTCTAACGCC..A..GACTGTACTGCACCCATTttgtt              <  1:8398735/36‑1 (MQ=255)
                       aaaTCTAACGCC..A..GACTGTACTGCACCCATTttgtt              >  1:2436293/1‑36 (MQ=255)
                        aaTCTAACGCCAGA..GACTGTACTGCACCCATTttgt               <  1:6215601/36‑1 (MQ=17)
                         gggtcagtaC..A..GACTGTACTGCACCCATTTTGTTgg            <  1:2817926‑M1/27‑1 (MQ=255)
                         aTCTAACGCC..A..GACTGTACTGCACCCATTTTGTTgg            >  1:667134/1‑36 (MQ=255)
                         aTCTAACGCC..A..GACTGTACTGCACCCATTTTGTTgg            <  1:232999/36‑1 (MQ=255)
                         aTCTAACGCC..A..GACTGTACTGCACCCATTTTGTTgg            >  1:4831787/1‑36 (MQ=255)
                          tCTAACGCC..A..GACTGTACTGCACCCATTTTGTTGGa           >  1:5388018/1‑36 (MQ=255)
                          tCTAACGCC..A..GACTGTACTGCACCCATTTTGTTGGa           <  1:4853880/36‑1 (MQ=255)
                          tCTAACGCC..A..GACTGTACTGCACCCATTTTGTTGGa           >  1:3107208/1‑36 (MQ=255)
                          tCTAACGCC..A..GACTGTACTGCACCCATTTTGTTGGa           >  1:6850058/1‑36 (MQ=255)
                           cTAATGCC..A..GACTGTACTGCACCCATTTTGTTGGAc          <  1:8263360/36‑1 (MQ=255)
                           cTAACGCC..A..GACTGTACTGCGCCCATTTTGTTGGAc          <  1:356076/36‑1 (MQ=255)
                             aaCGCC..A..GACTGTACTGCACCCATTTTGTTGGACGa        <  1:3519170/36‑1 (MQ=255)
                              aCGCC..A..GACTGTACTGCACCCGCTTTGTTGGACGAt       <  1:8332064/36‑1 (MQ=18)
                              aCGCC..A..GACTGTACTGCACCCATTTTGTTGGACGAt       <  1:1074364/36‑1 (MQ=255)
                               cGCC..AG.GACTGTACTGCACCCATTTTGTTGGACGAt       >  1:3684187/1‑36 (MQ=11)
                                gCC..A..GACTGTACTGCACCCATTTTGTTGGACGATGa     <  1:7201701/36‑1 (MQ=255)
                                gCC..A..GACTGTAATGCCCCCATTTTGTTGGACGATGa     <  1:7533339/36‑1 (MQ=17)
                                ccc..A..GACTGTACTACACCCATTTTGTTGGACGATGa     <  1:7802520/35‑1 (MQ=255)
                                 cc..A..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:5247368/1‑36 (MQ=255)
                                 cc..A..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:5055145/1‑36 (MQ=255)
                                 cc..A..GACTGTACTGCACCCATTTTGTTGGACGATGaa    <  1:4780760/36‑1 (MQ=255)
                                 cc..A..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:3602261/1‑36 (MQ=255)
                                 cc..A..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:3096512/1‑36 (MQ=255)
                                 cc..A..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:1724356/1‑36 (MQ=255)
                                  cagA..GACTGTACTGCACCCATTTTGTTGGACGATGa     <  1:4123569‑M1/33‑1 (MQ=255)
                                  C..A..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:2892839‑M1/1‑36 (MQ=255)
                                  C..A..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:1535233‑M1/1‑36 (MQ=255)
                                     A..GACTGTACTGCACCCATTTTGTTGGACGATTAAAt  >  1:2878227‑M1/1‑36 (MQ=255)
                                      |                                    
TATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..A..GACTGTACTGCACCCATTTTGTTGGACGATGAAAT  >  REL606/3893515‑3893585

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG