Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 313,785 C→A 13.6% S4R (AGC→AGA yahJ → deaminase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606313,7850C→A13.6% 3.0 ‑3.1 ‑0.03612yahJdeaminase
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (18/13):  total (18/18)
Fisher's exact test strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.23e-01

TTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAA  >  REL606/313751‑313820
                                  |                                   
ttattaAAACGTAAGGACGTAAGATGAAAGAAAGCa                                    <  1:6518884/36‑1 (MQ=255)
ttattaAAACGTAAGGACGTAAGATGAAAGAAAGCa                                    >  1:8370618/1‑36 (MQ=255)
 tattaAAACGTAAGGACGTAAGATGAAAGAAAGCaa                                   >  1:7954558/1‑36 (MQ=255)
  attaAAACGTAAGGACGTAAGATGAAAGAAAGCAAt                                  >  1:2695606/1‑36 (MQ=255)
  attaAAACGTAAGGACGTAAGATGAAAGAAAGCAAt                                  <  1:2957051/36‑1 (MQ=255)
   ttaAAGCGTAAGGACGTAATATTAAAGAAAGCAATa                                 <  1:4033246/36‑1 (MQ=37)
   ttaAAACGTAAGGACGTAAGATGAAAGAAAGCAATa                                 >  1:5120497/1‑36 (MQ=255)
   ttaAAACGTAAGGACGTAAGATGAAAGAAAGCAATa                                 >  1:3504936/1‑36 (MQ=255)
    taAAACGTAAGGACGTAAGATGAAAGAAAGCAATAg                                <  1:7211388/36‑1 (MQ=255)
    taAAACGTAAGGACGTAAGATGAAAGAAAGCAATAg                                >  1:4799205/1‑36 (MQ=255)
    taAAACGTAAGGACGTAAGATGAAAGAAAGCAATAg                                <  1:2486947/36‑1 (MQ=255)
      aaaCGTAAGGACGTAAGATGAAAGAAAGCAATAgcc                              <  1:1765141/36‑1 (MQ=255)
      aaaCGTAAGGACGTAAGATGAAAGAAAGCAATAgcc                              <  1:4173288/36‑1 (MQ=255)
       aaCGTAAGGACGTAAGATGAAAGAAAGCAATAgccg                             >  1:6813167/1‑36 (MQ=255)
         cGTAAGGACGTAAGATGAAAGAAAGCAATAgccgcc                           <  1:6069032/36‑1 (MQ=255)
          gTAAGGACGTAAGATGAAAGAAAGCAATAGccgccg                          >  1:365795/1‑36 (MQ=255)
          gTAAGGACGTAAGATGAAAGAAAGCAATAGccgccg                          >  1:6203551/1‑36 (MQ=255)
            aaGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTg                        >  1:2530769/1‑36 (MQ=255)
            aaGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTg                        >  1:1421338/1‑36 (MQ=255)
            aaGGACGTAAAATGAAAGAAAGAAATAGCCGCCGTg                        <  1:2232628/36‑1 (MQ=38)
              ggACGTAAGATGAAAGAAAGCAATAGCCGCCGTGaa                      <  1:4465092/36‑1 (MQ=255)
               gACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAAt                     >  1:2041340/1‑36 (MQ=255)
                  gaaaGATGAAAGAAAGCAATAGCCGACGTGAATTTc                  >  1:4138735/3‑36 (MQ=38)
                   tAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCt                 >  1:2824305/1‑36 (MQ=255)
                    aaGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTg                >  1:7819480/1‑36 (MQ=255)
                    aaGATGAAAAAAAACAATAGCCGCCGTGAATTTCTg                <  1:1033813/36‑1 (MQ=37)
                     aGATAAAAGAAAGCAATAGCCGCCGTGAATTTCTGa               <  1:1215095/36‑1 (MQ=255)
                     aGATAAAAAAAAGAAATAGCCGCCGTGAATTTCTGa               <  1:5297102/36‑1 (MQ=25)
                      gATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAg              >  1:7182060/1‑36 (MQ=255)
                       ataaaagaaaaCAATAGCCGCCGTGAATTTCTGAGc             <  1:8399256/33‑1 (MQ=37)
                       ataaaaaaaagAAATAGCCGCCGTGAATTTCTGAGc             <  1:2513742/33‑1 (MQ=25)
                       aaaaaagaaagAAATAGCCGCCGTGAATTTCTGAGc             <  1:5001084/33‑1 (MQ=37)
                       aTGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGc             >  1:3812653/1‑36 (MQ=255)
                           aaaaaaaCAATAGCCGCCGTGAATTTCTGAGCCaga         <  1:4955526/33‑1 (MQ=37)
                             gaaagCAATAGCCGCCGTGAATTTCTGAGCCAGAGc       >  1:7688893/1‑36 (MQ=255)
                               aagAAATACCCGCCGTGAATTTCTGAGCCAGAGCgg     <  1:4459232/36‑1 (MQ=37)
                                  cAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTaa  >  1:7214932/1‑36 (MQ=255)
                                  cAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTaa  >  1:7812054/1‑36 (MQ=255)
                                  cAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTaa  >  1:1816116/1‑36 (MQ=255)
                                  |                                   
TTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAA  >  REL606/313751‑313820

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