Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,440,342 +G 5.4% coding (1308/1374 nt) mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,440,3421.→G 5.4% 4.0 ‑2.7 ‑0.137coding (1308/1374 nt)mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (8/27):  total (10/27)
Fisher's exact test strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.79e-01

CAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC.GGTTTTGGTGGGATCCATCAGAAACTGCTGGATGG  >  REL606/4440307‑4440377
                                    |                                   
cAGCCTGGCGACCATATTCTGGTGATGAGCAAcggc                                      <  1:4525394/36‑1 (MQ=255)
cAGCCTGGCGACCATATTCTGGTGATGAGCAAcggc                                      >  1:653693/1‑36 (MQ=255)
 aGCCTGGCGACCATATTCTGGTGATGAGCAAcggc.g                                    <  1:3229390/36‑1 (MQ=255)
 aGCCTGGCGACCATATTCTGGTGATGAGCAAcggc.g                                    <  1:2429385/36‑1 (MQ=255)
 aGCCTGGCGACCATATTCTGGTGATGAGCAAcggc.g                                    >  1:6626210/1‑36 (MQ=255)
  gCCTGGCGACCATATTCTGGTGATGAGCAAcggc.gg                                   >  1:6492328/1‑36 (MQ=255)
  gCCTGGCGACCATATTCTGGTGATGAGCAAcggc.gg                                   >  1:774547/1‑36 (MQ=255)
   ccTGGCGACCATATTCTGGTGATGAGCAACGGC.GGt                                  <  1:2063630/36‑1 (MQ=255)
    cTGGCGACCATATTCTGGTGATGAGCAACGGC.GGtt                                 <  1:7180724/36‑1 (MQ=255)
    cTGGCGACCATATTCTGGTGATGAGCAACGGC.GGtt                                 >  1:6281798/1‑36 (MQ=255)
    cTGGCGACCATATTCTGGTGATGAGCAACGGC.GGtt                                 <  1:303977/36‑1 (MQ=255)
     tGGCGACCATATTCTGGTGATGAGCAACGGC.GGttt                                <  1:4156664/36‑1 (MQ=255)
     tGGCGACCATATTCTGGTGATGAGCAACGGC.GGgtt                                >  1:5607983/1‑36 (MQ=255)
     tGGCGAACATATTCTGGTGATGAGCAACGGC.GGttt                                <  1:7706209/36‑1 (MQ=255)
      ggCGACCATATTCTGGTGATGAGCAACGGC.GGtttt                               <  1:4583478/36‑1 (MQ=255)
       gCGACCATATTCTGGTGATGAGCAACGGC.GGTTTtg                              <  1:283859/36‑1 (MQ=255)
        cGACCATATTCTGGTGATGAGCAACGGC.GGgttttg                             >  1:721052/1‑34 (MQ=37)
        cGACCATATTCTGGTGATGAGCAACGGC.GGTTTtgg                             <  1:1408789/36‑1 (MQ=255)
        cGACCATATTCTGGTGATGAGCAACGGC.GGTTTtgg                             <  1:7705179/36‑1 (MQ=255)
        cGACCATATTCTGGTGATGAGCAACGGC.GGTTTtgg                             <  1:2986212/36‑1 (MQ=255)
         gACCATATTCTGGTGATGAGCAACGGC.GGTTTtggt                            <  1:1665054/36‑1 (MQ=255)
         gACCATATTCTGGTGATGAGCAACGGC.GGGTTttgt                            >  1:103000/1‑36 (MQ=25)
          aacATATTCTGGTGATGAGCAACGGC.GGTTTttgtg                           >  1:5894191/3‑36 (MQ=37)
           ccATATTCTGGTGATGAGCAACGGC.GGTTTtggtgg                          <  1:2116424/36‑1 (MQ=255)
            cATATTCTGGTGATGAGCAACGGC.GGTTTTGGTggg                         <  1:1255753/36‑1 (MQ=255)
             atatTCTGGTGATGAGCAACGGCGGGTTTTGTTggg                         >  1:1055326/1‑36 (MQ=21)
             atatTCTGGTGATGAGCAACGGC.GGTTTTGTTGGGa                        <  1:1887090/36‑1 (MQ=255)
               atTCTGGTGATGAGCAACGGC.GGTTTTTTTgggggg                      >  1:7111667/1‑33 (MQ=21)
               atTCTGGTGATGAGCAACGGC.GGTTTTGGTGGGATc                      <  1:3117370/36‑1 (MQ=255)
               atTCGGGTGATGAGCAACGGC.GGGTTTTTTGGGGTc                      >  1:5690770/1‑36 (MQ=12)
                 tCTGGTGATGAGCAACGGC.GGTTTTGGTGGGATCCa                    <  1:1724540/36‑1 (MQ=255)
                 tCTGGTGATGAGCAACGGC.GGTTTTGGTGGGATCCa                    <  1:5781305/36‑1 (MQ=255)
                  cTGGTGATGAGCAACGGC.GGTTTTGGTGGGATCCAt                   <  1:3330239/36‑1 (MQ=255)
                   tGGTGATGAGCAACGGCGGGTTTTGTTGGGTTcctt                   >  1:5719002/1‑34 (MQ=12)
                   tGGTGATGAGCAACGGC.GGTTTTGGTGGGATCCATc                  <  1:5514048/36‑1 (MQ=255)
                   tGGTGATGAGCAACGGC.GGTTTTGGTGGGATCCATc                  <  1:4685505/36‑1 (MQ=255)
                    ggTGATGAGCAACGGC.GGTTTTGGTGGGATCCATCa                 <  1:5690195/36‑1 (MQ=255)
                    ggTGATGAGCAACGGC.GGTTTTGGTGGGATCCATCa                 <  1:2923576/36‑1 (MQ=255)
                    ggTGATGAGCAACGGC.GGTTTTGGTGGGATCCATCa                 <  1:725878/36‑1 (MQ=255)
                    ggTGATGAGCAACGGC.GGTTTTGGTGGGATCCATCa                 <  1:730517/36‑1 (MQ=255)
                       gatgaGCAACGGC.GGTTTTGGTGGGATCCATCAGaa              <  1:1165837/36‑1 (MQ=255)
                        atgaGCAACGGC.GGTTTTGGTGGGATCCATCAGaaa             <  1:5037935/36‑1 (MQ=255)
                          gaGCAACGGC.GGTTTTGGTGGGATCCATCAGAAAct           <  1:4816062/36‑1 (MQ=255)
                          gaGCAACGGC.GGGTTTTGTTGGGTTCATCAGaaaat           >  1:2655300/1‑34 (MQ=9)
                                cggc.ggTTTTGGTGGGATCCATCAGAAACTGCTGGa     <  1:2756361/36‑1 (MQ=255)
                                cggc.ggTTTTGGTGGGATCCATCAGAAACTGCTGGa     <  1:6837305/36‑1 (MQ=255)
                                 ggc.ggTTTTGGTGGGATCCATCAGAAACTGCTGGAt    <  1:4825546/36‑1 (MQ=255)
                                   c.ggTTTTGGTGGGATCCATCAGAAACTGCTGGATgg  >  1:5630508/1‑36 (MQ=255)
                                    |                                   
CAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC.GGTTTTGGTGGGATCCATCAGAAACTGCTGGATGG  >  REL606/4440307‑4440377

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG