Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,772,153 A→G 43.8% G201G (GGT→GGC tyrA ← fused chorismate mutase T and prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,772,1530AG43.8% 25.1 / 83.5 105G201G (GGT→GGCtyrAfused chorismate mutase T and prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (24/35);  new base G (23/23);  total (47/58)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 105 total aligned reads displayed.

ATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTA  >  minE/1772089‑1772211
                                                                |                                                          
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:1948344/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:211112/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:1577054/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:1147132/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:2647918/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:2618950/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:3291274/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:2347984/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:2431253/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:2742385/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:251690/1‑68 (MQ=255)
aTGCTTCCGGCTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCg                                                         >  1:1465653/1‑68 (MQ=255)
 tGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGaa                                                       >  1:656392/1‑69 (MQ=255)
 tGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGaa                                                       >  1:910954/1‑69 (MQ=255)
 tGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGaa                                                       >  1:3372704/1‑69 (MQ=255)
 tGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGaa                                                       >  1:1004708/1‑69 (MQ=255)
 tGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGaa                                                       >  1:2669944/1‑69 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACa                                                     >  1:3275065/1‑59 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCggg                                                >  1:764173/1‑64 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGgt                                               >  1:2806379/1‑65 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGgt                                               >  1:2232815/1‑65 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGgt                                               >  1:939235/1‑65 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGgt                                               >  1:3226993/1‑65 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGgt                                               >  1:2989606/1‑65 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGgt                                               >  1:824340/1‑65 (MQ=255)
             aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTa                                            >  1:2271311/1‑68 (MQ=255)
                tctaTCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGc                                          >  1:2077593/4‑67 (MQ=255)
                tctaTCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGc                                          >  1:1313898/4‑67 (MQ=255)
                tctaTCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGc                                          >  1:1998716/4‑67 (MQ=255)
                tctaTCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGc                                          >  1:2576615/4‑67 (MQ=255)
                tctaTCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGc                                          >  1:1460516/4‑67 (MQ=255)
                tctaTCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGc                                          >  1:3366322/4‑67 (MQ=255)
                tctaTCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGc                                          >  1:2319222/4‑67 (MQ=255)
                tctaTCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGc                                          >  1:532095/4‑67 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:1194806/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:924697/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:2794989/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:2459937/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:331230/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:2389688/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:2383374/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:2184917/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:1601804/1‑55 (MQ=255)
                      acacCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCggg                                                >  1:1917000/1‑55 (MQ=255)
                             aCCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCTGACca                           <  1:710080/69‑1 (MQ=255)
                             aCCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACca                           <  1:683580/69‑1 (MQ=255)
                             aCCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACca                           <  1:552695/69‑1 (MQ=255)
                             aCCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACca                           <  1:2688194/69‑1 (MQ=255)
                             aCCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACca                           <  1:200725/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCTAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:2750167/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:2958695/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:2283866/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:481967/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:489120/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:2499912/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:544210/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:1858361/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:1556602/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:1538974/69‑1 (MQ=255)
                              ccACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACcat                          <  1:1051633/69‑1 (MQ=255)
                                  aaCTTGCTTTGCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:1120963/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACGGGGGCATCGTgc                     <  1:2938567/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:513977/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:2660029/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:3316413/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:91963/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:3177792/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:3113727/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:1437926/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:3053328/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:2873956/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:2871743/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:2613945/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:2272174/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:2489912/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:2523424/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:2586997/69‑1 (MQ=255)
                                   aCTTGCTTTGCCAGGCTACCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:258946/69‑1 (MQ=255)
                                    cTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATg                      >  1:1550591/1‑67 (MQ=255)
                                     ttGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATgcgc                   <  1:3134176/69‑1 (MQ=255)
                                     ttGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATgcgc                   <  1:1237605/69‑1 (MQ=255)
                                     ttGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATgcgc                   <  1:1458122/69‑1 (MQ=255)
                                     ttGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATgcgc                   <  1:2186042/69‑1 (MQ=255)
                                     ttGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATgcgc                   <  1:340171/69‑1 (MQ=255)
                                     ttGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATgcgc                   <  1:339484/69‑1 (MQ=255)
                                     ttGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATgcgc                   <  1:467946/69‑1 (MQ=255)
                                       gCTTTGCCAGTCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACca                           <  1:407763/59‑1 (MQ=255)
                                                    aCCGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:698291/52‑1 (MQ=255)
                                                     ccGCTGTCCGGGCCGAACATCGGGTGTAACCCCAGCACCGGGCCATCGTgc                     <  1:3000300/51‑1 (MQ=255)
                                                       gCTTTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCATCATGGCCTGt   <  1:2797062/67‑1 (MQ=255)
                                                       gCTGTCCTGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCCCCACCAGCATGGCCTGTa  <  1:192885/68‑1 (MQ=255)
                                                       gCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGt   <  1:1584121/67‑1 (MQ=255)
                                                       gCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGt   <  1:1516406/67‑1 (MQ=255)
                                                       gCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGt   <  1:3010793/67‑1 (MQ=255)
                                                       gCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGt   <  1:893702/67‑1 (MQ=255)
                                                       gCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGt   <  1:3106528/67‑1 (MQ=255)
                                                       gCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGt   <  1:1248279/67‑1 (MQ=255)
                                                       gCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGt   <  1:3170568/67‑1 (MQ=255)
                                                       gCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGt   <  1:1037294/67‑1 (MQ=255)
                                                           tCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACca                           <  1:577076/39‑1 (MQ=255)
                                                                |                                                          
ATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTA  >  minE/1772089‑1772211

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: