Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 976,493 C→T 11.9% S197F (TCT→TTT)  aroA → 3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606976,4930C→T11.9% 2.1 ‑1.9 ‑0.042590aroA3‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (22/13):  total (28/14)
Fisher's exact test strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.35e-01

AGATACGGTGATTCGTATTAAAGGCGATCT.GGTTTCTAAACCTTATATCGACATCACACTCAATCTGAT  >  REL606/976458‑976526
                                    |                                 
aGATACGGTGATTCGTATTAAAGGCGATCT.GGTTTc                                   >  1:4423515/1‑36 (MQ=255)
 gATACGGTGATTCGTATTAAAGGCGATCT.GGTTTCt                                  >  1:5697228/1‑36 (MQ=255)
  aTACGGCGATTCGTATTAAAGGCGATCT.GGTTTCTa                                 >  1:2697156/1‑36 (MQ=255)
    aCGGTGATTCGTATTAAAGGCGATCT.GGTTTCTaaa                               <  1:1344915/36‑1 (MQ=255)
      ggTGATTCGTATTAAAGGCGATTTGGGTTTTTaaat                              >  1:121593/1‑35 (MQ=16)
      ggTGATTCGTATTAAAGGCGATCT.TGTTTCTAAAcc                             >  1:4949228/1‑36 (MQ=255)
      ggTGATTCGTATTAAAGGCGATCT.GGTTTCTAAAcc                             >  1:165185/1‑36 (MQ=255)
      ggTGATTCGTATTAAAGGCGATCT.GGTTTCTAAAcc                             >  1:569765/1‑36 (MQ=255)
       gTGATTCGTATTAAAGGCGATCT.GGTTTCTAAACCt                            >  1:3668661/1‑36 (MQ=255)
        tGATTCGTATTAAAGGCGATCT.GGTTTTTAAACCtt                           >  1:2521342/1‑36 (MQ=255)
         gATTCGTATTAAAGGCGATTT.GTTTTCTAAACCTta                          <  1:1295666/36‑1 (MQ=37)
          aTTCGTATTAAAGGCGATCT.GGTTTTTAAACCTtat                         >  1:3508699/1‑36 (MQ=255)
           ttCGTATTAAAGGCGATCT.GGTTTTTAAACCTtata                        >  1:4675299/1‑36 (MQ=255)
           ttCGTATTAAAGGCGATCG.GGTTTTTAAACCTtata                        >  1:2666004/1‑36 (MQ=38)
            tCGTATTAAAGGCGATCT.GGTTTCTAAACCTTatat                       >  1:1795862/1‑36 (MQ=255)
             cGTATTAAAGGCGATCT.GGTTTCTAACCCTTATATc                      <  1:5399746/36‑1 (MQ=255)
             cGTATTAAAGGCGATCT.GGTTTCTAAACCTTATATc                      >  1:5837720/1‑36 (MQ=255)
             cGTATTAAAGGCGATCT.GGTTTCTAAACCTTATATc                      >  1:2608094/1‑36 (MQ=255)
             cGTATTAAAGGCGATCT.GGTTTCTAAACCTTATATc                      >  1:1541802/1‑36 (MQ=255)
              gTATTAAAGGCGATCT.GGTTTCTAAACCTTATATCg                     >  1:892943/1‑36 (MQ=255)
              gTATTAAAGGCGATCT.GGTTTCTAAACCTTATATCg                     >  1:2108669/1‑36 (MQ=255)
                aTTAAAGGCGATCT.GGTTTCTAAACCTTATATCGAc                   <  1:1313065/36‑1 (MQ=255)
                aTTAAAGGCGATCT.GGGTTCTAAACCTTATATCGAc                   >  1:3787909/1‑36 (MQ=255)
                 ttAAAGGCGATCT.GGTTTTTAAACCTTATATCGACa                  >  1:4664333/1‑36 (MQ=255)
                 ttAAAGGCGATCT.GGTTTCTAAACCTTATATCGACa                  <  1:3947202/36‑1 (MQ=255)
                  tAAAGGCGATCT.GGTTTCTAAACCTTATATCGACAt                 >  1:3523843/1‑36 (MQ=255)
                   aaaGGTGATCT.GGTTTCTAAACCTTATATCGACATc                <  1:5668228/36‑1 (MQ=255)
                   aaaGGCGATCT.GGTTTCTAAACCTTATATCGACATc                >  1:225476/1‑36 (MQ=255)
                    aaGGCGATCT.GGTTTCTAAACCTTATATCGACATca               <  1:4126055/36‑1 (MQ=255)
                     aGGCGATCT.GGTTTCTAAACCTTATATCGCCATcac              <  1:5662726/36‑1 (MQ=255)
                     aGGCGATCT.GGTTTCTAAACCTTATATCGACATcac              >  1:2315588/1‑36 (MQ=255)
                     aGGCGATCC.GGTTTCTAAACCTTATATCGACATcac              >  1:4239461/1‑36 (MQ=255)
                      ggCGGTCN.NNTTTGTTAGCCTGATATCGACATcaca             <  1:2028189/36‑1 (MQ=9)
                       gCGATCT.GGTTTCTAAACCTTATATCGACATCacac            >  1:3795334/1‑36 (MQ=255)
                         gATCT.GGTTTCTAAACCTTATATCGACATCACACTc          <  1:5445613/36‑1 (MQ=255)
                         gATCT.GGTTTCTAAACCTTATATCGACATCACACTc          >  1:291881/1‑36 (MQ=255)
                          aTCT.GGTTTCTAAACCTTATATCGACATCACACTCa         <  1:5144641/36‑1 (MQ=255)
                            cT.GGTTTCTAAACCTTATATCGACATCACACTCAAt       <  1:2563348/36‑1 (MQ=255)
                             t.GGTTTCTAAACCTTATATCGACATCACACTCAATc      <  1:4751803/36‑1 (MQ=255)
                             t.GGTTTCTAAACCTTATATCGACATCACACTCAATc      >  1:563709/1‑36 (MQ=255)
                               ggTTTCTAAACCTTATATCGACATCACACTCAATCt     <  1:5610882/36‑1 (MQ=255)
                                 tttCTAAACCTTATATCGACATCACACTCAATCtga   >  1:3765355/1‑36 (MQ=255)
                                  ttCTAAACCTTATATCGAAATCACACTCAATCtgat  >  1:2470848/1‑36 (MQ=255)
                                    |                                 
AGATACGGTGATTCGTATTAAAGGCGATCT.GGTTTCTAAACCTTATATCGACATCACACTCAATCTGAT  >  REL606/976458‑976526

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG