Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,282	G→A	8.2%	Y324Y (TAC→TAT) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,282	0	G	A	8.2%	255.4 / 9.4	85	Y324Y (TAC→TAT) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (41/37);  new base A (5/2);  total (46/39)
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.92e-01

TACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTT  >  minE/609212‑609352                                                                       |                                                                       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)           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)                   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)                        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:573522/46‑1 (MQ=255)                                           aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                       <  1:3788985/46‑1 (MQ=255)                                           aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                >  1:2126738/1‑69 (MQ=255)                                           aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:2917603/1‑70 (MQ=255)                                           aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:1925454/1‑70 (MQ=255)                                           aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                               >  1:307915/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: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: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: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)                                                                       |                                                                       TACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTT  >  minE/609212‑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: ‑