Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,291 A→T 7.7% T321T (ACT→ACA)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,2910AT 7.7% 253.5 / 11.4 91T321T (ACT→ACA) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (40/44);  new base T (5/2);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCG  >  minE/609222‑609360
                                                                     |                                                                     
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:2931341/1‑71 (MQ=255)
       aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                               >  1:3649890/1‑71 (MQ=255)
        tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                              <  1:3539381/71‑1 (MQ=255)
        tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                              <  1:3259992/71‑1 (MQ=255)
         accaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                              >  1:909603/1‑70 (MQ=255)
          ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                                       >  1:2663843/1‑60 (MQ=255)
           caccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAg                                                           >  1:3445086/1‑71 (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:541514/1‑69 (MQ=255)
                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                                    >  1:1964512/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:348863/1‑63 (MQ=255)
                                         cGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                     >  1:3909273/1‑63 (MQ=255)
                                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                            >  1:3845756/1‑71 (MQ=255)
                                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                            >  1:1762852/1‑71 (MQ=255)
                                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtgg                             <  1:2784738/70‑1 (MQ=255)
                                           ccAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTc                           <  1:1947220/71‑1 (MQ=255)
                                           ccAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTc                           <  1:3130669/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:2072885/1‑53 (MQ=255)
                                                 aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                       >  1:492928/1‑53 (MQ=255)
                                                 aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                   <  1:2626022/57‑1 (MQ=255)
                                                 aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                   <  1:1532569/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)
                                                               ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt       <  1:3454476/71‑1 (MQ=255)
                                                                tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                   >  1:1953983/1‑42 (MQ=255)
                                                                 ccGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTg      <  1:53838/70‑1 (MQ=255)
                                                                  cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt    <  1:472621/71‑1 (MQ=255)
                                                                   ggAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGct           <  1:1354607/63‑1 (MQ=255)
                                                                    gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCg  <  1:2907411/71‑1 (MQ=255)
                                                                    gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCg  <  1:495117/71‑1 (MQ=255)
                                                                    gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCg  >  1:3039745/1‑71 (MQ=255)
                                                                     |                                                                     
GAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCG  >  minE/609222‑609360

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