Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,240 A→C 9.7% S338S (TCT→TCG) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,2400AC 9.7% 173.3 / 12.1 72S338S (TCT→TCG) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base C (2/5);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCA  >  minE/609171‑609310
                                                                     |                                                                      
cAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATa                                                                        <  1:4373535/70‑1 (MQ=255)
  gggTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATAACACCGATAGAt                                                                     <  1:4219930/71‑1 (MQ=255)
       ttCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGt                                                                >  1:146097/1‑71 (MQ=255)
          aGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                              <  1:1815458/70‑1 (MQ=255)
          aGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGc                                                             <  1:3223435/71‑1 (MQ=255)
              aaCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAg                                                         <  1:5120824/71‑1 (MQ=255)
               aCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGa                                                        <  1:1618778/71‑1 (MQ=255)
                cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                                            >  1:2858751/1‑66 (MQ=255)
                cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGAt                                                       <  1:1468218/71‑1 (MQ=255)
                 ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                                            >  1:4980736/1‑65 (MQ=255)
                    aGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                                    <  1:2476638/70‑1 (MQ=255)
                     ggagTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGt                                                  <  1:2964678/68‑1 (MQ=255)
                        gTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCgg                                                                 >  1:4654558/1‑53 (MQ=255)
                         tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTggg                                                   >  1:4588029/1‑66 (MQ=255)
                         tttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                               <  1:151642/70‑1 (MQ=255)
                          ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                               >  1:1515688/1‑69 (MQ=255)
                          ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                               >  1:2138134/1‑69 (MQ=255)
                           tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCg                                                     >  1:2000489/1‑62 (MQ=255)
                           tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                                              >  1:2712682/1‑69 (MQ=255)
                           tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                                              >  1:2474570/1‑69 (MQ=255)
                            ttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                            <  1:1269378/70‑3 (MQ=255)
                              ggTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCagag                                          >  1:1969865/1‑70 (MQ=255)
                               gtcaCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGaa                                        <  1:4104932/68‑1 (MQ=255)
                               gTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                                         <  1:2512797/70‑1 (MQ=255)
                                  aCCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                                     >  1:845085/1‑71 (MQ=255)
                                  aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGcc                                              <  1:765651/62‑1 (MQ=255)
                                  aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAg                                     <  1:3508966/71‑1 (MQ=255)
                                    cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGtt                                   >  1:4221212/1‑71 (MQ=255)
                                    cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                    >  1:2647042/1‑70 (MQ=255)
                                    cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                    >  1:3660695/1‑70 (MQ=255)
                                    cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                            >  1:2796722/1‑60 (MQ=255)
                                       ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                                <  1:1293650/71‑1 (MQ=255)
                                         tACGACCGTCGAGGCCCATACCACCGATAGATTCgg                                                                 >  1:4678409/1‑36 (MQ=255)
                                         tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                              <  1:2079290/71‑1 (MQ=255)
                                         tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                              <  1:3889893/71‑1 (MQ=255)
                                           cGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                              <  1:3607041/37‑1 (MQ=255)
                                           cGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                              <  1:3613668/37‑1 (MQ=255)
                                             aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                          >  1:1698968/1‑71 (MQ=255)
                                             aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAt                           >  1:624818/1‑70 (MQ=255)
                                              ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                         >  1:572776/1‑71 (MQ=255)
                                              ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                         >  1:4124624/1‑71 (MQ=255)
                                                gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                                    <  1:2518763/42‑1 (MQ=255)
                                                gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                       <  1:850734/71‑1 (MQ=255)
                                                gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                       <  1:1691527/71‑1 (MQ=255)
                                                 tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                      <  1:4986670/55‑1 (MQ=255)
                                                 tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                      <  1:3899602/55‑1 (MQ=255)
                                                 tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                      <  1:3307660/55‑1 (MQ=255)
                                                 tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                      <  1:4247526/55‑1 (MQ=255)
                                                 tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                       >  1:4205076/1‑70 (MQ=255)
                                                    aGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                    >  1:4514107/1‑70 (MQ=255)
                                                        ccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTGCGGAGTACg                >  1:3484677/1‑70 (MQ=255)
                                                         cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                       <  1:3358088/46‑1 (MQ=255)
                                                         cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                       <  1:2868013/46‑1 (MQ=255)
                                                         cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                       <  1:4291754/46‑1 (MQ=255)
                                                          aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                       <  1:4092524/45‑1 (MQ=255)
                                                          aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg             >  1:4629968/1‑71 (MQ=255)
                                                           tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                       >  1:2804630/1‑44 (MQ=255)
                                                           tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                                     >  1:797551/1‑46 (MQ=255)
                                                           tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg            <  1:1712763/71‑1 (MQ=255)
                                                           tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg            <  1:3820930/71‑1 (MQ=255)
                                                             ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa           <  1:821195/70‑1 (MQ=255)
                                                             ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa           <  1:3361070/70‑1 (MQ=255)
                                                               accGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                       >  1:378396/1‑56 (MQ=255)
                                                                ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                         >  1:5616054/1‑53 (MQ=255)
                                                                 cGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                                         >  1:2719212/1‑36 (MQ=255)
                                                                 cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGaa                                        >  1:4196245/1‑37 (MQ=255)
                                                                 cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg       <  1:1657633/70‑1 (MQ=255)
                                                                 cGATAGATTCGGTTGCCCAGATCGGGTAGCCAGAGaa                                        >  1:2396641/1‑37 (MQ=255)
                                                                  gATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATc                               >  1:503729/1‑45 (MQ=255)
                                                                   aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTc                                  <  1:5490046/41‑1 (MQ=255)
                                                                   aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc    >  1:4952182/1‑71 (MQ=255)
                                                                    tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                         <  1:4114301/49‑1 (MQ=255)
                                                                    tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa      <  1:2988351/68‑1 (MQ=255)
                                                                     aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc        >  1:1460356/1‑65 (MQ=255)
                                                                     aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc   >  1:5282517/1‑70 (MQ=255)
                                                                     aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc    >  1:4756260/1‑69 (MQ=255)
                                                                     aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa  >  1:399609/1‑71 (MQ=255)
                                                                     |                                                                      
CAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCA  >  minE/609171‑609310

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: