Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,602,585 T→C 22.1% intergenic (+121/+119) yfeW → / ← yfeX predicted periplasmic esterase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,602,5850TC22.1% 241.1 / 67.4 131intergenic (+121/+119)yfeW/yfeXpredicted periplasmic esterase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (46/56);  new base C (18/11);  total (64/67)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

Only 100 of 138 total aligned reads displayed.

CGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGG  >  minE/1602515‑1602654
                                                                      |                                                                     
 gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGAttt                                                                      <  1:1279961/71‑1 (MQ=255)
 gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGAttt                                                                      <  1:2185713/71‑1 (MQ=255)
   aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCg                                                                    <  1:2659121/71‑1 (MQ=255)
   aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCg                                                                    <  1:3419870/71‑1 (MQ=255)
   aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCg                                                                    <  1:1513567/71‑1 (MQ=255)
    aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGAttt                                                                      >  1:867772/1‑68 (MQ=255)
      aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGAtt                                                                       >  1:1973801/1‑65 (MQ=255)
       cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAgg                                                                >  1:5153268/1‑71 (MQ=255)
          ccTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGcc                                                              <  1:641564/70‑1 (MQ=255)
          ccTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCg                                                             <  1:2982001/71‑1 (MQ=255)
              aTGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGAt                                                          <  1:2846181/70‑1 (MQ=255)
              aTGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGAt                                                          <  1:614292/70‑1 (MQ=255)
              aTGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATa                                                         <  1:5031299/71‑1 (MQ=255)
               tGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATa                                                         <  1:1430137/70‑1 (MQ=255)
                  gtggcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGc                                                     >  1:3298991‑M1/53‑71 (MQ=255)
                   tgtcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCg                                                    >  1:1139451‑M1/52‑71 (MQ=255)
                   tggcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCg                                                    >  1:2012125‑M1/52‑71 (MQ=255)
                     gcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCg                                                                    >  1:2927651‑M1/50‑53 (MQ=255)
                     gcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTa                                                                  <  1:3622469‑M1/6‑1 (MQ=255)
                     gcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGtt                                                  <  1:3459302‑M1/22‑1 (MQ=255)
                     acaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATa                                                         >  1:2427606‑M1/50‑64 (MQ=255)
                     aCGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGc                                                               >  1:616019/1‑58 (MQ=255)
                      cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGc                                                               <  1:2024242/57‑1 (MQ=255)
                      cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGc                                                               <  1:4181072/57‑1 (MQ=255)
                      cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                  <  1:3969903/70‑1 (MQ=255)
                       aaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCa                                                <  1:4205089‑M1/24‑1 (MQ=255)
                       aaatgcctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGgcgtcca                                                <  1:4713640‑M1/24‑4 (MQ=255)
                        cTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGcc                                                              >  1:4264186/1‑56 (MQ=255)
                        cTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCg                                                             >  1:3117800/1‑57 (MQ=255)
                         atgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCTTTCAcg                                              >  1:9523‑M1/46‑71 (MQ=255)
                         atgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCAc                                               >  1:2198300‑M1/46‑70 (MQ=255)
                           gcctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATaa                                                        <  1:545159‑M1/16‑1 (MQ=255)
                            tCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccg                                           >  1:4105072/1‑71 (MQ=255)
                            cctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCAcgccg                                           >  1:4313221‑M1/43‑71 (MQ=255)
                            cctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCAcgccg                                           >  1:5410137‑M1/43‑71 (MQ=255)
                             cAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccg                                           >  1:3341650/1‑70 (MQ=255)
                             cAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccg                                           >  1:4164196/1‑70 (MQ=255)
                              tgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTa                                                                  <  1:921575‑M1/6‑1 (MQ=255)
                              ctgatgcgctacgcttatcaggcctacaattctcgcgccaCTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         >  1:4208134‑M1/41‑71 (MQ=255)
                              aGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccgc                                          <  1:5149331/70‑1 (MQ=255)
                               gatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                        <  1:3601042‑M1/32‑1 (MQ=255)
                               gatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                        <  1:3279022‑M1/32‑1 (MQ=255)
                                atgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                       >  1:1533701‑M1/39‑71 (MQ=255)
                                atgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCa                                         >  1:884991‑M1/39‑69 (MQ=255)
                                 atgcgctacgcttatcaggcctacaattctcgcgccaCTCGTa                                                                  <  1:2265395‑M1/6‑1 (MQ=255)
                                  gcgctacgcttatcaggcctacagattcttgcgccaTTCGTa                                                                  <  1:3864489‑M1/6‑1 (MQ=255)
                                  gcgctacgcttatcaggcctacagattcttgcgccaTTCGTa                                                                  <  1:1525568‑M1/6‑1 (MQ=255)
                                  gcgctacgcttatcaggcctacagattcttgcgccaTTCGTa                                                                  <  1:2108340‑M1/6‑1 (MQ=255)
                                  gcgctacgcttatcaggcctacagattcttgcgccaTTCGTa                                                                  <  1:535221‑M1/6‑1 (MQ=255)
                                  gcgctacgcttatcaggcctacagattctcgcgccaCTCg                                                                    >  1:1097183‑M1/37‑40 (MQ=255)
                                  gcgctacgcttatcaggcctacagattctcgcgccaCTCg                                                                    >  1:4861843‑M1/37‑40 (MQ=255)
                                  gcgctacgcttatcaggcctacagattctcgcgccaCTCg                                                                    >  1:1371043‑M1/37‑40 (MQ=255)
                                  gcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCg                                     <  1:534969‑M1/35‑1 (MQ=255)
                                   tACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                  <  1:2807604/57‑1 (MQ=255)
                                   tACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                    >  1:4429880/1‑71 (MQ=255)
                                    cgctacgcttatcaggcctacaattctcgcgccaCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                    <  1:3979161‑M1/36‑1 (MQ=255)
                                     gctacgcttatcaggcctacaattctcgcgccaCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                   >  1:5620472‑M1/34‑70 (MQ=255)
                                     ctacgcttatcaggcctacagattcttgcgccaTTc                                                                     >  1:2081101‑M1/34‑36 (MQ=255)
                                     ctacgcttatcaggcctacagattcttgcgccaTTc                                                                     >  1:2413192‑M1/34‑36 (MQ=255)
                                     cGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                >  1:2535100/1‑57 (MQ=255)
                                     cGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                    >  1:4481191/1‑69 (MQ=255)
                                      tacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCg                                                             >  1:5606946‑M1/33‑43 (MQ=255)
                                       tGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtg                                <  1:536118/71‑1 (MQ=255)
                                        cgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGgca                                  <  1:73329‑M1/38‑2 (MQ=255)
                                         gAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtgt                               >  1:2940791/1‑70 (MQ=255)
                                            tatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAcat                           <  1:768055‑M1/45‑2 (MQ=255)
                                            tatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATcc                                      >  1:1181507‑M1/27‑60 (MQ=255)
                                            tatcaggcctacaaattcccgcgccaTCCGTAGGCCGGATAAGGCGTTCAcgccgc                                          >  1:2548375‑M1/27‑56 (MQ=255)
                                             atcaggcctacagattcttgcgccaTTCGTAGGCCgg                                                            <  1:4443203‑M1/12‑1 (MQ=255)
                                              tcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAacatg                          >  1:769987‑M1/25‑68 (MQ=255)
                                              tcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCg                                     <  1:4204713‑M1/35‑1 (MQ=255)
                                              tcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCg                                     <  1:2144493‑M1/35‑1 (MQ=255)
                                              tcaggcctacaaattcccgcaccaTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                   >  1:4522534‑M1/25‑61 (MQ=255)
                                              cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCg                                                    <  1:2953306/44‑1 (MQ=255)
                                              cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCg                                                    <  1:3730246/44‑1 (MQ=255)
                                              cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCg                                                    <  1:5560262/44‑1 (MQ=255)
                                                caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCg                                                    <  1:4705983/42‑1 (MQ=38)
                                                caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCg                                                    <  1:2927814/42‑1 (MQ=38)
                                                caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                  >  1:5003444/1‑44 (MQ=38)
                                                caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                  >  1:1857770/1‑44 (MQ=38)
                                                caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                  >  1:841473/1‑44 (MQ=38)
                                                  atatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTc                                                 >  1:2608166/1‑43 (MQ=255)
                                                   tatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTc                                                 >  1:2734752/1‑42 (MQ=255)
                                                   tatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCAc                    >  1:5068593/1‑71 (MQ=255)
                                                     tatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                >  1:1222205/1‑41 (MQ=255)
                                                     tatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  1:1007690/48‑1 (MQ=255)
                                                     tacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATc                                       >  1:2597924‑M1/18‑50 (MQ=255)
                                                      atTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgc                                             <  1:195599/43‑1 (MQ=255)
                                                       ttGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                     >  1:801883/1‑50 (MQ=255)
                                                       ttGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                     >  1:280150/1‑50 (MQ=255)
                                                       cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCa                                                <  1:3717949‑M1/24‑1 (MQ=255)
                                                       cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  1:4790849‑M1/31‑1 (MQ=255)
                                                       cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  1:315120‑M1/31‑1 (MQ=255)
                                                       cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  1:3579973‑M1/31‑1 (MQ=255)
                                                         gattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCAcgccg                                           <  1:870899‑M1/29‑1 (MQ=255)
                                                         gAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                >  1:1994417/1‑37 (MQ=255)
                                                           ttctcacgccaTCCGTAGGCCGGATAAGGCGTTCAcgcc                                            <  1:2525628‑M1/28‑1 (MQ=255)
                                                              ttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  1:2060699‑M1/31‑1 (MQ=255)
                                                                 aCGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCt      >  1:2755046/1‑71 (MQ=255)
                                                                     tttCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAgg  <  1:2886075/71‑1 (MQ=255)
                                                                      |                                                                     
CGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGG  >  minE/1602515‑1602654

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: