Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,099,912 A→G 19.4% K380E (AAA→GAA)  yegQ → predicted peptidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,099,9120A→G19.4% 3.3 ‑2.9 ‑0.0301138yegQpredicted peptidase
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (12/12):  total (18/12)
Fisher's exact test strand distribution p-value = 5.68e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.83e-01

GTGAGCGCAAGGGGGACCTCGCGGCGG.TAGCGGTGAAAAATAAATTCTCCGTTGGCGACAGCCTTGAGCTG  >  REL606/2099877‑2099947
                                    |                                   
gTGAGCGCAAGGGGGACCTCGCGGCGG.TAGCGGTGa                                     <  1:466610/36‑1 (MQ=255)
  gAGCGCAAGGGGGACCTCGCGGCGG.TAACGGGGaaa                                   >  1:412141/1‑36 (MQ=37)
   aGCGCAAGGGGGACCTCGCGGCGG.TAGCGGTGaaaa                                  <  1:5253516/36‑1 (MQ=255)
     cgcAAGGGGGACCTCGCGGCGG.TAGCGGTGAAaaat                                <  1:2140285/36‑1 (MQ=255)
      gcAAGGGGGACCTCGCGGCGG.GAGGGGTGGAgaata                               >  1:858866/1‑36 (MQ=19)
        aaGGGGGACCTCGCGGCGG.TCGCGGGGGAaaataaa                             >  1:415206/1‑36 (MQ=37)
         aGGGGGACCTCGCGGCGG.GGGCGGGGGAaaataaaa                            >  1:3103592/1‑35 (MQ=19)
         aGGGGGACCTCGCGGCGG.GAGCGGTGGAaaataaat                            >  1:3739832/1‑36 (MQ=37)
           ggggACCTCGCGGCGGCTGGGGGGGAAaaaaaaatc                           >  1:3975721/1‑35 (MQ=9)
            gggACCTCGCGGCGG.TAGCGGTGAAAAATAAATtct                         <  1:3288123/36‑1 (MQ=255)
            gggACCTCGCGGCGG.TAGCGGTGAAAAATAAATtct                         <  1:3694509/36‑1 (MQ=255)
             ggACCTCGCGGCGG.TAGCGGGGGAAAAAAAAAtctc                        >  1:1857502/1‑36 (MQ=19)
              gACCTCGCGGCGG.TAGCGGTGAAAAATAAATTCTcc                       <  1:3728908/36‑1 (MQ=255)
               gccTCGAGTCGG.TAGCGGTGAAAAATAAATTCTCCg                      <  1:4152412/35‑1 (MQ=37)
                ccTCGCGGCGG.TTATGGCGGAAAAAAAAATCTCCGt                     >  1:4747251/1‑36 (MQ=0)
                ccTCGCGGCGG.TAGCGGTGAAAAATAAATTCTCCGt                     <  1:2620841/36‑1 (MQ=255)
                ccTCGCGGCGG.TAGCGGTGAAAAATAAATTCTCCGt                     <  1:5501338/36‑1 (MQ=255)
                 cTCGCGGCGG.TAGCGGTGAAAAATAAATTCTCCGtt                    <  1:3057410/36‑1 (MQ=255)
                  tCGCGGCAG.TAGCGGTGAAAAATAAATTCTCCGTTg                   <  1:2357419/36‑1 (MQ=255)
                  tCGCGGAGG.TAGCGGTGAAAAATAAATTCTCCGTTg                   <  1:3280749/36‑1 (MQ=255)
                   cgcgGCGG.TAGCGGTGAAAAATAAATTCTCCGTTgg                  <  1:4031696/36‑1 (MQ=255)
                     cggcgg.TATCGGTGAAAAATAAATTCTCCGTTGGCg                <  1:4392127/36‑1 (MQ=255)
                         gg.TAGCGGTGAAAAATAAATTCTCCGTTGGCGACAg            >  1:1649064/1‑36 (MQ=255)
                         gg.TAGCGGTGAAAAATAAATTCTCCGTTGGCGACAg            >  1:3233/1‑36 (MQ=255)
                           gcaGCGGTGAAAAATAAATTCTCCGCTGGCGACAGc           >  1:2889853/3‑36 (MQ=38)
                                ggTGAAAAATAAATTCTCCGTTGGCGACAGCCTTGa      >  1:604242/1‑36 (MQ=255)
                                ggTGAAAAATAAATTCTCCGTCGGCGACAGCCTTGa      >  1:5052742/1‑36 (MQ=255)
                                 gcgAAAAATAAATTCTCCGTTGGCGACAGCCTTGAg     >  1:3802312/3‑36 (MQ=255)
                                 gTGAAAAATAAATTCTCCGTTGGCGACAGCCTTGAg     >  1:2422544/1‑36 (MQ=255)
                                 gTGAAAAATAAATTCTCCGTTGGCGACAGCCTTGAg     >  1:851059/1‑36 (MQ=255)
                                  tGAAAAATAAATTCTCCGTTGGCGACAGCCTTGAGc    >  1:1419596/1‑36 (MQ=255)
                                   gAAAAATAAATTCTCCGTTGGCGACAGCCTTGAGCt   >  1:3545587/1‑36 (MQ=255)
                                   gAAAAATAAATTCTCCGTTGGCGACAGCCTTGAGCt   >  1:2979931/1‑36 (MQ=255)
                                    aaaaaTAAATTCTCCGTCGGCGACAGCCTTGAGCtg  >  1:5922700/1‑36 (MQ=255)
                                    |                                   
GTGAGCGCAAGGGGGACCTCGCGGCGG.TAGCGGTGAAAAATAAATTCTCCGTTGGCGACAGCCTTGAGCTG  >  REL606/2099877‑2099947

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