Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 318,020 C→A 18.3% R1006S (CGT→AGT)  kefA → fused mechanosensitive channel proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE318,0200CA18.3% 202.4 / 28.3 98R1006S (CGT→AGT) kefAfused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base C (53/27);  new base A (13/5);  total (66/32)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

GTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGA  >  minE/317959‑318083
                                                             |                                                               
gTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTTATCCGTCTCgg                                                          <  1:3408611/69‑1 (MQ=255)
gTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCgg                                                          <  1:2824269/69‑1 (MQ=255)
gTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCgg                                                          <  1:30068/69‑1 (MQ=255)
   tGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCg                                                        >  1:2096330/1‑68 (MQ=255)
   tGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCg                                                        >  1:246237/1‑68 (MQ=255)
   tGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGt                                                       >  1:2711941/1‑69 (MQ=255)
   tGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGt                                                       >  1:2952089/1‑69 (MQ=255)
   tGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGt                                                       >  1:3385996/1‑69 (MQ=255)
       aCCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGcc                                                   <  1:2805123/69‑1 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGtctc                                                            >  1:1052792/1‑50 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTcc                                          >  1:1187303/1‑68 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTcc                                          >  1:79314/1‑68 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:3285628/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:440128/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:470243/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:2797753/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:860661/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:2422422/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:2399784/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:983103/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:1870321/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:105508/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:1058494/1‑69 (MQ=255)
                 tGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCg                                         >  1:1465617/1‑69 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGCTc                                      >  1:759444/1‑68 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATc                                      >  1:1488045/1‑68 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATc                                      >  1:766963/1‑68 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATc                                      >  1:1705665/1‑68 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATc                                      >  1:2419653/1‑68 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATc                                      >  1:1099054/1‑68 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATc                                      >  1:3011911/1‑68 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATc                                      >  1:1438541/1‑68 (MQ=255)
                     cAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATc                                      >  1:2239972/1‑68 (MQ=255)
                             cGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAg                             <  1:2237404/69‑1 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAg                             <  1:2714390/68‑1 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:2906074/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:1664306/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:3141189/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:1293916/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:2813422/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:160606/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:3323189/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:3393341/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:1208380/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:369362/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:386917/1‑69 (MQ=255)
                              gTTGACTGACACTACTACGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                            >  1:533284/1‑69 (MQ=255)
                                    tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:3305630/68‑1 (MQ=255)
                                    tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:1875517/68‑1 (MQ=255)
                                    tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:3371527/68‑1 (MQ=255)
                                    tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:1325123/68‑1 (MQ=255)
                                    tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:28166/68‑1 (MQ=255)
                                    tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:1014403/68‑1 (MQ=255)
                                    tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:1562134/68‑1 (MQ=255)
                                    tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:2639433/68‑1 (MQ=255)
                                         ctactaCGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAgtgt                 <  1:2775049/69‑1 (MQ=255)
                                         ctactaCGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAgtgt                 <  1:2372144/69‑1 (MQ=255)
                                         ctactaCGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAgtgt                 <  1:2563559/69‑1 (MQ=255)
                                         ctactaCGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAgtgt                 <  1:3255049/69‑1 (MQ=255)
                                          tactacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                <  1:393057/69‑1 (MQ=255)
                                          tactacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                <  1:2406459/69‑1 (MQ=255)
                                          tactacGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                <  1:2103086/69‑1 (MQ=255)
                                          tactacGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                <  1:453513/69‑1 (MQ=255)
                                          tactacGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                <  1:381387/69‑1 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaa                                >  1:3101522/1‑52 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                >  1:2813329/1‑68 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:2000654/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:2089856/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:537574/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:3095041/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:1847519/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:972599/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:804095/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:2845106/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:93945/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:2257120/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTa               >  1:2190660/1‑69 (MQ=255)
                                           actacGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                <  1:815244/68‑1 (MQ=255)
                                             tacGCGTCTGGTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                <  1:2254366/66‑1 (MQ=255)
                                                gcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:3218217/56‑1 (MQ=255)
                                                       gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGgcggcg   >  1:2627167/1‑69 (MQ=255)
                                                       gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGgcggcg   >  1:982414/1‑69 (MQ=255)
                                                       gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGgcggcg   >  1:2943948/1‑69 (MQ=255)
                                                       gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGgcggcg   >  1:2614423/1‑69 (MQ=255)
                                                       gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGgcggcg   >  1:3150697/1‑69 (MQ=255)
                                                       gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGgcggcg   >  1:2608060/1‑69 (MQ=255)
                                                       gTGATCAGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGgcggcg   >  1:2616028/1‑69 (MQ=255)
                                                         gATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAGgtgt                 <  1:74115/53‑1 (MQ=255)
                                                         gATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAgtgt                 <  1:187903/53‑1 (MQ=255)
                                                         gATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAgtgt                 <  1:1546825/53‑1 (MQ=255)
                                                         gATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAgtgt                 <  1:1900032/53‑1 (MQ=255)
                                                         gATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAgtgt                 <  1:1623182/53‑1 (MQ=255)
                                                          aTCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGtt                >  1:2259608/1‑53 (MQ=255)
                                                          aTCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGGGa  >  1:3077762/1‑67 (MQ=255)
                                                          aTCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGa  >  1:1064626/1‑67 (MQ=255)
                                                          aTCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGa  >  1:1285305/1‑67 (MQ=255)
                                                          aTCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGa  >  1:3097408/1‑67 (MQ=255)
                                                           tCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                       <  1:2128075/45‑1 (MQ=255)
                                                             |                                                               
GTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGA  >  minE/317959‑318083

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: