Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 523,473 A→G 16.1% E290G (GAA→GGA)  arcC → predicted carbamate kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606523,4730A→G16.1% 2.0 ‑4.4 0.039869arcCpredicted carbamate kinase
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (15/17):  total (22/17)
Fisher's exact test strand distribution p-value = 1.24e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGC  >  REL606/523438‑523508
                                   |                                   
ggCGTTATCGTGAATTGAAGAGACGCTGGCGGGCGa                                     <  1:6421005/36‑1 (MQ=255)
ggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGa                                     >  1:6084402/1‑36 (MQ=255)
 gCGTTNTCGCGAATTGAAGAGACGCTGGCGGGCGaa                                    <  1:6112200/36‑1 (MQ=255)
 gCGTTATCGCGAATTGAAGAGACGCTGGCGGGGGaa                                    >  1:7261482/1‑36 (MQ=255)
   gTTATCGCGAATTGAAGAGACGCTGGCGGGCGGAgg                                  >  1:6549052/1‑35 (MQ=255)
   gTTATCGCGAATTGAAGAGACGCTGGCGGGCGGAgg                                  >  1:1563098/1‑35 (MQ=255)
   gTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGc                                  >  1:3145773/1‑36 (MQ=255)
   gTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGc                                  <  1:2164905/36‑1 (MQ=255)
    ttATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCg                                 >  1:3528152/1‑36 (MQ=255)
    ttATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCg                                 <  1:2901597/36‑1 (MQ=255)
     tATCGCGAATTGAAGAGATGCTGGCGGGGGGAGCgg                                >  1:4736253/1‑36 (MQ=21)
      aTCGCGAATTGAAGAGACGCTGGCGGGCGGAGCggg                               >  1:2210764/1‑36 (MQ=255)
      aTCGCGAATTGAAGAGACGCTGGCGGGCGAAGCggg                               <  1:5624819/36‑1 (MQ=255)
       tCGCGAATTGAAGAGACGCTGGCGGGCGAAGCgggg                              <  1:1816741/36‑1 (MQ=255)
         gcgAATTGAAGAGACGCTGGCGGGGGGAGCgggggg                            >  1:1208187/1‑34 (MQ=25)
         gcgAATTGAAGAGACGCTGGCGGGCGAAGCGGGGAc                            <  1:7173336/36‑1 (MQ=255)
           gAATTGAAGAGACGCTGTCGGGCGAAGCGGGGACCt                          <  1:7378840/36‑1 (MQ=255)
           gAATTGAAGAGACGCTGGCGGGCGAAGCGGGGGCCt                          >  1:3553048/1‑36 (MQ=255)
            aaTTGAAGAGACGCTGGCGGGCGAAGCGGGGGcccg                         >  1:7225434/1‑34 (MQ=37)
            aaTTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTg                         <  1:196061/36‑1 (MQ=255)
            aaTTGAAGAGACGCTGGCGGGCGAAACGGGGACGTg                         >  1:6371594/1‑36 (MQ=37)
             aTTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGt                        <  1:627332/36‑1 (MQ=255)
              ttGAAGAGACGCTGGCGGGCGGAGCGGGGGCCTGTa                       >  1:3755593/1‑36 (MQ=37)
              ttGAAGAGACGCTGGCGGGCGAAGCGGGGGCCGGtt                       >  1:664768/1‑35 (MQ=25)
              ttGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTa                       <  1:7511681/36‑1 (MQ=255)
               tGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTAt                      <  1:2124384/36‑1 (MQ=255)
                gAAGAGACGCTGTCGGGCGAAGCGGGGACCTGTAtt                     <  1:6617081/36‑1 (MQ=255)
                 aaGAGACGATGGCGGGCGGGGCGGGGACCTGGAttt                    >  1:6849810/1‑36 (MQ=19)
                  agagACGCTGGCGGGCGAAGCGGGGACCTGTATTTc                   <  1:5937555/36‑1 (MQ=255)
                    agaCGCTGGCGGGCGAAACGGGGACCTGTATTTCGc                 >  1:6067710/1‑36 (MQ=255)
                         cTGGCGGGCGAATCGGGGTCCTGTATTTCGCTGTAg            <  1:3037870/36‑1 (MQ=37)
                         cTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAg            <  1:4013316/36‑1 (MQ=255)
                         cTGGCGGGCGAAGCGGGGACATGTATTTCGCTGTAg            >  1:475927/1‑36 (MQ=255)
                          tGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGt           >  1:5048363/1‑36 (MQ=255)
                          tGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGt           >  1:1964682/1‑36 (MQ=255)
                             cgggcgAAGCGGGGACCTGTATTTCGCTGTAGtgat        >  1:4833235/1‑36 (MQ=255)
                             cgggcgAAGCGGGGACCTGTATTTCGCTGTAGtgat        <  1:1353112/36‑1 (MQ=255)
                               ggcgAAGCGGGGCCCTGTATTTCGCTGTAGtgatga      >  1:112350/1‑36 (MQ=255)
                               ggcgAAGCGGGGACCTGTATTTCGCTGTAGtgatga      <  1:664120/36‑1 (MQ=255)
                                 cgAAGCGGGGACCTGTATTTCGCTGTAGTGATgaga    >  1:65038/1‑36 (MQ=255)
                                 cgAAGCGGGGACCTGTATTTCGCTGTAGTGATgaga    >  1:4576279/1‑36 (MQ=255)
                                 cgAAGCGGGGACCTGTATTTCGCTGTAGTGATgaga    <  1:3277039/36‑1 (MQ=255)
                                  gAAGCGGGGACCTGTATTTCGCTGTAGTGATGagag   >  1:1631608/1‑36 (MQ=255)
                                   aaGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGc  >  1:3258485/1‑36 (MQ=255)
                                   |                                   
GGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTGATGAGAGC  >  REL606/523438‑523508

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