Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,277	T→G	7.9%	E326A (GAA→GCA) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,277	0	T	G	7.9%	222.5 / 13.5	89	E326A (GAA→GCA) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (43/39);  new base G (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 = 8.84e-01

GGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCT  >  minE/609209‑609346                                                                     |                                                                      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:2517411/1‑71 (MQ=255)                     aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                 >  1:2931341/1‑71 (MQ=255)                     aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                 >  1:3649890/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: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:307915/1‑70 (MQ=255)                                              aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                         >  1:2917603/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: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: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:632035/71‑1 (MQ=255)                                                              gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa        <  1:2784164/71‑1 (MQ=255)                                                               aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                         >  1:492928/1‑53 (MQ=255)                                                               aaCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                         >  1:2072885/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)                                                                     |                                                                      GGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCT  >  minE/609209‑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: ‑