Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,118	G→A	13.4%	N712N (AAC→AAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,118	0	G	A	13.4%	219.4 / 20.1	82	N712N (AAC→AAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (38/33);  new base A (2/9);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 4.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCA  >  minE/608049‑608188                                                                      |                                                                       gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                                                       <  1:3744660/71‑1 (MQ=255)    aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                                                    <  1:385533/71‑1 (MQ=255)    aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                                                    >  1:396873/1‑71 (MQ=255)     gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                                   <  1:803100/68‑1 (MQ=9)      aTACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                                   >  1:3380607/1‑70 (MQ=255)        acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                >  1:1149639/1‑71 (MQ=255)         cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                <  1:251882/70‑1 (MQ=255)            ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg                                                            <  1:1812077/71‑2 (MQ=255)            ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                             >  1:2612834/1‑70 (MQ=255)                  tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                      <  1:928150/71‑1 (MQ=255)                   cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                       >  1:254837/1‑69 (MQ=255)                    aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                    <  1:3241491/71‑1 (MQ=255)                    aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                    >  1:3155460/1‑71 (MQ=255)                    aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                    >  1:1010417/1‑71 (MQ=255)                     gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                                  >  1:2916952/1‑56 (MQ=255)                     gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                     <  1:1912623/69‑1 (MQ=255)                      cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                  >  1:40621/1‑71 (MQ=255)                      cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                  <  1:3583132/71‑1 (MQ=255)                      cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                  <  1:2538351/71‑1 (MQ=255)                        taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                                >  1:1934025/3‑70 (MQ=255)                        gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                                                >  1:3690568/1‑71 (MQ=255)                           gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                             <  1:1929460/71‑1 (MQ=255)                           gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                                             >  1:616967/1‑71 (MQ=12)                              atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGt                                                              <  1:1328372/51‑1 (MQ=255)                              atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                          >  1:564393/1‑71 (MQ=255)                               tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                                         <  1:2149654/71‑1 (MQ=255)                                ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                        >  1:3251582/1‑71 (MQ=255)                                ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                        >  1:855116/1‑71 (MQ=255)                                ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                        >  1:534188/1‑71 (MQ=255)                                 tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                        <  1:3278536/70‑1 (MQ=255)                                     cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                   <  1:868636/71‑1 (MQ=255)                                     cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                   <  1:1278827/71‑1 (MQ=255)                                      gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                                   <  1:558239/38‑1 (MQ=255)                                      gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat                                  <  1:1627687/71‑1 (MQ=255)                                        gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                 >  1:1553947/1‑70 (MQ=255)                                         ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                    <  1:237535/50‑1 (MQ=255)                                         ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                    <  1:1953432/50‑1 (MQ=255)                                         ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                               >  1:1659397/1‑71 (MQ=255)                                           ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                              >  1:1977977/1‑70 (MQ=255)                                            tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                            >  1:3689004/1‑71 (MQ=255)                                             ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTGAGGTGCTGACCACCTTCGATGGATGc                            <  1:3419300/70‑1 (MQ=255)                                             ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                            <  1:1213134/70‑1 (MQ=255)                                              cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAgg                                                   <  1:904983/46‑1 (MQ=255)                                              cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCtt                          >  1:3171065/1‑71 (MQ=255)                                               aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                      >  1:1568133/1‑58 (MQ=255)                                               aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                      >  1:1577473/1‑58 (MQ=255)                                               aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                         >  1:350886/1‑71 (MQ=255)                                                tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                        >  1:136118/1‑71 (MQ=255)                                                 cgcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                    >  1:1356558/1‑42 (MQ=255)                                                  gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                 >  1:3904503/1‑44 (MQ=255)                                                  gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                 >  1:2517574/1‑44 (MQ=255)                                                  gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                 >  1:1872647/1‑44 (MQ=255)                                                    gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                             >  1:3597212/1‑46 (MQ=255)                                                    gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                    >  1:2158204/1‑71 (MQ=255)                                                     tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                   <  1:2830291/55‑1 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                               >  1:3094586/1‑41 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                         <  1:37952/63‑1 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                  <  1:1361846/70‑1 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                  >  1:1551415/1‑70 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                  <  1:1775066/70‑1 (MQ=255)                                                       gAGCAGCAATTCACGGTTCATCACGTTAACGTTCAg                                                    >  1:1227772/1‑36 (MQ=255)                                                       cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                  <  1:3562883/69‑1 (MQ=255)                                                       cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                 <  1:1099635/70‑1 (MQ=255)                                                       cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                 <  1:1307701/70‑1 (MQ=255)                                                        agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                        <  1:1502184/47‑1 (MQ=255)                                                        agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaa                <  1:2845961/71‑1 (MQ=12)                                                        agcagcATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                             >  1:2861022/1‑42 (MQ=255)                                                         gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa               <  1:2209198/71‑2 (MQ=255)                                                          cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGt              <  1:410139/71‑1 (MQ=255)                                                            gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa            >  1:690240/1‑71 (MQ=255)                                                            gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa            >  1:686513/1‑71 (MQ=255)                                                              aTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc          <  1:3187811/71‑1 (MQ=255)                                                               tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                    <  1:3859320/60‑1 (MQ=255)                                                               tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca         >  1:653763/1‑71 (MQ=255)                                                                ttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat        >  1:1077683/1‑71 (MQ=255)                                                                 tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                              <  1:2052940/48‑1 (MQ=255)                                                                  cACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                    >  1:923323/1‑57 (MQ=255)                                                                   aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc      <  1:597418/70‑1 (MQ=255)                                                                   aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc      <  1:1859894/70‑1 (MQ=255)                                                                    cGGTTCATCACGTTAACGTTCAGGTGCTTACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat    <  1:3467065/71‑1 (MQ=255)                                                                    cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat    <  1:1082883/71‑1 (MQ=255)                                                                    cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat    <  1:103763/71‑1 (MQ=255)                                                                      gTTCATCACGTTAACGTTCAGGTTCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa  <  1:1865490/71‑1 (MQ=255)                                                                      gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat    >  1:2550594/1‑69 (MQ=255)                                                                      gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa  <  1:3007658/71‑1 (MQ=255)                                                                      gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa  <  1:743175/71‑1 (MQ=255)                                                                      gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa  <  1:1093349/71‑1 (MQ=255)                                                                      |                                                                       GCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCA  >  minE/608049‑608188

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑