Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,276 T→C 7.8% E326E (GAA→GAG) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,2760TC 7.8% 247.4 / 9.9 89E326E (GAA→GAG) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (43/39);  new base C (5/2);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

CAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCT  >  minE/609207‑609346
                                                                     |                                                                      
cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                                                        >  1:2670185/1‑70 (MQ=255)
  gggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCa                                                                     <  1:3512495/71‑1 (MQ=255)
   ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                                                                    >  1:3673025/1‑71 (MQ=255)
     tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATc                                                                   >  1:1815092/1‑70 (MQ=255)
     tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                                                  <  1:1732713/71‑1 (MQ=255)
         aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTa                                                                <  1:2181355/69‑1 (MQ=255)
         aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                                                              <  1:1443009/71‑1 (MQ=255)
            gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                           <  1:2675117/71‑1 (MQ=255)
               gAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                                                                    >  1:3282865/1‑59 (MQ=255)
               gAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                                                        <  1:3213131/71‑1 (MQ=255)
                  gCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCa                                                                     <  1:628682/55‑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:2931341/1‑71 (MQ=255)
                      aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                 >  1:3649890/1‑71 (MQ=255)
                      aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                 >  1:2517411/1‑71 (MQ=255)
                       tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                <  1:3259992/71‑1 (MQ=255)
                       tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                <  1:3539381/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)
                          caccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                                                  <  1:3651723/50‑1 (MQ=255)
                          caccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAg                                             >  1:3445086/1‑71 (MQ=255)
                           accGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTc                                                                      >  1:873084/1‑45 (MQ=255)
                           accGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTc                                                                      >  1:2205307/1‑45 (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:3062336/1‑71 (MQ=255)
                                              cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                         <  1:3907715/71‑1 (MQ=255)
                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                 <  1:3788985/46‑1 (MQ=255)
                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                 <  1:573522/46‑1 (MQ=255)
                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                          >  1:2126738/1‑69 (MQ=255)
                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                         >  1:2917603/1‑70 (MQ=255)
                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                         >  1:307915/1‑70 (MQ=255)
                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                         >  1:1925454/1‑70 (MQ=255)
                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                         >  1:538969/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:1762852/1‑71 (MQ=255)
                                                         gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt              >  1:3845756/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)
                                                                     |                                                                      
CAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCT  >  minE/609207‑609346

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