Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,738,210 G→A 10.5% V172V (GTC→GTT tyrA ← fused chorismate mutase T and prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,738,2100GA10.5% 153.0 / 6.9 57V172V (GTC→GTTtyrAfused chorismate mutase T and prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (32/18);  new base A (3/3);  total (36/21)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACA  >  W3110S.gb/2738142‑2738279
                                                                    |                                                                     
cccAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGAc                                                                     >  1:897477/1‑71 (MQ=255)
 ccAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGAc                                                                     >  1:3434410/1‑70 (MQ=255)
   aGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAg                                                                  <  1:927560/71‑1 (MQ=255)
     cACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGaa                                                                >  1:852907/1‑71 (MQ=255)
         ggACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAc                                                             >  1:515292/1‑70 (MQ=255)
         ggACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACa                                                            <  1:2045227/71‑1 (MQ=255)
          gACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACaa                                                           <  1:1981646/71‑1 (MQ=255)
               tcatGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCttt                                                      >  1:1879397/1‑71 (MQ=255)
               tcatGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCtt                                                       >  1:2004257/1‑70 (MQ=255)
                  tGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCgg                                                   >  1:1592606/1‑71 (MQ=255)
                    cgcCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTa                                                 <  1:508584/71‑1 (MQ=255)
                      ccaccaGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTaaa                                               <  1:354400/71‑1 (MQ=255)
                      ccaccaGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTaaa                                               <  1:897839/71‑1 (MQ=255)
                          caGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCg                                           <  1:993141/71‑1 (MQ=255)
                           aGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTaaa                                               >  1:3481027/1‑66 (MQ=255)
                            gCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGt                                         <  1:923688/71‑1 (MQ=255)
                              aTGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTaa                                       <  1:2814175/71‑1 (MQ=255)
                              aTGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTaa                                       >  1:265418/1‑71 (MQ=255)
                               tGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCg                                                                       <  1:2040759/38‑1 (MQ=255)
                                 gCCTGTAATGGTCCATTTTTCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTaa                                       <  1:1344080/68‑1 (MQ=255)
                                 gCCTGTAATGGTCCATTTTTCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTaa                                       <  1:418385/68‑1 (MQ=255)
                                  ccTGTAATGGTCCATTTTTCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTaa                                       <  1:2123854/67‑1 (MQ=255)
                                      tAATGGCCCATTTTTCACTGATGCCAGATCGACCAGaa                                                                >  1:2252388/1‑38 (MQ=255)
                                      tAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTa                                                 <  1:276967/53‑1 (MQ=255)
                                      tAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCaa                               >  1:2295422/1‑71 (MQ=255)
                                      tAATGGCCCATTTTTAACTGATGCCAGATCGACCAGaa                                                                >  1:521677/1‑38 (MQ=255)
                                        aTGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATa                             <  1:1707769/71‑1 (MQ=255)
                                               aTTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCt                      >  1:1807054/1‑71 (MQ=255)
                                                tttttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGt                                         <  1:2186047/51‑1 (MQ=255)
                                                tttttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCt                      >  1:647371/1‑70 (MQ=255)
                                                tttttCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCGATAACTTGCt                      >  1:1897191/1‑70 (MQ=255)
                                                tttttCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCGATAACTTGCt                      >  1:1698111/1‑70 (MQ=255)
                                                 ttttCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCGATAACTTGCt                      >  1:3671202/1‑69 (MQ=255)
                                                  tttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAAt                                      >  1:1028129/1‑52 (MQ=255)
                                                  tttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAAt                                      >  1:327439/1‑52 (MQ=255)
                                                  tttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGc                                  >  1:2319650/1‑56 (MQ=255)
                                                  tttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATaa                            >  1:1119042/1‑62 (MQ=255)
                                                  tttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAg                   >  1:539052/1‑71 (MQ=255)
                                                  tttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAg                   >  1:775886/1‑71 (MQ=255)
                                                  tttCACTGATGCCAGATCGACCAGAATACAATCTTTAGGTAAAGGCGGTAAtt                                     >  1:1717312/1‑53 (MQ=255)
                                                   ttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCaa                               >  1:3842340/1‑58 (MQ=255)
                                                       cTGATGCCAGATCGACCAGAATACAATCTTTCGGTa                                                 <  1:2714307/36‑1 (MQ=255)
                                                       cTGATGCCAGATCGACCAGAATACAATCTTTCGGTa                                                 <  1:468168/36‑1 (MQ=255)
                                                        tGATGCCAGATCGACCAGAATACAATCTTTCGGTaa                                                >  1:3351226/1‑36 (MQ=255)
                                                        tGATGCCAGATCGACCAGAATACAATCTTTCGGTaa                                                >  1:2220064/1‑36 (MQ=255)
                                                          aTGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAAt                                      >  1:841864/1‑44 (MQ=255)
                                                          aTGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTgg           <  1:2563785/71‑1 (MQ=255)
                                                          aTGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTg            <  1:2639930/70‑1 (MQ=255)
                                                           tGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTg            <  1:2544157/69‑1 (MQ=255)
                                                           tGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGa          >  1:2667351/1‑71 (MQ=255)
                                                             ccAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGAtt        >  1:241234/1‑71 (MQ=255)
                                                                gATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCt                      <  1:2642858/54‑1 (MQ=255)
                                                                 aTCTACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGc     >  1:251242/1‑70 (MQ=255)
                                                                 aTCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGca    >  1:2885045/1‑71 (MQ=255)
                                                                 aTCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGca    >  1:1052114/1‑71 (MQ=255)
                                                                  tCGACCAGAATACAATCTTTCGGTAAAGGCGGTAAt                                      >  1:805555/1‑36 (MQ=255)
                                                                   cGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGcaca  >  1:2722932/1‑71 (MQ=255)
                                                                   cGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGcaca  >  1:2084596/1‑71 (MQ=255)
                                                                    |                                                                     
CCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACA  >  W3110S.gb/2738142‑2738279

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: