Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,812,851 A→T 15.0% D64E (GAT→GAA yqcC ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,812,8510A→T15.0% 3.0 ‑4.1 0.048192yqcChypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (21/18):  total (30/18)
Fisher's exact test strand distribution p-value = 1.67e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCA  >  REL606/2812818‑2812886
                                 |                                   
aaCGGCAAAAGCGCCCGGTAACGGCTGTTTTTTTTc                                   >  1:410259/1‑36 (MQ=37)
aaCGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATc                                   <  1:7012994/36‑1 (MQ=255)
   ggCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTaa                                <  1:1656705/36‑1 (MQ=255)
   ggCAAAAGCGCCCGGTAACGGCTGGTTGTTTTCTat                                >  1:5917036/1‑35 (MQ=25)
     cAAAAGCGCCCGGTAACGGCTGTTTTTTATCTAACa                              >  1:61903/1‑36 (MQ=255)
     cAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACa                              <  1:4100425/36‑1 (MQ=255)
     cAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACa                              <  1:5216655/36‑1 (MQ=255)
      aaaaGCGCCCGGTAACGGCTGGTTGTTATTTTACaa                             >  1:4694210/1‑36 (MQ=21)
         aGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATc                          <  1:5991794/36‑1 (MQ=255)
          gcgcCCGGTAACGGCTGTTTTTTTTCTAACAAATcc                         >  1:1794847/1‑35 (MQ=25)
          gcgcCCGGTAACGGCTGTTTGTTTTCTTACCAATCg                         >  1:6007926/1‑36 (MQ=25)
           cgcCCGGTAACGGCTGTTTGTTATCTAACAAATCGt                        <  1:6367689/36‑1 (MQ=255)
            gcCCGGTAACGGCTGTTTTTTTTCTAACAAAACGTg                       >  1:6481737/1‑36 (MQ=25)
             cccGGTAACGGCTGTTTGTTATCTAACAAATCGTGc                      >  1:5160281/1‑36 (MQ=255)
              ccGGTAACGGCTGTTTTTTTTCTAACAAATCGTGCa                     >  1:5079850/1‑36 (MQ=37)
               cGGTAACGGCTGTTTGTTTTCTAACAAATCGCGCAt                    >  1:1534487/1‑36 (MQ=37)
               cGGTAACGGCTGTTTGTTATCTAACAAATCGTGCAt                    <  1:6613714/36‑1 (MQ=255)
               cGGTAACGGCTGTTTGTTATCTAACAAATCGTGCAt                    >  1:6903842/1‑36 (MQ=255)
               cGGGAACGGCTGTTTCTTCGTTTATAAATCGTGCAt                    >  1:2119130/1‑36 (MQ=2)
                ggTAACGGCTGTTTTTTTTCTAACAAATCGTGCAtt                   >  1:7241635/1‑36 (MQ=37)
                 gTAACGGCTGTTTGTTATCTAACAAATCGTGCATTc                  <  1:5204389/36‑1 (MQ=255)
                 gTAACGGCTGTTTGTTATCTAACAAATCGTGCATTc                  <  1:4557988/36‑1 (MQ=255)
                  taccGGCTGTTTGTTATCTAACAAATCGTGCATTcg                 <  1:1036690/33‑1 (MQ=255)
                  tAACGGCTGTTTGTTATCTAACAAATCGTGCATTcg                 <  1:884913/36‑1 (MQ=255)
                   aaCGGCTGTTTGTTATCTAACAAATCGTGCATTcgc                <  1:1942591/36‑1 (MQ=255)
                    aCGGCTGTTTGTTATCTAACAAATCGTGCATTcgcg               >  1:2703572/1‑36 (MQ=255)
                    aCGGCTGTTTGTTATCTAACAAATCGTGCATTcgcg               >  1:3337988/1‑36 (MQ=255)
                    aCGGCTGTTTGTTATCTAACAAATCGGGCATTcgcg               >  1:5885813/1‑36 (MQ=255)
                     cGGCTGTTTGTTATCTAACAAATCGTGCATTCGCgg              <  1:5310911/36‑1 (MQ=255)
                      ggCTGTTTGTTATCTAACAAATCGTGCATTCGCggg             >  1:6159831/1‑36 (MQ=255)
                        cTGTTTGTTATCTAACAAATCGTGCATTCGCGGGAt           <  1:7298612/36‑1 (MQ=255)
                        cTGTTTGTTATCTAACAAATCGTGCATTCGCGGGAt           >  1:2798254/1‑36 (MQ=255)
                         tgtttgttATCTAACAAATCGTGCATTCGCGGGATc          >  1:2628109/1‑36 (MQ=255)
                         tgtttgttATCTAACAAATCGTGCATTCGCGGGATc          >  1:4972034/1‑36 (MQ=255)
                         tgtttgttATCTAACAAATCGTGCATTCGCGGGATc          >  1:6746392/1‑36 (MQ=255)
                          gtttgttATCTAACAAATCGTGCATTCGCGGGATCa         >  1:3558692/1‑36 (MQ=255)
                          gtttgttATCTAACAAATCGTGCATTCGCGGGATCa         >  1:1512176/1‑36 (MQ=255)
                           tttgttATCTAACAAATCGTGCATTCGCGGGATCAg        >  1:3848358/1‑36 (MQ=255)
                            ttgttATCTAACAAATCGTGCATTCGCGGGATCAGc       >  1:597634/1‑36 (MQ=255)
                            ttgttATCTAACAAATCGTGCATTCGCGGGATCAGc       <  1:6944550/36‑1 (MQ=255)
                             tgttATCTAACAAATCGTGCATTCGCGGGATCAGCa      >  1:4488769/1‑36 (MQ=255)
                             tgttATCTAACAAATCGTGCATTCGCGGGATCAGCa      <  1:7040649/36‑1 (MQ=255)
                             tgttATCTAACAAATCGTGCATTCGCGGGATCAGCa      <  1:1971689/36‑1 (MQ=255)
                              gttATCTAACAAATCGTGCATTCGCGGGATCAGCAc     >  1:7386047/1‑36 (MQ=255)
                               ttATCTAACAAATCGTGCATTCGCGGGATCAGCAcc    >  1:44882/1‑36 (MQ=255)
                               ttATCTAACAAATCGTGCATTCGCGGGATCAGCAcc    >  1:3441547/1‑36 (MQ=255)
                                tATCTAACAAATCGTGCATTCGTGGGATCAGCAccc   <  1:4308801/36‑1 (MQ=255)
                                tATCTAACAAATCGTGCATTCGCGGGATCAGCAccc   >  1:2506338/1‑36 (MQ=255)
                                 aTCTAACAAATCGTGCATTCGTGGGATAAGCACCCa  >  1:847112/1‑36 (MQ=38)
                                 aTCTAACAAATCGTGCATTCGCGGGATCAGCACCCa  >  1:7230419/1‑36 (MQ=255)
                                 aTCTAACAAATCGTGCATTCGCGGGATCAGCACCCa  >  1:4701293/1‑36 (MQ=255)
                                 aTCTAACAAATCGTGCATTCGCGGGATCAGCACCCa  >  1:5342480/1‑36 (MQ=255)
                                 |                                   
AACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCA  >  REL606/2812818‑2812886

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG