Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,229,323 T→G 4.6% G111G (GGT→GGG ycgN → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,229,3230T→G 4.6% 2.8 ‑1.4 ‑0.665333ycgNhypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (31/31):  total (34/31)
Fisher's exact test strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.62e-01

AAAGATTTACCTGCGTGGCATCCGCTACTTACGGGTTCGAAAGCGGCAATGCATGGTGAACGTATCTCTGT  >  REL606/1229288‑1229358
                                   |                                   
aaaGATTTACCTGCGTGGCATCCGCTACTTACGGGt                                     <  1:7233699/36‑1 (MQ=255)
 aaGATTTACCTGCGTGGCATCCGCTACTTACGGGtt                                    >  1:1163616/1‑36 (MQ=255)
 aaGATTTACCTGCGTGGCATCCGCTACTTACGGGtt                                    >  1:4085708/1‑36 (MQ=255)
 aaGATTTACCTGCGTGGCATCCGCTACTTACGGGtt                                    <  1:7413564/36‑1 (MQ=255)
 aaGATTTACCTGCGTGGCATCCGCTACTTACGGGtt                                    <  1:6725122/36‑1 (MQ=255)
  aGATTTACCTGCGTGGCATCCGCTACTTACGGTTTc                                   <  1:6072401/36‑1 (MQ=255)
   gATTTACCTGCGTGGCATCCGCTACTTACGGGTTCg                                  <  1:1485868/36‑1 (MQ=255)
   gATTTACCTGCGTGGCATCCGCTACTTACGGGTTCg                                  >  1:7562615/1‑36 (MQ=255)
    aTTTACCTGCGTGGCATCCGCTACTTACGGGTTCGa                                 >  1:3630133/1‑36 (MQ=255)
    aTTTACCTGCGTGGCATCCGCTACTTACGGGTTCGa                                 <  1:4609127/36‑1 (MQ=255)
     tttACCTGCGTGGCATCCGCTACTTACGGGTTCGaa                                >  1:6423565/1‑36 (MQ=255)
     tttACCTGCGTGGCATCCGCTACTTACGGGTTCGaa                                <  1:2233020/36‑1 (MQ=255)
      ttACCTGCGTGGCATCCGCTACTTACGGGTTCGaaa                               >  1:6246781/1‑36 (MQ=255)
      ttACCTGCGTGGCATCCGCTACTTACGGGTTCGaaa                               >  1:2141783/1‑36 (MQ=255)
      ttACCTGCGTGGCATCCGCTACTTACGGGTTCGaaa                               <  1:396134/36‑1 (MQ=255)
      ttACCTGCGTGGCATCCGCTACTTACGGGGTCGaaa                               >  1:4273689/1‑36 (MQ=255)
       tACCTGCGTGGCATCCGCTACTTACGGGTTCGAAAg                              <  1:2764587/36‑1 (MQ=255)
        aCCTGCGTGGCATCCGCTACTTACGGGTTCGAAAGc                             >  1:3049689/1‑36 (MQ=255)
        aCCTGCGTGGCATCCGCTACTTACGGGTTCGAAAGc                             >  1:249911/1‑36 (MQ=255)
        aCCTGCGTGGCATCCGCTACTTACGGGGTCGAAAGc                             >  1:3202769/1‑36 (MQ=255)
          cTGCGTGGCATCCGCTACTTACGGGTTCGAAAGCgg                           >  1:5834205/1‑36 (MQ=255)
           tGCGTGGCATCCGCTACTTACGGGTTCGAAAGCGGc                          <  1:2945010/36‑1 (MQ=255)
           tGCGTGGCATCCGCTACTTACGGGTTCGAAAGCGGc                          <  1:7688964/36‑1 (MQ=255)
           tGCGTGGCATCCGCTACTTACGGGTTCGAAAGCGGc                          <  1:120646/36‑1 (MQ=255)
            gCGTGGCATCCGCTACTTACGGGTTCGAAAGCGGCa                         >  1:5897005/1‑36 (MQ=255)
            gCGTGGCATCCGCTACTTACGGGTTCGAAAGCGGCa                         <  1:7389720/36‑1 (MQ=255)
            gCGTGGCATCCGCTACTTACGGGTTCGAAAGCGGCa                         >  1:859064/1‑36 (MQ=255)
             cGTGGCATCCGCTACTTACGGGTTCGAAAGCGGCaa                        >  1:6623251/1‑36 (MQ=255)
              gTGGCATCCGCTACTTACGGGTTCGAAAGCGGCAAt                       <  1:2004001/36‑1 (MQ=255)
              gTGGCATCCGCTACTTACGGGTTCGAAAGCGGCAAt                       >  1:7679914/1‑36 (MQ=255)
              gTGGCATCCGCTACTTACGGGTTCGAAAGCGGCAAt                       <  1:7656078/36‑1 (MQ=255)
              gTGGCATCCGCTACTTACGGGTTCGAAAGCGGCAAt                       >  1:1529069/1‑36 (MQ=255)
               tGGCATCCGCTACTTACGGGTTCGAAAGCGGCAATg                      >  1:4292623/1‑36 (MQ=255)
               tGGCATCCGCTACTTACGGGTTCGAAAGCGGCAATg                      <  1:4459083/36‑1 (MQ=255)
               tGGCATCCGCTACTTACGGGTTCGAAAGCGGCAATg                      <  1:2164490/36‑1 (MQ=255)
               tGGCATCCGCTACTTACGGGTTCGAAAGCGGCAATg                      >  1:4018726/1‑36 (MQ=255)
                ggCATCCGCTACTTACGGGTTCGAAAGCGGCAATGc                     <  1:8011084/36‑1 (MQ=255)
                ggCATCCGCTACTTACGGGTTCGAAAGCGGCAATGc                     <  1:5642502/36‑1 (MQ=255)
                   aTCCGCTACTTACGGGTTCGAAAGCGGCAATGCATg                  >  1:7263483/1‑36 (MQ=255)
                   aTCCGCTACTTACGGGTTCGAAAGCGGCAATGCATg                  <  1:237909/36‑1 (MQ=255)
                   aTCCGCTACTTACGGGTTCGAAAGCGGCAATGCATg                  >  1:3717744/1‑36 (MQ=255)
                    tCCGCTACTTACGGGTTCGAAAGCGGCAATGCATgg                 >  1:1976405/1‑36 (MQ=255)
                     ccGCTACTTACGGGTTCGAAAGCGGCAATGCATGGt                <  1:7064516/36‑1 (MQ=255)
                     ccGCTACTTACGGGTTCGAAAGCGGCAATGCATGGt                <  1:1735344/36‑1 (MQ=255)
                     ccGCTACTTACGGGTTCGAAAGCGGCAATGCATGGt                >  1:3602706/1‑36 (MQ=255)
                     ccGCGACTTACGGGTGCGAAAGCGGCAATGCATGGt                >  1:2419803/1‑36 (MQ=37)
                       tcTACTTACGGGTTCGAAAGCGGCAATGCATGGTGa              <  1:1899176/35‑1 (MQ=255)
                       gCTACTTACGGGTTCGAAATCGGCAATGCATGGTGa              <  1:7048689/36‑1 (MQ=255)
                       gCTACTTACGGGTTCGAAAGCGGCAATGCATGGTGa              >  1:7354123/1‑36 (MQ=255)
                       gCTACTTACGGGTTCGAAAGCGGCAATGCATGGTGa              <  1:3750889/36‑1 (MQ=255)
                       gCTACTTACGGGTTCGAAAGCGGCAATGCATGGTGa              <  1:3342779/36‑1 (MQ=255)
                        cTACTTACGGGTTCGAAAGCGGCAATGCATGGTGaa             <  1:4005139/36‑1 (MQ=255)
                        cTACTTACGGGTTCGAAAGCGGCAATGCATGGTGaa             >  1:2832932/1‑36 (MQ=255)
                         tACTTCCGGGTTCGAAAGCGGCAATGCATGGTGAAc            <  1:1690529/36‑1 (MQ=255)
                         tACTTACGGGTTCGAAAGCGGCAATGCATGGTGAAc            >  1:5394523/1‑36 (MQ=255)
                         tACTTACGGGTTCGAAAGCGGCAATGCATGGTGAAc            <  1:3873925/36‑1 (MQ=255)
                         tACTTACGGGTTCGAAAGCGGCAATGCATGGTGAAc            >  1:8370475/1‑36 (MQ=255)
                         tACTTACGGGGTCGAAAGCGGCAATGCATGGGGaag            >  1:3148075/1‑35 (MQ=37)
                          aCTTACGGGTTCGAACGCGGCAATGCATGGTGAACg           <  1:2764315/36‑1 (MQ=255)
                          aCTTACGGGTTCGAAAGCGGCAATGCATGGTGAACg           >  1:1250240/1‑36 (MQ=255)
                            ttACGGGTTCGAAAGGGGCAATGCATGGTTAACGTa         <  1:2041660/36‑1 (MQ=38)
                            ttACGGGTTCGAAAGCGGCAATGCATGGTGAACGTa         >  1:1548586/1‑36 (MQ=255)
                            ttACGGGTTCGAAAGCGGCAATGCATGGTGAACGTa         >  1:2475792/1‑36 (MQ=255)
                            ttACGGGTTCGAAAGCGGCAATGCATGGTGAACGTa         <  1:1113292/36‑1 (MQ=255)
                             tACGGGTTCGAAAGCGGCAATGCATGGTGAGCGTAt        <  1:3209815/36‑1 (MQ=255)
                             tACGGGTTCGAAAGCGGCAATGCATGGTGAACGTAt        >  1:6680951/1‑36 (MQ=255)
                             tACGGGTTCGAAAGCGGCAATGCATGGTGAACGTAt        >  1:5375417/1‑36 (MQ=255)
                                gggTTCGAAAGCGGCAATGCATGGTGAACGTAtctc     >  1:4225640/1‑36 (MQ=255)
                                gggTTCGAAAGCGGCAATGCATGGTGAACGTAtctc     <  1:3711548/36‑1 (MQ=255)
                                  gTTCGAAAGCGGCAATGCATGGTGAACGTATCTCtg   >  1:443126/1‑36 (MQ=255)
                                   ttCGAAAGCGGCAATGCATGGTGAACGTATCTCtgt  >  1:6687271/1‑36 (MQ=255)
                                   ttCGAAAGCGGCAATGCATGGTGAACGTATCTCtgt  <  1:1509983/36‑1 (MQ=255)
                                   |                                   
AAAGATTTACCTGCGTGGCATCCGCTACTTACGGGTTCGAAAGCGGCAATGCATGGTGAACGTATCTCTGT  >  REL606/1229288‑1229358

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG