Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 33,798 +G 6.3% coding (77/1149 nt) carA → carbamoyl‑phosphate synthase small subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL60633,7981.→G 6.3% 4.1 ‑2.6 ‑0.532coding (77/1149 nt)carAcarbamoyl‑phosphate synthase small subunit
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (7/23):  total (9/23)
Fisher's exact test strand distribution p-value = 7.26e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.00e-01

TTTCACGGTCGGGCCATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTCGTTTTCAATACTTCAATGACCG  >  REL606/33764‑33833
                                   |                                   
tttCACGGTCGGGCCATAGGGGCAACAGGTTcggc.g                                    <  1:5455001/36‑1 (MQ=255)
tttCACGGTCGGGCCATAGGGGCAACAGGTTcggc.g                                    >  1:5116191/1‑36 (MQ=255)
tttCACGGTCGGGCCATAGGGGCAACAGGTTcggc.g                                    >  1:6762067/1‑36 (MQ=255)
 ttCACGGTCGGGCCATAGGGGCAACAGGTTcggc.gg                                   >  1:7471824/1‑36 (MQ=255)
  tCACGGTCGGGCCATAGGGGCAACAGGTTcggc.ggg                                  >  1:266758/1‑35 (MQ=255)
  tCACGGTCGGGCCATAGGGGCAACAGGTTCGGC.GGt                                  <  1:5472760/36‑1 (MQ=255)
  tCACGGTCGGGCCATAGGGGCAACAGGTTCGGC.GGt                                  <  1:157270/36‑1 (MQ=255)
      ggTCGGGCCATAGGGGCAACAGGTTCGGC.GGTTggg                              <  1:6034701/36‑1 (MQ=255)
       gacGGGCCATAGGGGCAACAGGTTCGGC.GGTTgggg                             <  1:5882232/34‑1 (MQ=255)
       gTCGGGCCATAGGGGCAACAGGTTCGGC.GGTTgggg                             <  1:804128/36‑1 (MQ=255)
         cGGGCCATAGGGGCAACAGGTTCGGC.GGTTGGGGaa                           <  1:796681/36‑1 (MQ=255)
          gggCCATAGGGGCAACAGGTTCGGCGGGTTgggggg                           >  1:3290531/1‑34 (MQ=37)
          gggCCATAGGGGCAACAGGTTCGGC.GGGTggggggg                          >  1:4275957/1‑33 (MQ=38)
           ggCCATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGt                         <  1:4413329/36‑1 (MQ=255)
           ggCCATAGGGGCAACAGGTTCGGC.GGGTGGGGGAGt                         >  1:6239029/1‑36 (MQ=38)
            gCCATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTc                        <  1:7880220/36‑1 (MQ=255)
            gCCATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTc                        <  1:97140/36‑1 (MQ=255)
            gCCATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTc                        <  1:4433249/36‑1 (MQ=255)
             ccATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTCg                       <  1:5545742/36‑1 (MQ=255)
             ccATAGGGGCAACAGGTTCGGC.GGGTGGGGGAGGCg                       >  1:6685883/1‑36 (MQ=37)
              cATAGGGGCAACAGGTTCGGCGGGTTGGGGAGGtgt                       >  1:2904965/1‑34 (MQ=25)
              cATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTCGt                      <  1:5555427/36‑1 (MQ=255)
              cATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTCGt                      <  1:5817390/36‑1 (MQ=255)
               aTAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTCGtt                     <  1:7789367/36‑1 (MQ=255)
                  ggggCAACAGGTTCGGC.GGTTGGGGAAGTCGTTTTc                  <  1:2253913/36‑1 (MQ=255)
                  ggggCAACAGGTTCGGC.GGTTGGGGAAGTCGTTTTc                  <  1:7900548/36‑1 (MQ=255)
                   gggCAACAGGTTCGGC.GGTTTGGGGAGGCGgttttt                 >  1:5585123/1‑34 (MQ=16)
                    ggCAACAGGTTCGGC.GGTTGGGGAAGTCGTTTTCaa                <  1:5360109/36‑1 (MQ=255)
                     gCAACAGGTTCGGC.GGTTGGGGAAGTCGTTTTCAAt               <  1:4613460/36‑1 (MQ=255)
                     gCAACAGGTTCGGC.GGTTGGGGAAGTCGTTTTCAAt               <  1:3194976/36‑1 (MQ=255)
                       aaCAGGTTCGGC.GGTTTGGGGAGTCGTTTTCCATag             >  1:7097183/1‑35 (MQ=25)
                       aaCAGGTTCGGC.GGTTGGGGAAGTCGTTTTCAATAc             <  1:6399142/36‑1 (MQ=255)
                       aaCAGGTTCGGC.GGTTGGGGAAGTCGTTTTCAATAc             <  1:4929549/36‑1 (MQ=255)
                         cAGGTTCGGC.GGTTGGGGAAGTCGTTTTCAATACtt           <  1:3583195/36‑1 (MQ=255)
                           ggTTCGGC.GGTTGGGGGAGGGGTTTTTTATACttta         >  1:6717059/1‑34 (MQ=9)
                            gTTCGGC.GGTTGGGGGAGTCGTTTTTAATGCTTCaa        >  1:6165037/1‑36 (MQ=25)
                              tCGGC.GGTTGGGGAAGTCGTTTTCAATACTTCAATg      <  1:944819/36‑1 (MQ=255)
                                ggc.ggTTGGGGAAGTCGTTTTCAATACTTCAATTAc    >  1:1765408/1‑36 (MQ=255)
                                ggc.ggTTGGGGAAGTCGTTTTCAATACTTCAATGAc    <  1:3605062/36‑1 (MQ=255)
                                ggc.ggTTGGGGAAGTCGTTTTCAATACTTCAATGAc    <  1:92902/36‑1 (MQ=255)
                                  c.ggTTGGGGAAGTCGTTTTCAATACTTCAATGACCg  >  1:6828871/1‑36 (MQ=255)
                                  c.ggTTGGGGAAGTCGTTTTCAATACTTCAATGACCg  <  1:3075716/36‑1 (MQ=255)
                                  c.ggTTGGGGAAGTCGTTTTCAATACTTCAATGACCg  >  1:5121958/1‑36 (MQ=255)
                                   |                                   
TTTCACGGTCGGGCCATAGGGGCAACAGGTTCGGC.GGTTGGGGAAGTCGTTTTCAATACTTCAATGACCG  >  REL606/33764‑33833

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG