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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606523,4810A→G14.1% 2.1 ‑1.0 0.036877arcCpredicted carbamate kinase
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (18/12):  total (23/13)
Fisher's exact test strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCATCCGGCA  >  REL606/523450‑523516
                               |                                   
aaTTGAAGAGACGCTGGCGGGCGGAACGGGGACGTg                                 >  1:5727718/1‑36 (MQ=21)
   tGAAGAGACGCTGGCGGGCGGAGCGGGGGGCTGtgt                              >  1:2458192/1‑34 (MQ=18)
     aaGAGACGCTGGAGGGCCAAGGGGGGGGGTGTAttt                            >  1:5058497/1‑36 (MQ=9)
       gagaCGCTGGCGGGGGGAGCGGGGGCCGGTGttttg                          >  1:1753433/1‑34 (MQ=12)
       gagaCGCTGGCGGGGGAAGCGGGGGCCGGGATTTCg                          >  1:307483/1‑36 (MQ=16)
         gaCGCTGGGGGGCGAAGCGGGGACCTGTATTTCGCt                        <  1:502616/36‑1 (MQ=255)
          aCGCTGGCGGGCGAAGCGGGGGCCTGTATTTCGCTg                       >  1:2912044/1‑36 (MQ=255)
           cGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGt                      <  1:3689080/36‑1 (MQ=255)
              tGGCGGGCGAAGCGGGGAACCGTATTTAGCTGTTGt                   >  1:420421/1‑36 (MQ=16)
                gcgggcgAAGCGGGGACCTGTATTTCGCTGTAGtga                 <  1:5882695/36‑1 (MQ=255)
                gcgggcgAAGCGGGGACCTGTATTTCGCTGTAGtga                 >  1:4410636/1‑36 (MQ=255)
                gcgggcgAAGCGGGGACCTGTATTTCGCTGTAGtga                 <  1:4076086/36‑1 (MQ=255)
                  gggcgAAGCGGGGACCTGTATTTCGCTGTAGtgatg               >  1:4848637/1‑36 (MQ=255)
                  gggcgAAGCGGGGACCTGTATTTCGCTGTAGtgatg               >  1:3617569/1‑36 (MQ=255)
                     cgAAGCGGGGACCTGTATTTCGCTGTAGTGATgaga            <  1:2867347/36‑1 (MQ=255)
                     cgAAGCGGGGACCTGTATTTCGCTGTAGTGATgaga            >  1:5202782/1‑36 (MQ=255)
                     cgAAGCGGGGACCTGTATTTCGCTGTAGTGATgaga            <  1:475555/36‑1 (MQ=255)
                      gAAGCGGGGACCTGTATTTCGCTGTAGTGATGagag           >  1:5680669/1‑36 (MQ=255)
                      gAAGCGGGGACCTGTATTTCGCTGTAGTGATGagag           >  1:785965/1‑36 (MQ=255)
                      gAAGCGGGGACCTGTATTTCGCTGTAGTGATGagag           >  1:3217238/1‑36 (MQ=255)
                      gAAGCGGGGACCTGTATTTCGATGTAGTGATGagag           >  1:1939025/1‑36 (MQ=255)
                       aaGCGGGGGCCTGTATTTCGCTGTAGTGATGAGAGc          <  1:5777379/36‑1 (MQ=255)
                       aaGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGc          <  1:499248/36‑1 (MQ=255)
                       aaGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGc          <  1:5136468/36‑1 (MQ=255)
                       aaGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGc          <  1:5235158/36‑1 (MQ=255)
                        aGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCa         <  1:1535339/36‑1 (MQ=255)
                         gCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCAt        <  1:3169493/36‑1 (MQ=255)
                          cGGGGACCTGTATTTTGCTGTAGTGGTGAGAGAATc       >  1:2877879/1‑36 (MQ=25)
                           ggggACCTGTATTTCGCTGTAGTGATGAGAGCATcc      >  1:1130125/1‑36 (MQ=255)
                            gggACCTGTATTTCGCTGTAGTGATTAGAGCATcag     >  1:651913/1‑34 (MQ=38)
                            gggACCTGTATTTCGCTGTAGTGATGAGAGCATCCg     <  1:1474342/36‑1 (MQ=255)
                             ggACCTGTATTTCGCTGTAGTGATGAGAGCATCCgg    >  1:5476971/1‑36 (MQ=255)
                             ggACCTGTATTTCGCTGTAGTGATGAGAGCATCCgg    >  1:2183661/1‑36 (MQ=255)
                             ggACCTGTATTTCGCTGTAGTGATGAGAGCATCCgg    >  1:3231709/1‑36 (MQ=255)
                              gACCTGTATTTCGCTGTAGTGATGAGAGCATCCGGc   >  1:4576957/1‑36 (MQ=255)
                              gACCTGTATTTCGCTGTAGTGATGAGAGCATCCGGc   >  1:3809808/1‑36 (MQ=255)
                               aCCTGTATTTCGCTGTAGTGATGAGAGCATCCGGCa  >  1:5337692/1‑36 (MQ=255)
                               |                                   
AATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGCATCCGGCA  >  REL606/523450‑523516

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