Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 734,169 C→T 57.1% P268L (CCG→CTG)  yceG → predicted aminodeoxychorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE734,1690CT57.1% 20.2 / 65.7 56P268L (CCG→CTG) yceGpredicted aminodeoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base T (19/13);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGAT  >  minE/734099‑734236
                                                                      |                                                                   
gacccgaccGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGcc                                                                     <  1:690693/71‑1 (MQ=255)
 acccgaccGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCg                                                                    <  1:4587881/71‑1 (MQ=255)
    cgaccGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGa                                                                   >  1:2829949/1‑69 (MQ=255)
          tGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGtat                                                           <  1:1273130/71‑1 (MQ=255)
          tGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGtat                                                           <  1:4698550/71‑1 (MQ=255)
            aTTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATaa                                                         >  1:2664537/1‑71 (MQ=255)
               tACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAacac                                                      >  1:3505327/1‑71 (MQ=255)
               tACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAacac                                                      >  1:1554344/1‑71 (MQ=255)
               tACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAaca                                                       <  1:3444783/70‑1 (MQ=255)
                aCGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAacac                                                      <  1:3356899/70‑1 (MQ=255)
                 cGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACAcc                                                     <  1:4143050/70‑1 (MQ=255)
                 cGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACAcc                                                     <  1:4501887/70‑1 (MQ=255)
                    gATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCtata                                                 >  1:585861/1‑71 (MQ=255)
                      tGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATAcc                                               <  1:1040205/71‑1 (MQ=255)
                             gCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCg                                        >  1:908610/1‑71 (MQ=255)
                              cGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGa                                                                   <  1:4104734/43‑1 (MQ=255)
                              cGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCgg                                       <  1:4267289/71‑1 (MQ=255)
                               gTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGt                                      <  1:3383676/71‑1 (MQ=255)
                               gTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCgg                                       <  1:2715431/70‑1 (MQ=255)
                                ttATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCg                                                              <  1:2514866/46‑1 (MQ=255)
                                ttATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCta                                                   <  1:4022107/57‑1 (MQ=255)
                                 tataATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCAtt                                            <  1:1285307/63‑1 (MQ=255)
                                  ataATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTg                                   >  1:1547959/1‑71 (MQ=255)
                                   taATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCg                                                              <  1:3368786/43‑1 (MQ=255)
                                   taATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATaa                                                         >  1:2373207/1‑48 (MQ=255)
                                   taATGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTgc                                  >  1:445463/1‑71 (MQ=255)
                                    aaTGGCAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATaa                                                         >  1:2268715/1‑47 (MQ=255)
                                         cAAACTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAaca                                                       <  1:706478/44‑1 (MQ=255)
                                          aaaCTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAgg                           >  1:3315368/1‑71 (MQ=255)
                                          aaaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGtat                                                           <  1:4249573/39‑1 (MQ=255)
                                          aaaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATaa                                                         >  1:4035512/1‑41 (MQ=255)
                                          aaaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAacac                                                      <  1:2498731/44‑1 (MQ=255)
                                          aaaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAgg                           >  1:1176816/1‑71 (MQ=255)
                                           aaCTTTCTCGTGCAGACCTGGAAACGCTGACAGCg                                                              <  1:2437303/35‑1 (MQ=255)
                                           aaCTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTc                                     <  1:470655/60‑1 (MQ=255)
                                           aaCTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGt                          >  1:1531795/1‑71 (MQ=255)
                                            aCTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCta                                                   >  1:434898/1‑45 (MQ=255)
                                            aCTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGt                          >  1:966555/1‑70 (MQ=255)
                                             cTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCa                             >  1:758004/1‑66 (MQ=255)
                                             cTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGc                        <  1:3748249/71‑1 (MQ=255)
                                             cTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGc                        >  1:3632077/1‑71 (MQ=255)
                                             cTTTCTCGTGCAGACCTGGAAACGCTGACAGCGTATAACACATATACCATTACCGGTCTg                                   >  1:2555848/1‑60 (MQ=255)
                                                   cGTGCAGACCTGGAAACGCCGACAGCGTATAACACCta                                                   <  1:2020773/38‑1 (MQ=255)
                                                    gTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAg                   <  1:1614360/69‑1 (MQ=255)
                                                         gACCTGGAAACGCCGACAGCGTATAACACCTATAcc                                               <  1:3471750/36‑1 (MQ=255)
                                                         gACCTGGAAACGCCGACAGCGTATAACACCTATACCATTAGCGGTCTGCCGCCAGGTGCGATAGCGACGcc            <  1:4379464/71‑1 (MQ=255)
                                                            cTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAg                            >  1:545697/1‑52 (MQ=255)
                                                            cTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAg                            >  1:4094906/1‑52 (MQ=255)
                                                            cTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAg                            >  1:148166/1‑52 (MQ=255)
                                                            cTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCggg         >  1:1921048/1‑71 (MQ=255)
                                                            cTGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCggg         >  1:4021059/1‑71 (MQ=255)
                                                            cTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCgc          >  1:236956/1‑69 (MQ=255)
                                                             tGGAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCg                       <  1:491585/56‑1 (MQ=255)
                                                              ggAAACGCTGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCggggg       >  1:2494966/1‑71 (MQ=255)
                                                                aaaCGCCGACAGCGTATAACACCTATACCATTACCTGTCTGCCGCCAg                            >  1:2170813/1‑48 (MQ=255)
                                                                   cGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGAt  >  1:4561714/1‑71 (MQ=255)
                                                                    gCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCgggg        >  1:1422867/1‑64 (MQ=255)
                                                                      |                                                                   
GACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGAT  >  minE/734099‑734236

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