Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 523,481 A→G 16.2% T293A (ACC→GCC)  arcC → predicted carbamate kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606523,4810A→G16.2% 7.2 ‑3.1 ‑0.053877arcCpredicted carbamate kinase
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (27/17):  total (36/17)
Fisher's exact test strand distribution p-value = 4.42e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01

CGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCATCCGGCA  >  REL606/523446‑523516
                                   |                                   
cgcgAATTGAAGAGACGCTGGCGGGCGAAGCGGGGa                                     <  1:7632284/36‑1 (MQ=255)
 gcgAATTGAAGAGACGCTGGCGGGCGAAGCGGGGAc                                    <  1:5483248/36‑1 (MQ=255)
   gAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCt                                  <  1:3653592/36‑1 (MQ=255)
   gAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCt                                  <  1:5365619/36‑1 (MQ=255)
     aTTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGt                                <  1:3824386/36‑1 (MQ=255)
      ttGAAGAGACGCTGGCGGGGGGAGCGGGGGCCTGtg                               >  1:5346412/1‑35 (MQ=21)
      ttGAAGAGACGCTGGCGGGCGAAGCGGGGGGCTGtg                               >  1:6085966/1‑35 (MQ=25)
       tGAAGAGACGCTGGCGGGCGAAGCGGGGGGCTGtgg                              >  1:2741591/1‑34 (MQ=25)
       tGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTAt                              <  1:7406939/36‑1 (MQ=255)
        gAAGAGACGCTGGCGGGCGGAAAGGGGGGCTGTAtt                             >  1:2953359/1‑36 (MQ=12)
         aaGAGACGCTGGCGGGCGAAACGGGGGCCTGTAttt                            >  1:7346812/1‑36 (MQ=37)
          agagACGCTGGCGGGCGCAGCGGGGGCATGTCtttt                           >  1:6944666/1‑35 (MQ=16)
          agagAAGCTGGCGGGAGAAGGGGGGAACTGTAtttt                           >  1:7395331/1‑35 (MQ=16)
           gagaCGCTTGCGGGCGAAGCGGGGACCTGTATTTCg                          <  1:7558609/36‑1 (MQ=255)
            agaCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGc                         <  1:5584509/36‑1 (MQ=255)
             gaCGCTGGCGGGCGAAGCGGGGGGCTGTATTTCGCt                        >  1:178563/1‑36 (MQ=37)
             gaCGCTGGCGGGCGAAGCGGGGGCCTGTATTTCggt                        >  1:2489521/1‑34 (MQ=38)
             gaCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCt                        <  1:4123757/36‑1 (MQ=255)
             gaCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCt                        <  1:7897241/36‑1 (MQ=255)
                 cTGGCGTGCGAAGCGGGGACCTGTATTTCGCTGtcg                    >  1:3536813/1‑34 (MQ=38)
                 cTGGCGGGCGAAGCGGGGATCTGTATTTCGCTGttg                    >  1:139660/1‑34 (MQ=38)
                 cTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAg                    >  1:3697839/1‑36 (MQ=255)
                 cTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAg                    <  1:7860560/36‑1 (MQ=255)
                  tGGCGGGCGGAGCGGGGGACTGTTTTTTGCTTTAGt                   >  1:2186773/1‑36 (MQ=9)
                  tGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGt                   <  1:3875114/36‑1 (MQ=255)
                  tGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGt                   <  1:5443076/36‑1 (MQ=255)
                   ggcgggcgAAGCGGGGACCTGTATTTCGCTGTAGtg                  <  1:4855287/36‑1 (MQ=255)
                   ggcgggcgAAGCGGGGACCTGTATTTCGCTGTAGtg                  >  1:6633482/1‑36 (MQ=255)
                    gcgggcgAAGCGGGGACCTGTATTTCGCTGTAGtga                 >  1:3643547/1‑36 (MQ=255)
                    gcgggcgAAGCGGGGACCTGTATTTCGCTGTAGtga                 >  1:5392972/1‑36 (MQ=255)
                     cgggcgAAGCGGGGACCTGTATTTCGCTGTAGtgat                >  1:1801235/1‑36 (MQ=255)
                      gggcgAAGCGGGGACCTGTATTTCGCTGTAGtgatg               <  1:6668535/36‑1 (MQ=255)
                       ggcgAAGCGGGGACCTGTATTGCGCTGTAGtgatga              <  1:2182395/36‑1 (MQ=255)
                        gcgAAGCGGGGACCTGTATTTCGCTGTAGtgatggg             >  1:6385518/1‑34 (MQ=255)
                         cgAAGCGGGGACCTGTATTTCGCTGTAGTGATgaga            >  1:5769312/1‑36 (MQ=255)
                          gAAGCGGGGACCTGTGTTTCGCTGTAGTGATGagag           >  1:7203515/1‑36 (MQ=255)
                          gAAGCGGGGACCTGTATTTCGCTGTAGTGATGagag           >  1:7746257/1‑36 (MQ=255)
                          gAAGCGGGGACCTGTATTTCGCTGTAGTGATGagag           >  1:7202978/1‑36 (MQ=255)
                          gAAGCGGGAACCTGTATATAGCTGTAGTGATGagag           >  1:6729063/1‑36 (MQ=37)
                           aaGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGc          >  1:1984201/1‑36 (MQ=255)
                            aGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCa         >  1:1207146/1‑36 (MQ=255)
                            aGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCa         >  1:4960233/1‑36 (MQ=255)
                             gCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCAt        >  1:2801336/1‑36 (MQ=255)
                              cGGGGACCTGTATTTCGCTGTAGTGATGAGAGCATc       <  1:2803196/36‑1 (MQ=255)
                              cGGGGACCTGTATTTCGCTGTAGTGATGAGAGCATc       >  1:645286/1‑36 (MQ=255)
                               ggggACCTGTATTTCGCTGTAGTGATGAGAGCATcc      >  1:1533584/1‑36 (MQ=255)
                               ggggACCTGTATTTCGCTGTAGTGATGAGAGCATcc      >  1:4724986/1‑36 (MQ=255)
                                gggACCTGTATTTCGCTGTAGTGATGAGAGCATCCg     <  1:7006969/36‑1 (MQ=255)
                                gggACCTGTATTTCGCTGTAGTGATGAGAGCATCCg     >  1:1918717/1‑36 (MQ=255)
                                gggACCTGTATTTCGCTGTAGTGATGAGAGCATCCg     >  1:2710933/1‑36 (MQ=255)
                                gggACCTGTATTTCGCTGTAGTGATGAGAGCATCCg     >  1:7587098/1‑36 (MQ=255)
                                gggACCTGTATTTATCTGGAGTGATGTGAGCATGCg     >  1:5896674/1‑36 (MQ=16)
                                 ggACCTGTATTTCGCTGTAGGGATGAGAGCATCCgg    >  1:3504323/1‑36 (MQ=255)
                                  gACCTGTATTTCGCTGTAGTGATGAGAGCATCCGGc   >  1:383983/1‑36 (MQ=255)
                                   aCCGGTATTTCGCTGGAGTGATGAGAGCATCCGGCa  >  1:6130434/1‑36 (MQ=37)
                                   |                                   
CGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCATCCGGCA  >  REL606/523446‑523516

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG