Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,738,198 T→C 11.3% S176S (TCA→TCG 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,1980TC11.3% 132.2 / 8.1 53S176S (TCA→TCGtyrAfused chorismate mutase T and prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (27/20);  new base C (3/3);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAG  >  W3110S.gb/2738133‑2738262
                                                                 |                                                                
gggTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGc                                                             <  1:1781256/71‑1 (MQ=255)
       gCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATc                                                       <  1:2267241/70‑1 (MQ=255)
       gCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATc                                                       >  1:1207094/1‑70 (MQ=255)
       gCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATc                                                       >  1:1310543/1‑70 (MQ=255)
         cccAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGAc                                                    >  1:897477/1‑71 (MQ=255)
          ccAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGAc                                                    >  1:3434410/1‑70 (MQ=255)
            aGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAg                                                 <  1:927560/71‑1 (MQ=255)
             gCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGa                                                                >  1:1545369/1‑55 (MQ=255)
             gCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGa                                                                >  1:3782638/1‑55 (MQ=255)
             gCACCGGACCATAATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGa                                                                >  1:420584/1‑55 (MQ=255)
              cACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACt                                                                  <  1:1022055/52‑1 (MQ=255)
              cACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACt                                                                  <  1:2568449/52‑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)
                      catcatGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGAt                                                               <  1:1848155/47‑1 (MQ=255)
                      catcatGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGAt                                                               <  1:3361721/47‑1 (MQ=255)
                      catcatGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGAt                                                               <  1:1160580/47‑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:265418/1‑71 (MQ=255)
                                       aTGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTaa                      <  1:2814175/71‑1 (MQ=255)
                                        tGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCg                                                      <  1:2040759/38‑1 (MQ=255)
                                          gCCTGTAATGGTCCATTTTTCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTaa                      <  1:418385/68‑1 (MQ=255)
                                          gCCTGTAATGGTCCATTTTTCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTaa                      <  1:1344080/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:1698111/1‑70 (MQ=255)
                                                         tttttCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCGATAACTTGCt     >  1:1897191/1‑70 (MQ=255)
                                                          ttttCACCGATGCCAGATCAACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCGATAACTTGCt     >  1:3671202/1‑69 (MQ=255)
                                                           tttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAAt                     >  1:327439/1‑52 (MQ=255)
                                                           tttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAAt                     >  1:1028129/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:468168/36‑1 (MQ=255)
                                                                cTGATGCCAGATCGACCAGAATACAATCTTTCGGTa                                <  1:2714307/36‑1 (MQ=255)
                                                                 tGATGCCAGATCGACCAGAATACAATCTTTCGGTaa                               >  1:2220064/1‑36 (MQ=255)
                                                                 tGATGCCAGATCGACCAGAATACAATCTTTCGGTaa                               >  1:3351226/1‑36 (MQ=255)
                                                                 |                                                                
GGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAG  >  W3110S.gb/2738133‑2738262

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