Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,771,424 C→T P332S (CCG→TCG)  pheA → fused chorismate mutase P and prephenate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,771,4240CT100.0% 275.2 / NA 88P332S (CCG→TCG) pheAfused chorismate mutase P and prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (50/38);  total (50/38)

GCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCCCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAATCAGCG  >  minE/1771358‑1771489
                                                                  |                                                                 
gcgcTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGAt                                                               >  1:2756243/1‑71 (MQ=255)
gcgcTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGAt                                                               >  1:2058698/1‑71 (MQ=255)
  gcTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTc                                                             <  1:1848163/71‑1 (MQ=255)
   cTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTc                                                             <  1:313975/70‑1 (MQ=255)
    tGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCAc                                                           >  1:1550651/1‑71 (MQ=255)
     ggTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACg                                                          >  1:1912218/1‑71 (MQ=255)
      gTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCAc                                                           <  1:4077166/69‑1 (MQ=255)
      gTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCAc                                                           <  1:4000405/69‑1 (MQ=255)
       ttGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCAc                                                           <  1:3068698/68‑1 (MQ=255)
       ttGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGt                                                        <  1:3876351/71‑1 (MQ=255)
        tGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGt                                                        >  1:1199744/1‑70 (MQ=255)
         gAATCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTaa                                                      >  1:3898575/1‑71 (MQ=255)
         gAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTaa                                                      >  1:51968/1‑71 (MQ=255)
         gAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTa                                                       >  1:652825/1‑70 (MQ=255)
          aaGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACgg                                                         >  1:440827/1‑67 (MQ=255)
          aaGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGGCTGGAATCACGCTCGATTCACGGTAAt                                                     <  1:1102556/71‑1 (MQ=255)
             cGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCa                                                  <  1:1010085/71‑1 (MQ=255)
              gTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCAt                                                 <  1:265978/71‑1 (MQ=255)
               ttGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGa                                                                >  1:3486373/1‑55 (MQ=255)
                tGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCAt                                                 >  1:2180894/1‑69 (MQ=255)
                tGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATgg                                               >  1:507748/1‑71 (MQ=255)
                 gCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCAt                                                 >  1:64432/1‑68 (MQ=255)
                  cTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGa                                             <  1:636161/71‑1 (MQ=255)
                   tGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCAt                                                 <  1:1239808/66‑1 (MQ=255)
                   tGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCAt                                                 <  1:2302904/66‑1 (MQ=255)
                      tACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAagag                                         <  1:603241/71‑1 (MQ=255)
                      tACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAagag                                         <  1:2767971/71‑1 (MQ=255)
                       aCTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGaa                                            >  1:2509608/1‑67 (MQ=255)
                       aCTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGaa                                            >  1:2606917/1‑67 (MQ=255)
                          gcgcAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGAtt                                                              >  1:1944305/1‑46 (MQ=255)
                          gcgcAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAaga                                          >  1:2861974/1‑66 (MQ=255)
                          gcgcAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAaga                                          >  1:1809916/1‑66 (MQ=255)
                           cgcAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTaa                                                      >  1:3150689/1‑53 (MQ=255)
                           cgcAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTa                                                       <  1:677097/52‑1 (MQ=255)
                           cgcAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATg                                      <  1:748206/69‑1 (MQ=255)
                             cAACCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCt                                  <  1:2673921/71‑1 (MQ=255)
                              aaCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCAc                                                           <  1:1176827/45‑1 (MQ=255)
                              aaCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTa                                 <  1:405813/71‑1 (MQ=255)
                              aaCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTa                                 <  1:3544411/71‑1 (MQ=255)
                              aaCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTa                                 <  1:2941808/71‑1 (MQ=255)
                               aCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTa                                                       <  1:783040/48‑1 (MQ=255)
                               aCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTa                                 >  1:1127714/1‑70 (MQ=255)
                               aCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTa                                 >  1:1231758/1‑70 (MQ=255)
                               aCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTAt                                >  1:3093791/1‑71 (MQ=255)
                               aCCACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTAt                                >  1:1448619/1‑71 (MQ=255)
                                ccACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATc                               >  1:3437759/1‑71 (MQ=255)
                                ccACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATc                               >  1:3469261/1‑71 (MQ=255)
                                ccACAATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATc                               >  1:3534763/1‑71 (MQ=255)
                                 cacaATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCt                              <  1:590989/71‑1 (MQ=255)
                                 cacaATCTGATTATGACCAGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCt                              >  1:2964114/1‑71 (MQ=255)
                                   caATCTGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTg                             >  1:2056527/1‑70 (MQ=255)
                                     aTCTGATTATGACCCGTCTGGAATCACGCTCGATTCa                                                            >  1:3590358/1‑37 (MQ=255)
                                      tCTGATTATGACCCGTCTGGAATCACGCTCGATTCa                                                            >  1:864204/1‑36 (MQ=255)
                                        tGATTATGACCCGTCTGGAATCACGCTCGATTCACTGTAATCCATGGGAAGAGATGTTCTa                                 <  1:356902/61‑1 (MQ=255)
                                        tGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTa                                 <  1:2520869/61‑1 (MQ=255)
                                        tGATTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGatat                        >  1:266684/1‑70 (MQ=255)
                                          aTTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGatat                        <  1:356041/68‑1 (MQ=255)
                                          aTTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGatat                        <  1:3105451/68‑1 (MQ=255)
                                          aTTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGata                         >  1:203740/1‑67 (MQ=255)
                                          aTTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGata                         >  1:3747170/1‑67 (MQ=255)
                                          aTTATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCa                     <  1:3863038/71‑1 (MQ=255)
                                           ttATGACCCGTCTGGAATCCCGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAg                    >  1:2951799/1‑71 (MQ=255)
                                           ttATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAg                    <  1:1431887/71‑1 (MQ=255)
                                           ttATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAg                    >  1:1361737/1‑71 (MQ=255)
                                            tATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAag                                           >  1:2899194/1‑47 (MQ=255)
                                            tATGACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAgg                   >  1:2624258/1‑71 (MQ=255)
                                               gACCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGcc                 >  1:195999/1‑70 (MQ=255)
                                                aCCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCaa               >  1:1210578/1‑71 (MQ=255)
                                                aCCCGTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCa                <  1:590857/70‑1 (MQ=255)
                                                    gTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCtt           >  1:1739180/1‑71 (MQ=255)
                                                    gTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCtt           <  1:1109942/71‑1 (MQ=255)
                                                    gTCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCtt           >  1:349351/1‑71 (MQ=255)
                                                     tCTGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGa                           <  1:3541442/54‑1 (MQ=255)
                                                       tGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGata                         >  1:3490160/1‑54 (MQ=255)
                                                       tGGAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGaa        >  1:716514/1‑71 (MQ=255)
                                                         gAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTa                                 <  1:1354636/44‑1 (MQ=255)
                                                         gAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTa                                 <  1:1532331/44‑1 (MQ=255)
                                                         gAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAAt       >  1:987164/1‑70 (MQ=255)
                                                         gAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAAt       >  1:458758/1‑70 (MQ=255)
                                                         gAATCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAAt       >  1:3116059/1‑70 (MQ=255)
                                                          aaTCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGa         <  1:850583/67‑1 (MQ=255)
                                                          aaTCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAATc      <  1:733788/70‑1 (MQ=255)
                                                          aaTCACGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAATCa     <  1:619851/71‑1 (MQ=255)
                                                              aCGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTc                      <  1:3102445/50‑1 (MQ=255)
                                                              aCGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGa         >  1:1774643/1‑63 (MQ=255)
                                                              aCGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAATCAGc   >  1:1242808/1‑69 (MQ=255)
                                                              aCGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAATCAGCg  >  1:1861765/1‑70 (MQ=255)
                                                              aCGCTCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAATCAGCg  >  1:1803884/1‑70 (MQ=255)
                                                                  |                                                                 
GCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATGACCCGTCTGGAATCACGCCCGATTCACGGTAATCCATGGGAAGAGATGTTCTATCTGGATATTCAGGCCAATCTTGAATCAGCG  >  minE/1771358‑1771489

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: ‑