Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 316,095 A→G 64.2% N364S (AAC→AGC)  kefA → fused mechanosensitive channel proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE316,0950AG64.2% 129.1 / 118.4 123N364S (AAC→AGC) kefAfused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base G (38/41);  total (59/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 128 total aligned reads displayed.

GGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCT  >  minE/316025‑316163
                                                                      |                                                                    
 gCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGaaaa                                                                      >  1:1313001/1‑70 (MQ=255)
  cAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCa                                                                    >  1:3567170/1‑71 (MQ=255)
  cAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCa                                                                    <  1:3853223/71‑1 (MQ=255)
   aGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCAt                                                                   <  1:648954/71‑1 (MQ=255)
   aGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACAt                                                                   <  1:1225611/71‑1 (MQ=255)
   aGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACAt                                                                   >  1:76848/1‑71 (MQ=255)
         ctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCaa                                                             >  1:1487130/1‑71 (MQ=255)
         ctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCa                                                              <  1:1721958/70‑1 (MQ=255)
         ctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCaa                                                             >  1:4059004/1‑71 (MQ=255)
               tctcGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCa                                                        <  1:3217399/70‑1 (MQ=255)
                 tcGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCa                                                              <  1:3907527/62‑1 (MQ=255)
                 tcGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATc                                                      >  1:881836/1‑70 (MQ=255)
                   gTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGa                                                                 >  1:882087/1‑57 (MQ=255)
                   gTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATcgcg                                                   >  1:30408/1‑71 (MQ=255)
                   gTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATcgcg                                                   >  1:3423386/1‑71 (MQ=255)
                     aTCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGa                                                 >  1:2863600/1‑71 (MQ=255)
                      tCCTTTACCAGCAACAACATACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGAt                                                >  1:734104/1‑71 (MQ=255)
                      tCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGTAAAGCATGACCAACCGCATCGCGGa                                                 <  1:1534478/70‑1 (MQ=255)
                      tCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGAt                                                >  1:137502/1‑71 (MQ=255)
                       ccTTTACCAGCAACGACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCAt                                                                   <  1:1518134/51‑1 (MQ=39)
                       ccTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCAt                                                                   <  1:1819054/51‑1 (MQ=255)
                        cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGAttt                                              >  1:3057854/1‑71 (MQ=255)
                        cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGAttt                                              <  1:280977/71‑1 (MQ=255)
                        cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGAtt                                               >  1:2762042/1‑70 (MQ=255)
                          ttACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGc                                            >  1:1303080/1‑71 (MQ=255)
                          ttACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATcgc                                                    >  1:3844424/1‑63 (MQ=255)
                           tACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCgg                                                  <  1:3290451/64‑1 (MQ=255)
                           tACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGc                                            <  1:2541866/70‑1 (MQ=255)
                           tACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCg                                           >  1:2939878/1‑71 (MQ=255)
                            aCCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCa                                                                    >  1:565092/1‑45 (MQ=255)
                            aCCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCg                                           <  1:783370/70‑1 (MQ=255)
                             ccAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGtc                                         <  1:3382205/71‑1 (MQ=255)
                             ccAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGtc                                         <  1:4291750/71‑1 (MQ=255)
                               aGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACAt                                                                   <  1:4331295/43‑1 (MQ=255)
                                gCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAAc                                                                     <  1:3213820/40‑1 (MQ=255)
                                gCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCg                                      <  1:2292388/71‑1 (MQ=255)
                                 caacaacaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCaa                                                             >  1:4649372/1‑47 (MQ=255)
                                 caacaacaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCg                                      >  1:2944986/1‑70 (MQ=255)
                                   acaacaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGc                                                         <  1:3366840/49‑1 (MQ=255)
                                   acaacaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATc                                                      >  1:2089822/1‑52 (MQ=255)
                                    caacaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCg                                      <  1:281242/67‑1 (MQ=255)
                                    caacaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCg                                      <  1:579565/67‑1 (MQ=255)
                                    caacaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAAc                                   <  1:2682893/70‑1 (MQ=255)
                                    caacaaACGCTGCCCTCGGCGGATGAACTGGAAAACAt                                                                   <  1:1871919/38‑1 (MQ=255)
                                     aacaaACGCTTCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGa                                                 <  1:338840/55‑1 (MQ=255)
                                     aacaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGa                                                 <  1:3026114/55‑1 (MQ=255)
                                     aacaaACGCTGCCCTCGGCGGATGAACTGGAAAACAt                                                                   >  1:476461/1‑37 (MQ=255)
                                     aacaaACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACat                                 >  1:117850/1‑70 (MQ=255)
                                     aacaaACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAg                                 >  1:911349/1‑71 (MQ=255)
                                      acaaACGCTGCCCTCGGCGGATGAACTGGAAAGCATg                                                                  >  1:2040273/1‑37 (MQ=255)
                                       caaACGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCa                                                              <  1:2613610/40‑1 (MQ=255)
                                       caaACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGtt                               <  1:4497018/71‑1 (MQ=255)
                                        aaaCGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGt                                >  1:1607362/1‑69 (MQ=255)
                                         aaCGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTg                             <  1:3649958/71‑1 (MQ=255)
                                          aCGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGa                            <  1:2094419/71‑1 (MQ=255)
                                          aCGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGa                            <  1:4478502/71‑1 (MQ=255)
                                          aCGCTGCCCTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGa                            <  1:4011523/71‑1 (MQ=255)
                                                cccTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTaa                      <  1:3902306/71‑1 (MQ=255)
                                                 ccTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAc                     <  1:1366847/71‑1 (MQ=255)
                                                 ccTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAc                     <  1:4079431/71‑1 (MQ=255)
                                                  cTCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGa                                                 <  1:1807890/42‑1 (MQ=255)
                                                  cTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAcc                    <  1:3451217/71‑1 (MQ=255)
                                                   tCGGCGGATGAACTGGAAAGCATGACCAACCGCATCGCGGAtt                                               <  1:2415291/43‑1 (MQ=255)
                                                   tCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAg                                 <  1:2181028/57‑1 (MQ=255)
                                                   tCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAg                                 <  1:519648/57‑1 (MQ=255)
                                                   tCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTg                             <  1:2760275/61‑1 (MQ=255)
                                                   tCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAc                     >  1:2029563/1‑69 (MQ=255)
                                                    cggcggATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACcag                  >  1:445645/1‑71 (MQ=255)
                                                     ggcggATGAACTGGAAAGCATGACCAACCGCATCGCGGAtt                                               >  1:3703121/1‑41 (MQ=255)
                                                     ggcggATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAc                     >  1:1354404/1‑67 (MQ=255)
                                                     ggcggATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACcagc                 <  1:4204648/71‑1 (MQ=255)
                                                     ggcggATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACcag                  >  1:2151783/1‑70 (MQ=255)
                                                     ggcggATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACcag                  >  1:1409215/1‑70 (MQ=255)
                                                      gcggATGAACTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACcag                  >  1:1209202/1‑69 (MQ=255)
                                                        ggATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCagcagc              <  1:3993147/71‑1 (MQ=255)
                                                          aTGAACTTGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGt            <  1:4558351/71‑1 (MQ=255)
                                                          aTGAACTTGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGt            <  1:1425328/71‑1 (MQ=255)
                                                           tGAACTTGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGaa                           >  1:518977/1‑55 (MQ=255)
                                                           tGAACTGGAAAGCATGACCAACCGCATCGCGGAtt                                               >  1:1790796/1‑35 (MQ=255)
                                                           tGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTg           <  1:53070/71‑1 (MQ=255)
                                                            gAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAcc                    >  1:2111505/1‑61 (MQ=255)
                                                             aaCTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGttt                              <  1:1185293/50‑1 (MQ=255)
                                                             aaCTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGAt         >  1:2348021/1‑71 (MQ=255)
                                                              aCTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGt            <  1:4035004/67‑1 (MQ=255)
                                                              aCTGGAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATg        <  1:2968208/71‑1 (MQ=255)
                                                              aCTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAc                     >  1:496120/1‑58 (MQ=255)
                                                              aCTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATg        <  1:4393004/71‑1 (MQ=255)
                                                              aCTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATg        <  1:2905319/71‑1 (MQ=255)
                                                               cTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATg        >  1:273161/1‑70 (MQ=255)
                                                                tGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGaa                           >  1:2526725/1‑50 (MQ=255)
                                                                tGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCa      >  1:539552/1‑71 (MQ=255)
                                                                 ggAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGtt                        <  1:4219195/52‑1 (MQ=255)
                                                                 ggAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGa          >  1:772507/1‑66 (MQ=255)
                                                                 ggAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCAc     >  1:3904842/1‑71 (MQ=255)
                                                                  gAAAGCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAGCAGCAGCGTGATGCAct    >  1:1805076/1‑71 (MQ=255)
                                                                  gAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGa                            >  1:3439494/1‑47 (MQ=255)
                                                                  gAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAAcc                    >  1:1675295/1‑55 (MQ=255)
                                                                  gAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCAct    <  1:3184427/71‑1 (MQ=255)
                                                                  gAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCAct    >  1:4093399/1‑71 (MQ=255)
                                                                     aaCATGACCAACCGCATCGCGGATTTGCGTCTCGAACAg                                 <  1:2404638/39‑1 (MQ=255)
                                                                      |                                                                    
GGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGATTTGCGTCTCGAACAGTTTGAAGTTAACCAGCAGCGTGATGCACTCT  >  minE/316025‑316163

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