Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,283 T→A 7.8% Y324F (TAC→TTC) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,2830TA 7.8% 251.8 / 11.3 90Y324F (TAC→TTC) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (42/41);  new base A (5/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

ACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTT  >  minE/609216‑609352
                                                                   |                                                                     
aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTa                                                                      <  1:2181355/69‑1 (MQ=255)
aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                                                                    <  1:1443009/71‑1 (MQ=255)
   gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                                 <  1:2675117/71‑1 (MQ=255)
      gAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                                                              <  1:3213131/71‑1 (MQ=255)
         gCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTAc                                                           >  1:672663/1‑71 (MQ=255)
          cccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACg                                                          <  1:1902754/71‑1 (MQ=255)
            cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCa                                                        <  1:3732657/71‑1 (MQ=255)
             aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                       >  1:2178600/1‑71 (MQ=255)
             aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                       >  1:2517411/1‑71 (MQ=255)
             aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                       >  1:3649890/1‑71 (MQ=255)
             aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                       >  1:2931341/1‑71 (MQ=255)
              tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                      <  1:3539381/71‑1 (MQ=255)
              tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                      <  1:3259992/71‑1 (MQ=255)
              tACCACCGATAGATTCGGTTGCCCAGATCGCGTCGCCAGAGAACAGTTCATCGTATTcc                                                                  >  1:3144426/1‑59 (MQ=255)
               accaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                      >  1:909603/1‑70 (MQ=255)
                ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                               >  1:2663843/1‑60 (MQ=255)
                 caccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAg                                                   >  1:3445086/1‑71 (MQ=255)
                   ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTa                                                                      >  1:2173449/1‑50 (MQ=255)
                   ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc                                                  >  1:3293635/1‑70 (MQ=255)
                   ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc                                                  >  1:2412562/1‑70 (MQ=255)
                   ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                                 <  1:3848922/71‑1 (MQ=255)
                     gATAGATTCGGTTGCCCAGATCGTGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                               <  1:3586135/71‑1 (MQ=255)
                      aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                              >  1:3537815/1‑71 (MQ=255)
                       tAGATTCGGTTGCCCAGATCGAGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                             <  1:1846673/71‑1 (MQ=255)
                        aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTa                                                            >  1:3160464/1‑55 (MQ=255)
                          aTTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                            <  1:3952719/69‑1 (MQ=255)
                           ttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATAc                                         >  1:1544739/1‑71 (MQ=255)
                           ttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATAc                                         >  1:684166/1‑71 (MQ=255)
                            tCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                               >  1:2291340/1‑48 (MQ=255)
                            tCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                               >  1:2695288/1‑64 (MQ=255)
                               gTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCAt                                           <  1:1053385/65‑1 (MQ=255)
                                ttGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGt                                    <  1:1267864/71‑1 (MQ=255)
                                 tGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATAc                                         >  1:1444793/1‑65 (MQ=255)
                                  gCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATTCGc                                       <  1:1298279/66‑1 (MQ=255)
                                    ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                  <  1:3358455/69‑1 (MQ=255)
                                    ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                <  1:1184321/71‑1 (MQ=255)
                                     cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               <  1:3907715/71‑1 (MQ=255)
                                     cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:3062336/1‑71 (MQ=255)
                                      aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                       <  1:573522/46‑1 (MQ=255)
                                      aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                       <  1:3788985/46‑1 (MQ=255)
                                      aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                >  1:2126738/1‑69 (MQ=255)
                                      aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:538969/1‑70 (MQ=255)
                                      aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:1925454/1‑70 (MQ=255)
                                      aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:307915/1‑70 (MQ=255)
                                      aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:2917603/1‑70 (MQ=255)
                                       gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGa                             >  1:2039893/1‑71 (MQ=255)
                                        aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAg                                                   <  1:1041417/48‑1 (MQ=255)
                                        aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:1932102/1‑68 (MQ=255)
                                        aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                            >  1:2522340/1‑71 (MQ=255)
                                         tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTGATGAcc                           <  1:2507/71‑1 (MQ=255)
                                          cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                       >  1:2111563/1‑42 (MQ=255)
                                          cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                            >  1:1964512/1‑69 (MQ=255)
                                          cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                            >  1:541514/1‑69 (MQ=255)
                                            ggTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCAt                                           <  1:2904565/52‑1 (MQ=255)
                                            ggTCGCCAGAGAACAGTTCATCGTATT‑CGGAGTACGCAGGAAGCGAACCa                                            >  1:2743476/1‑50 (MQ=255)
                                             gTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                            <  1:3201019/66‑1 (MQ=255)
                                              tCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCAt                                           >  1:522535/1‑50 (MQ=255)
                                              tCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtg                      >  1:1979002/1‑71 (MQ=255)
                                              tCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggt                       <  1:3696411/70‑1 (MQ=255)
                                               cGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                             >  1:3909273/1‑63 (MQ=255)
                                               cGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                             >  1:348863/1‑63 (MQ=255)
                                                gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                    >  1:1762852/1‑71 (MQ=255)
                                                gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                    >  1:3845756/1‑71 (MQ=255)
                                                gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtgg                     <  1:2784738/70‑1 (MQ=255)
                                                 ccAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTc                   <  1:3130669/71‑1 (MQ=255)
                                                 ccAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTc                   <  1:1947220/71‑1 (MQ=255)
                                                  cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                  <  1:3357569/55‑1 (MQ=255)
                                                  cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                  <  1:721465/71‑1 (MQ=255)
                                                   agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTc                   >  1:1095852/1‑69 (MQ=255)
                                                   agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCaa                 >  1:3883954/1‑71 (MQ=255)
                                                    gagaACAGCGAATCAAATTCCGGTGTACGCAGGAAGCg                                                 <  1:732359/38‑1 (MQ=12)
                                                      gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa              <  1:2784164/71‑1 (MQ=255)
                                                      gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa              <  1:632035/71‑1 (MQ=255)
                                                       aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:492928/1‑53 (MQ=255)
                                                       aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:2072885/1‑53 (MQ=255)
                                                       aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                           <  1:1532569/57‑1 (MQ=255)
                                                       aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                           <  1:2626022/57‑1 (MQ=255)
                                                       aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAAc                >  1:1105144/1‑68 (MQ=255)
                                                       aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTCCATGAcc                           <  1:1927333/57‑1 (MQ=255)
                                                        aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                                 <  1:3705663/50‑1 (MQ=255)
                                                            ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCt        <  1:1304144/71‑1 (MQ=255)
                                                              cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa              >  1:1928517/1‑63 (MQ=255)
                                                               aTCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgcg     <  1:1382219/71‑1 (MQ=255)
                                                                tCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggt                       >  1:1712425/1‑52 (MQ=255)
                                                                 cGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTcc         >  1:2304435/1‑65 (MQ=255)
                                                                  gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                  <  1:1573467/39‑1 (MQ=255)
                                                                  gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAg                         <  1:662108/48‑1 (MQ=255)
                                                                  gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAt             >  1:676049/1‑60 (MQ=255)
                                                                  gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt  <  1:1959041/71‑1 (MQ=255)
                                                                  gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGct   <  1:2196871/70‑1 (MQ=255)
                                                                   |                                                                     
ACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTT  >  minE/609216‑609352

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: ‑